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A  non-contact,  non-destructive  method  is  described  that  allows  the  user  to 
measure  pavement  texture  and  loaded  and  unloaded  profiles  from  a  moving  vehicle. 

The  device  uses  laser  distance-measuring  gages.  The  gages  are  arranged  such 
i’,uc’  that  their  readings  can  be  Interpreted  to  give  the  loaded  and  unloaded  pavement 
profiles  and  a  measure  of  its  texture.  Four  gages  are  mounted  on  a  rigid  beam 
that,  in  turn,  is  mounted  on  the  side  of  a  load  vehicle.  The  gages  read  the  distance 
to  the  pavement.  Three  of  the  gages  are  located  outside  the  range  of  influence  of 
the  load  wheel,  and  hence  the  pavement  below  these  gages  is  undeflected.  One  of  the 
gages  is  adjacent  to  Che  load  wheel  to  measure  the  Induced  deflection.  An  algorithm 
Is  presented  that  relates  all  measurements  to  a  common  datum. 

The  three  gages  yield  readings  that  are  Interpreted  to  give  the  profile  (do  the 
surveying).  The  gages  make  so  many  readings  and  so  quickly,  that  the  standard  devia¬ 
tion  of  the  readings,  taken  over  a  short  time  period,  gives  a  measure  or  the  texture. 

A  prototype  device  was  constructed  and  tested.  The  results  were  - 

1.  Pavement  profiles  were  calculated  using  the  new  system.  These  compare 
favorably,  though  not  completely,  with  the  transit  surveys. 

2.  Standard  deviations  were  recorded.  These  readings  were  taken  over¬ 
concrete,  asphalt  and  porous  friction  surface  pavements.  Differences 
were  noted  between  the  textures  of  concrete  (smooth)  and  asphalt 
(rough)  pavements.  No  correlation  between  this  measure  and  standard 
measures  was  made. 

3.  The  deflections  caused  by  prototype  loads  were  measured  and  recorded. 

This  was  done  from  the  load  vehicle,  while  the  vehicle  was  moving  across 
a  pavement  at  the  0.  S.  Army  Corps  Of  Engineers  Waterways  Experiment 
Station,  in  Vicksburg,  Mississippi.  The  scheme  for  obtaining  the  re¬ 
sulting  deflections  is  described  In  detail.  An  algorithm  for  accomodat¬ 
ing  the  misalignment  of  the  laser  gages  on  the  rigid  beam  Is  also 
presented. 

Suggestions  are  offered  for  improving  the  system. 
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I.  INTRODUCTION 


A.  Description  of  the  Problem 

This  study  addresses  the  problem  of  improving  the  abil¬ 
ity  of  airport  pavement  engineers  to  plan  the  maintainance 
of  asphalt  runways  and  taxiways.  Because  of  the  high  cost 
of  closing  a  runway  of  an  airport)  the  scheduling  of  mainte¬ 
nance  is  very  important.  This  problem  is  accentuated  by  the 
increasing  frequency  of  air  travel  and  of  the  weight  of  aii — 
craft.  Consequently,  any  system  that  will  allow  maintenance 
to  be  predicted  and  performed  before  a  failure  takes  place 
(requiring  an  unscheduled  closing)  would  be  most  welcome. 

At  present,  there  exist  several  methods  of  testing 
asphalt  pavements  and  of  predicting  when  they  will  require 
maintenance.  They  are  discussed  in  the  next  section,  the 
"Review  of  Literature".  They  all  suffer  from  one  or  more  of 
the  following  problems; 

1.  Being  very  time  consuming  (thus  expensive). 

?.  Only  being  practical  to  test  a  very  few  locations  on 
the  pavement. 
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3)  Requiring  the  destruction  of  a  piece  of  the  pavement 
in  o^der  to  get  the  necessary  parameters  for  predic¬ 
tions 

B.  Requirements  of  the  So  1 u t i on 


The  solution  to  the  problem  described  above  is  to 
establish  a  fast,  efficient  method  of  predicting  pavement 
performance.  The  method  should  be  such  that  large  areas  of 
a  pavement  can  be  traversed  and  tested.  The  method  should  be 
quick  and  easy  to  implement  so  that  3  minimum  of  time  is 
required  and  normal  airport  operation  can  be  resumed  as  fast 
as  possible.  It  would  be  preferable  to  have  a  test  that  did 
not  require  rebuilding  a  piece  of  the  pavement  after  the 
test  had  been  run.  Not  only  is  this  time  consuming,  but  it 
creates  discontinuities.  These  discontinuities  may  limit 
where  this  type  of  test  can  be  run.  It  would  be  best  if  the 
test  could  be  run  day  or  night,  so  that  testing  could  be 
done  whenever  a  lag  in  traffic  occurs. 

The  use  of  actual  wheel  loads  is  preferred  over  simu¬ 
lated  loadings,  as  they  more  accurately  portray  real  condi¬ 
tions  and  the  state  of  loading  under  which  the  pavement  will 
have  to  perform.  Full-scale  loadings,  using  prototype  air¬ 
craft.  would  be  the  most  desirable  The  load  vehicle  used 
should  be  one  that  is  available  when  needed,  so  that  test¬ 


ing  may  be  repeated,  as  desired 
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To  maximize  the  utility  of  a  method,  the  equipment 
should  require  little  or  no  training  to  operate  so  that 
there  will  be  no  great  expense  incurred  in  training  or  hir¬ 
ing  personnel.  Regular  maintenance  crews  should  be  able  to 
run  the  test. 

Lastly,  the  cost  of  the  system  should  not  be  prohibi¬ 
tive. 

C.  Proposed  So  1  u t i on 


The  proposed  solution  is  shown  schematically  in  Fi  e 
1-1.  Shown  is  a  load  vehicle  (here  represented  as  a  truck) 

carrying  the  system  that  will  analyze  the  pavement.  The 
system  consists  of  four  laser  distance  measuring  gages 
attached  to  a  rigid  beam  which,  in  turn,  is  attached  to  the 
side  of  the  load  vehicle.  One  gage  is  adjacent  to  the  load 
wheel,  which  is  used  to  measure  a  deflection  created  by  the 
load.  As  the  vehicle  traverses  the  pavement,  the  gages  meas¬ 
ure  the  distance  to  the  pavement  at  intervals  specified  by 
the  fifth  wheel 

In  the  process  of  measuring  deflections,  two  other 
important  pavement  characteristics  are  obtained:  the  unde¬ 
flected  profile  and  surface  texture.  The  texture  measure¬ 
ment  is  possible  because  the  gages  have  a  resolution  smaller 
than  the  asperities  of  the  asphalt  pa*eme-t 


The  resolution 


ec 


is  enhanced  by  tf?  £*ct  that  a  large  number  of  measurements 
can  be  obtained  in  the  near  vicinity  of  a  point  on  the  pave- 
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II.  REVIEW  OF  LITERATURE 


A  Introduction 


The  evaluation  of  asphalt  pavements  is  a  topic  of  much 
interest  and  discussion  among  highway  and  airport  engineers. 
The  need  to  assess  adequately  the  condition  of  pavements,  to 
predict  when  repair  will  be  necessary  and  to  carry  out  those 
repairs  in  an  economical  manner  is  a  pressing  one,  particu¬ 
larly  as  the  frequency  of  traffic  increases. 

This  review  of  literature  covers  methods  of  asphalt 
pavement  evaluation  that  have  developed  over  the  years. 

Par t i c u lar 1 y ,  objective  methods  are  detailed,  although  sub¬ 
jective  methods  exist  and  are  in  use.  The  objective  methods 
can  be  divided  into  two  categories  -  those  that  destroy  por¬ 
tions  of  the  pavement  in  the  process  of  evaluation,  and 
those  that  do  not.  Since  the  thrust  of  the  research  done 
for  this  report  deals  with  a  non-destructive  technique,  and 
because  it  is  felt  that  non-destructive  methods  will  see 


increasing  use  in  the  future,  destructive  techniques  are 
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Factors  that  affect  pavement  performance  are  considered 
in  this  review.  Par t i c u 1 ar 1 y ,  the  effect  of  weather  (sea¬ 
sonal  and  daily)  and  changes  in  material  parameters  due  to 
construction  or  pavement  history  are  examined. 

Of  the  several  parameters  sought  in  pavement  evalua¬ 
tion,  four  are  reviewed  in  this  Section  -  roughness,  sur¬ 
face  texture,  deflections  under  wheel  loads  and  modulus  of 
elasticity  of  the  subgrade.  Roughness  measures  changes  in 
tne  pavement  surface  elevation  that  can  affect  the  "ridea— 
bility"  of  the  pavement.  Texture  is  a  measure  of  the 
pavement's  ability  to  produce  friction.  Texture  is  fre¬ 
quently  used  to  evaluate  hydroplaning  potential.  Deflec¬ 
tion,  as  used  in  this  study,  refers  to  the  vertical  dis¬ 
placement  of  the  loaded  pavement  surface  from  its  unloaded 
position.  The  modulus  of  elasticity  of  the  subgrade  relates 
the  elastic  stress-strain  properties  of  the  subgrade.  This 
modulus  is  used  both  in  design  and  evaluation  of  pavements. 


Q.  Sub  i ec t i ve  and  Qb  i e c  t i ve  Methods  of  Pavement  Evaluation 

1.  Subjective  Methods 

The  most  common  subjective  measurement  of  pavement  ser¬ 
viceability  is  the  Present  Serviceability  Ratio  (PSR).  This 
method  described  by  Yoder  and  Witczak  (1975)  entails  driving 


over  the  pavement  in  question  and  rating  it  on  a  scale  of  0 
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to  5  (very  poor  to  very  good).  After  a  number  of  people 
have  done  this,  the  average  of  their  ratings  is  taken  and 
declared  to  be  the  PSR. 

Related  to  the  PSR  is  the  Present  Serviceability  Index 
(PS1).  The  PSI  is  a  number  derived  from  a  regression  equa¬ 
tion  relating  various  pavement  qualities  (roughness,  crack¬ 
ing.  area  of  patching,  etc.  ).  It  is  used  to  derive  results 
that  agree  with  the  PSR.  Evaluation  of  the  nature  and 
extent  of  the  parameters  of  the  equation  is  partially  sub¬ 
jective.  making  this  method  something  less  than  exact. 
Nevertheless,  both  the  PSI  and  PSR  play  useful  roles  in 
current  pavement  evaluation  (Yoder  and  Witczak  (1975)>. 

The  PSI  and  PSR  methods,  being  global,  have  the  disad¬ 
vantage  of  not  being  able  to  pinpoint  problem  areas.  This 
makes  repair  more  difficult.  Another  disadvantage  of  the 
methods  is  the  fact  the  that  they  reflect  variability  due  to 
human  judgment. 

No  subjective  methods  were  found  for  evaluating  the 
subg^ade.  texture  or  deflection  were  found. 

2.  Objective  Methods 

a.  Introduction.  Most  objective  methods  have  been  developed 
by  highway  engineers  to  obtain  measures  of  the  strength  of 
the  pavement  (as  a  indicator  of  its  remaining  life),  the 
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rideability  and/or  the  skid  resistance.  These  methods  usu¬ 
ally  involve  a  mechanical  or  electrical  device  that  make 
contact  with  the  pavement  and  make  measurements  that  reflect 
the  desired  quantity.  Moreover,  many  of  these  devices 
require  that  the  apparatus  be  stationary  with  respect  to  the 
pavement  at  the  time  of  measurement,  and  thus  yield  data  at 
only  one  location  per  set-up  of  the  apparatus. 

Deflection,  texture,  roughness  and  subgrade  modulus  are 
discussed  on  the  basis  of  whether  the  particular  test  is 
destructive  or  non-destructive,  contact  or  non-contact. 
Whether  or  not  the  test  is  continuous  or  discrete  follows 
from  the  description  of  the  device. 

b_.  Destructive  Tests  -  Contact  and  Non-contact.  Destructive 
tests  of  asphalt  pavements  have  been  practiced  for  some 
time.  One  of  the  most  popular  is  the  California  Bearing 
Ratio  (CBR).  This  test  entails  pushing  a  standard  cylinder 
into  the  base  (or  subbase  or  subgrade  of  a  pavement)  at  a 
prescribed  rate  and  measuring  the  resistance  required  to 
accomplish  this  task  This  resistance  is  then  correlated 
with  laboratory  tests  or  field  performance  data  {Baker 
(1975);  Bowles  (1970)>.  The  results  of  the  CBR  test  have 
been  correlated  with  the  modulus  of  elasticity  of  the 
subgrade  by  the  following  equation  {Asphalt  Institute 


(1973);  Yoder  and  Witczak  (1975>> 
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E  =  1500  C3R 

where  E  =  modulus  of  elasticity  (psi) 

CBR  =  California  Bearing  Ratio 

Voder  and  Witczak.  however,  note  that  "extreme  caution 
should  be  exerc ised. . . when  using  this  relation". 

The  modulus  of  elasticity  of  the  subgrade  can  also  be 
determined  from  the  laboratory  testing  of  a  sample.  The 
Asphalt  Institute  (1978)  describes  how  to  determine  the 
"resilient  modulus",  which  is  defined  as 


where 


KM  =  resilient  modulus 

*  deviator  stress  in  triaxial  cell 

e.  =  vertical  strain  of  sample 

V 

Others  have  described  slightly  more  subjective  ways  to 
determine  the  modulus.  The  Asphalt  Institute  (1973)  relates 
the  FAA  Soil  Classi f ication  to  the  modulus,  ranging  from  an 
F10  soil  with  a  5500  psi  modulus  to  an  Fa  soil  with  31.000 
psi  modulus.  Kezdi  (1975)  gives  approximate  ranges  of  the 
modulus  for  different  soils  and  conditions.  His  values 


range  from  50-400  psi  for  very  soft  clays  to  1 4» 000-28< 000 
psi  for  dense  sand  and  gravel. 

The  plate  bearing  test  is  another  common  destructive 
test  used  to  evaluate  the  soil  beneath  a  pavement.  It  con¬ 
sists  of  digging  a  pit  large  enough  to  accommodate  the  plate 
and  loading  the  plate  while  measuring  the  corresponding  set¬ 
tlement.  Usually  the  test  is  run  to  determine  the  modulus 
of  subgrade  reaction,  which  is  used  in  the  analysis  of  the 
pavement  system.  The  general  procedure  is  described  by 
Yoder  and  Witczak  (1975).  McLeod  (1957)  used  the  plate  load 
test  to  determine  a  relation  between  deflection  and  settle¬ 
ment. 


The  modulus  of  subgrade  reaction  is  defined  as 


k 


>7 

6 


where 


It  =  modulus  of  subgrade  reaction 


*7  =  st  res*  on  subor sde 


£  =  deflection  of  plate 


Often  this  test  is  run  before  the  pavement  is  constructed  so 

of  base. 


that  the  proper  thicknesses 
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subbases  and  surface  courses  can  be  computed. 

Vesic  and  Saxena  (1970)  studied  rigid  pavements  and  the 
effect  of  subgrade  reaction  on  the  AASHD  Road  Test  pave¬ 
ments.  They  found  that  it  was  very  difficult  to  find  a  sin¬ 
gle  value  of  the  modulus  that  would  predict  deflection, 
shear  stresses,  moments  and  contact  pressures  at  the  same 
time.  They  noted  that  the  modulus  was  a  function  of  how  the 
test  was  run  and  the  size  of  the  plate  used.  For  shallow 
depth  subgrades  (i.  e.  those  with  bedrock  near  the  surface), 
they  were  able  to  get  a  single  value  that  satisfied  all  the 
statical  parameters  (deflection,  moments,  etc. ). 

Terzaghi.  in  1955.  also  noted  that  the  modulus  of 
subgrade  reaction  varied  with  the  width  of  the  plate  (or 
footing)  resting  on  the  soil.  He  proposed  the  following 
corrections  for  the  modulus  of  subgrade  reaction  - 

on  clay 


M 


'?>< 


E 

B'> 


and  on  sand 


M"  =  M 


= 


■IV 


where  M  =  modulus  of  subgrade  reaction 


cr  =  stress  on  subgr  sde  caused  by  a  beam  of  width  E 


6  -  deflect  ’.on  of  subgrade  caused  by  a  beam  of 


width  B 

B  =  width  of  beam  on  subgrade 

B'  =  width  of  beam  to  be  used  in  design 

B''  =  width  of  beam  used  in  design 

M'  =  adjusted  modulus  of  subgrade  reaction 

M''  =  known  modulus  of  subgrade  reaction  for  one 
foot  wide  footing 


The  above  methods  were  developed  primarily  for  highway 
use.  While  they  are  common  tests  for  design,  they  are  not 
common  tests  for  pavement  evaluation  because  they  are  expen¬ 
sive,  time  consuming  and  interfere  with  traffic.  Other, 
more  conventional,  methods  of  testing  the  soil  below  the 
pavement  exist  (these  belong  to  the  province  of  soil  mechan¬ 
ics:  the  interested  reader  is  referred  to  Terzaghi  and  Peck 
(1967)),  Most  conventional  methods  of  soil  testing  require 
a  samp le  and  thus  are  destructive. 

As  far  as  could  be  determined  ,  there  are  no  destruc¬ 
tive  non-contact  methods  of  pavement  evaluation.  Thus,  this 
review  proceeds  to  examine  non-destructive  contact  and  non¬ 
contact  methods  of  deflection,  roughness,  texture  and 


sufcgrade  modulus  measurements 
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c..  Non-destructive  Tests  —  Contact  and  Non-contact. 

1.  Deflection.  The  deflection  of  the  surface  of  a  pave¬ 
ment  under  load  is  probably  the  most  obvious  and  frequent 
non—  destructive  method  of  evaluating  the  adequacy  of  a 
pavement.  Many  researchers  use  deflections  to  aid  in  the 
calculation  of  remaining  life  in  a  pavement.  Consequently, 
many  different  apparatuses  have  been  developed  to  measure 
pavement  deflections. 

Probably  the  most  widely  used  (and  widely  acclaimed) 
device  is  the  Benkelman  beam  {Yoder  and  Witczak  (1975)>. 

This  device,  which  sits  on  the  ground,  is  a  long  simple 
lever  with  a  dial  gage  placed  at  one  end.  The  long  side  of 
the  lever  is  usually  placed  between  the  dual  tires  of  a 
parked  highway  truck.  As  the  truck  moves  away,  the  rebound 
of  the  pavement  is  measured  by  the  movement  of  the  beam. 

The  number  of  studies  using  the  beam  are  too  numerous  to 
list.  Rather,  the  reader  is  referred  to  the  following 
representative  studies  -  Nichols  <1963),  Kondner  and  Krizek 
(1966).  Scrivner  and  Michalak  (1969).  Beca,  et  al.  (1974), 
and  Moore,  Hanson  and  Hall  (197S). 

The  advantages  of  the  Benkelman  device  (which  probably 
account  for  its  popularity)  are  its  ease  of  use,  simplicity 
of  construction  and  its  durability.  The  d i sadvantages  are 
the  length  of  time  it  takes  to  set  up.  the  fact  that  only 
one  deflection  measurement  is  made  at  each  site,  deflection 
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beneath  the  tire  cannot  b e  measured,  and  that  the  rebound  of 
ttie  pavement  is  measured  instead  of  the  deflection  caused  by 
a  moving  wheel  load.  This  last  qualification  is  important 
because  the  rebound  is  often  less  than  the  deflection  (at 
least  for  short  measurement  times)  since  the  pavement 
requires  time  to  return  to  its  original  position  (see  Led¬ 
better  (1976).  p92;  Harr  and  Ng-A-Gui  (1977),  p7B;  Boyer 

(1972),  p 136>. 

The  problem  of  making  only  one  measurement  per  set-up 
has  been  somewhat  circumvented  by  the  LaCroix  De f 1 ec tograp h 
•C(Beca,  et  al.  (1974)>.  This  device,  in  essence,  is  a 
true k-mounted  series  of  B en k e 1  man ’ b earns.  As  the  truck  moves 
forward,  one  beam  measures  the  deflection  of  a  wheel  while 
another  beam  is  positioned  in  front  of  it.  After  the  first 
beam  has  made  its  measurement,  it  is  picked  up  and  moved  in 
front  of  the  other  beam  and  makes  another  measurement,  and 
so  on.  This  permits  a  series  of  points  to  be  measured.  The 
State  of  California  has  developed  a  similar  device  called 
the  traveling  d  ef  1  ec  tog  r  ap  h  -CBeca,  et  al.  (1974);  Yoder  and 
witciak  (1975);  Moore,  Hanson  and  Hall  (1978);  and  Zube 
<196o)>.  This  device  travels  very  slowly  (about  1/2  mph) 
and  can  make  measurements  at  about  twenty  foot  intervals. 

Another  very  popular  method  of  making  deflection  meas¬ 
urements  uses  steady-state  vibrators.  These  devices  all 
operate  on  the  same  basic  principle:  a  plate  is  placed  on 

the  pavement  surface  and  is  c-xcited  vertically  with  an 
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eccentric  u/eight  vibrator  or  by  high  capacity 
hydraulic  means.  Measurements  of  the  surface 
waves  produced  are  then  made  by  transducers  placed  on  the 
plate/pavement.  Several  vibrators  are  available  commer¬ 
cially.  notably,  the  Dynaflect.  Shell  Vibrator  and  the  Road 
Rater  {(Beca.  et  al.  (1974);  Voder  and  Witczak  (1975). 

Baker  (1975);  Public  Worts  (1973)>.  Variables  that  affect 
the  parameters  obtained  from  vibrators  are  the  static  weight 
of  the  vibrator,  the  frequency  of  vibration,  the  diameter  of 
the  plate  and  the  load  induced  during  vibration  {(Green  and 
Hall  (1975)>.  The  size  of  the  vibrators  range  from  portable 
models  to  the  trailer  mounted  lo-kip  Waterways  Experiment 
Station  model. 

The  parameter  most  commonly  determined  by  vibrators  is 
the  elastic  modulus  of  the  subsurface  soils  .  This  is  done 
by  relating  the  wavelength  produced  by  the  vibrator  to  the 
velocity  of  wave  propagation.  This  is  then  related  to  the 
elastic  modulus.  Comparisons  of  laboratory  and  field  values 
show  that  it  is  difficult  to  get  agreement  between  the  two 
{(Moore.  Hanson  and  Hall  (1978)>. 

Anderson  (1976)  reports  that  the  State  of  Utah  used  the 
Dynaflect  as  part  of  a  regular  pavement  evaluation  program. 
Peterson  and  Shepherd  (1972)  used  the  deflection  basin 
caused  by  a  Dynaflect  to  determine  parameters  that  indicated 
which  layer  beneath  the  pavement  is  ac  fault  when  excessive 
deflections  are  found  He  also  notes  that  Dynaflect  deflec¬ 
tions  cannot  be  related  tc  3enkelman  beam  measurements. 


Poehl  and  Scnvner  (1971)  used  the  Dynaflect  to  i.'easure 
changes  in  deflections  with  the  seasons  This  topic  is 
covered  in  a  later  part  of  this  review  Lastly,  Vang 

(1977)  uses  a  vibrator  that  varies  the  force  and  frequency 
applied  to  the  pavement  with  their  frequency  sweep  vibrator 
( Eng ineer i no  News  Record  (1978);  Vang  (1977)>. 

Advantages  of  the  vibrators  are  that  they  are  (gen¬ 
erally)  easily  transported,  and  are  relatively  quick  to  use. 
They  require  somewhat  more  skill  than  does  the  Benkelman 
beam,  and  they  are  considerably  more  expensive.  Because  the 
vibrators  do  not  simulate  the  actual  pavement  loading  condi¬ 
tion  that  the  pavement  experiences  during  its  life,  the 
interpretation  of  the  test  results  is  more  complicated  than, 
say.  a  static  deflection  measurement.  Furthermore,  the 
deflection  basin  created  by  a  vibrator  is  not  the  same  as 
that  created  by  an  actual  wheel  load  {(Kennedy  (1978)>.  The 
difference  in  shape  and  magnitude  of  the  vibrator  deflection 
basin  adds  to  the  confusion  surrounding  the  interpretation 
of  results  from  this  test.  The  vibrators  also  have  the 
disadvantage  of  measuring  only  one  point  on  the  pavement  per 
set-up  (i.  e.  ,  it's  not  a  continuous  measur ement ) .  Conse¬ 
quently,  interference  with  traffic  is  often  unavoidable. 

It  is  interesting  to  iote  that  Ved~os  and  Barker  (1977) 
used  both  she  Benkelman  bear  and  a  vibrascr  (lie  Road  Rater; 
on  two  similar  sections  of 
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pavement  in  Kentucky,  but  mere  unable  to  determine  why  one 
section  failed  and  the  other  did  not 

Another  wave  propagation  technique  is  the  impulse  test 
-CBeca.  et  al.  (1974);  Moore.  Hanson  and  Hall  <1978)>.  This 
method  entails  the  dropping  of  a  weight  onto  the  pavement 
and  measuring  the  deflections  produced.  Measurements  ar3 
made  uiith  velocity  or  displacement  transducers.  The  pro¬ 
perty  sought  is  generally  the  modulus  of  the  subsoils. 

Bohn,  et  al.  (1972)  describe  testing  with  the  French  Fal¬ 
ling  Weight  Def 1 ec tome  ter .  They  report  finding  a  good 
correlation  between  falling  weight  deflections  and  moving 
wheel  deflections.  Classen.  Valkering  and  Dimirsch  (1976) 
report  good  results  with  their  falling  weight  device.  Their 
elastic  parameters  and  layer  thicknesses  agree  closely  with 
core  samples.  The  impulse  test  has  about  the  same  advan¬ 
tages  and  d i sadvantages  as  the  steady— state  vibrator. 

Measures  of  permanent  deflection  are  also  used  to 
evaluate  highway  pavements.  Par t ic u lar  1  y .  the  depth  of  the 

ruts  made  by  the  tires  is  used  as  a  measure  of  the  pavement 
condition.  The  procedure  is  to  lay  a  straight-edge  across 
the  rut  and  measure  the  distance  from  the  straight-edge  to 
the  deepest  part  of  the  rut.  This  device  is  sometimes 
called  a  curvature  meter  (Moore,  Hanson  and  Hall  (1973)>. 

The  va lue  of  this  measure  is  somewhat  questionable,  although 
it  is  a  good  measure  of  the  serviceability  {(Kondner  and 


Kriiek  (1 966 )>.  Huang  (1971)  measured  the  curvature  with  a 
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straight-edge  twice  the  radios  jf  the  tire  print  He  used 
this  value  to  derive  a  curvature  equation 

Displacement  transducers  have  been  used  to  measure 
deflection  under  and  near  actual  wheel  loads.  Moore,  Hanson 
and  Hall  (1978)  discuss  the  General  Electric  Travel  Gage,  a 
device  embedded  in  the  pavement,  anchored  at  depth,  that 
measured  how  the  surface  moved  It  was  introduced  in  1938 
but  was  never  used  widely.  More  recently,  Doyer  and  Harr 
(1972)  embedded  transducers  in  airport  pavements  in  order  to 
measure  the  surface  deflection  under  pr  near  moving  wheel 
loads.  Highter  and  Harr  (1975)  also  used  transducers  in 
airport  pavements  to  measure  deflections  and  the  shape  of 
the  deflection  basin.  They  attached  the  transducers  to  rods 
anchored  17.  5  feet  below  che  pavement  in  order  to  remove  the 
influence  of  the  wheel  loads.  Baladi  and  Harr  (1978)  used 
transducers  on  a  rigid  beam  to  measure  the  deflection  basin 
near  wheel  loads  This  device  had  the  advantage  of  being 
able  to  measure  the  deflection  basin  at  many  locations 
because  the  beam  was  easily  transported.  Ledbetter  (1976) 
used  displacement  transducers  in  airport  pavements  to  meas¬ 
ure  the  responses  to  maneuvering 

In-place  transducers  have  the  advantage  of  being  able 
to  measure  deflections  beneath  a  wheel  -  something  none  of 
the  other  devices  can  do.  However,  they  are  time  consuming 
to  install  and  cannot  be  moved  easily  unc?  installed.  Port¬ 
able  beams  with  transducers  on  them  can  maVe  measurements  at 
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different  locations  but  cannot  measure  beneath  the  wheel. 


There  exist  a  few  methods  of  measuring  deflections  with 
non-contact  devices.  All  involve  reflected  light. 

Baker  (1957)  measured  deflections  by  taking  pictures  of 
the  pavement  before  and  during  loading  An  accuracy  of 
0.0005  inches  was  achieved.  This  method,  although  accurate, 
is  slow  and.  in  the  near  future,  will  probably  be  reserved 
for  research  purposes. 

Still  and  Winnett  (1975)  used  lasers  and  charged  cou¬ 
pled  devices  (similar  in  operation  to  a  row  of  photocells) 
to  measure  deflection  of  pavements.  The  operation  of  the 
device  was  similar  to  that  used  bg  Harr  and  Ng-A-Qui 
described  below. 

Harr  and  Ng-A-Qui  (1977)  developed  a  non— contact  device 
using  light  emitting  diodes  and  a  linear  photocell.  The 
light  from  the  diode  shined  on  the  pavement  and  reflected  up 
onto  the  photocell.  As  the  position  of  the  pavement  changed 
(under  load),  the  position  of  the  light  spot  on  the  photo¬ 
cell  changed,  providing  a  measure  of  deflection.  A  series 
of  these  gages  was  mounted  on  a  rigid  portable  beam  and  was 
used  to  measure  deflections  adjacent  to  moving  wheel  loads. 

The  advantage  of  these  methods  lies  in  the  fact  that 


they  do  not  affect  the  pavement  when  making  the 


measurements  Transducers  on  a  portable  beam  afford  a  high 
degree  of  flexibility,  but  do  not  permit  measurements  to  be 
made  beneath  the  wheel.  Photographic  methods  also  suffer 
this  limitation.  Moreover,  present  photographic  methods  do 
not  allow  continuous  measurements. 

2.  Texture.  The  surface  texture  of  pavements  is  of 
interest  to  the  engineer  who  desires  to  evaluate  the  skid 
resistance  or  hydroplaning  potential  of  a  pavement. 

Surface  texture  of  asphalt  pavements  is  most  commonly 
measured  by  contact  methods.  Rose,  Hutchinson  and  Galloway 
(1973)  provide  an  interesting  review  of  methods  that  have 
been  used.  They  cite  the  following  methods  - 

1)  The  “patch"  method  -  Rub  a  given  quantity  of 
some  substance  onto  the  surface  of  the  road 
until  the  substance  forms  an  approx imately  cir¬ 
cular  area  whose  surface  is  at  the  height  of  the 
projections  of  the  asphalt  aggregate.  The  diam¬ 
eter  of  the  circle  is  a  measure  of  texture 
(large  diameter  implies  smooth  texture). 
Materials  used  include  sand,  grease  and  silicone 
putty. 

2)  The  direct  measurement  method  -  This  consists  of 


22 


drawing  a  feeler  needle  across  the  surface  and 
watching  its  variation  in  height  Large  varia¬ 
tion  in  needle  height  implies  rough  texture  {see 
also  Moore  (1966),  and  Goodman  (1970)>. 


3)  Mi  seel laneous  methods  - 

a)  Lay  a  piece  of  metal  foil  on  the  pavement, 
strike  it  with  a  standard  rubber  mallet  and 
count  the  number  of  holes  made  in  the  foil. 
This  method  has  the  advantage  of  finding  the 
sharp  points,  which  are  not  distinguished  by 
the  previous  methods. 


b)  Make  a  plaster  cast  of  the  surface,  smear 
it  with  paint  and  measure  the  percent  of 
area  that  took  paint. 

Moore  (1966)  developed  an  indirect  method  for  measuring 
pavement  texture.  His  device,  called  an  outflow  meter, 
works  thusly:  a  Lucite  cylinder  with  a  circular  rubber 
gasket  around  the  bottom  is  pressed  down  on  the  pavement  and 
the  cylinder  is  Filled  with  water.  The  time  required  for 
the  water  to  flow  out  of  the  cylinder  is  a  measure  of  the 
pavement  texture.  Rough  pavements  will  let  the  water  flow 
out  quickly,  while  smootn  ones  take  longer  because  a  better 


seal  is  made  between  the  gasket  and  the  pavement. 


There  are  some  non-cont jet  methods  of  evaluating  the 
texture  of  a  pavement  Except  for  one  {(Cooper  (1974>>/ 
they  are  strictly  research  devices  Basically,  all  depend 
on  the  principle  of  shining  light  on  the  pavement  and 
measuring  the  variation  ir  its  reflection. 

Rose.  Hutchinson  and  Galloway  (1973)  note  that  stereo¬ 
photographs  are  sometimes  used  to  measure  texture.  Keeping 
a  microscope  focused  on  a  sample  as  it  moves  slowly  past  the 
lens  is  another  method.  The  movement  of  the  lens'  barrel  up 
and  down  is  recorded,  providing  a  measure  of  texture. 

Goodman  (1970)  used  a  vertical  light  shining  on  the 
pavement  to  obtain  his  measure  of  texture  He  measured  the 
amount  of  light  scattered  to  both  sides  of  the  source  to 
determine  the  "mean  void  width".  He  also  measured  the  image 
cast  by  light  reflecting  off  the  pavement  to  get  a  measure 
of  the  drainage  depth.  The  measurements  correlated  well 
with  those  made  by  a  mechanical  stylus,  although  the  accu¬ 
racy  was  not  given.  His  attempt  to  mount  the  device  on  a 
moving  trailer  was  met  with  limited  success. 

Gee.  King  and  Hegmon  (1974)  described  a  preliminary 
study  they  did  using  lasers  to  measure  texture.  They  meas¬ 
ured  the  change  of  the  dimensions  of  an  ellipse  produced  on 
a  pavement  by  a  laser  as  a  measure  of  the  texture. 

Gee  (1973)  used  polarized  light  to  measure  texture.  The 
degree  to  which  the  light  was  depolarized  after  reflection 
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mas  his  measure  of  the  texture  This  system  mas  mounted  on 
a  vehicle  and  driven  over  test  roads  to  gather  data 


Cooper  (1974)  used  lasers  and  charged  coupled  devices 
(CCDs)  to  measure  texture.  He  shined  the  light  on  the  pave¬ 
ment  so  that  it  reflected  onto  the  CCD.  The  texture  mas 
determined  by  noting  the  changes  in  position  of  the  light  on 
the  CCD 


3.  Roughness.  Pavement  roughness  is  of  interest  to 
engineers  concerned  uiith  the  users  comfort  and  safety. 

Rough  pavements  make  vehicle  control  difficult  and  the  ride 
uncomfortable  At  present  there  are  tmo  basic  methods  of 
measuring  roughness  -  (1)  quantifying  the  change  in  the 
pavement  profile  (either  locally  or  globally)  and.  (2) 
measuring  vehicle  accelerations  caused  by  the  changes  in 
profile.  Of  the  tmo.  the  former  is  more  common. 

Yoder  and  Witczak  (1975)  provide  a  review  of  some  com¬ 
monly  used  contact  profilometers.  The  straight-edge  profi- 
lometer  is  the  most  common.  It  measures  the  distance  to  the 
pavement  from  a  straight-edge  supported  on  both  ends  by 
mheels.  The  accuracy  of  this  device  is  limited  by  the 
mheelbase  to  the  measurement  of  irregularities  shorter  than 
the  mheelbase.  The  "slope  prof i lometer"  is  an  improvement 
on  the  common  straight-edge  in  that  it  measures  differences 
in  surface  slope  frcm  a  fixed  horizontal  datum.  Its  accu¬ 
racy  is  still  limited  by  its  mheelbase  The  CHLOE 


profilometer  is  similar  to  the  slope  profilometer  except 
that  it  lacks  the  horizontal  reference  device.  The  GMR  pro¬ 
filometer  uses  a  truck  for  the  straight-edge.  An  extra 
uheel  is  attached  under  the  center  of  the  truck  and  its 
movements!  monitored  by  accelerometers  and  potentiometers, 
measure  the  roughness.  This  is  a  high  speed  device,  the 
advantage  of  which  is  obvious. 

The  roughometer  varies  in  concept  from  the  straight¬ 
edge  devices.  It  uses  a  tracking  wheel  mounted  in  a  heavy 
trailer  frame.  Presumably,  the  inertia  of  the  frame  pro¬ 
vides  the  measuring  device  with  a  datum  against  which  varia¬ 
tions  in  profile  are  determined.  This  is  also  a  relatively 
high  speed  device 

The  roadmeter  {(Yoder  et  al  <1973>>  is  different  from 
either  of  the  above  described  devices.  It  uses  a  cable 
attached  to  the  rear  axle  of  a  car  to  pull  on  a  sliding 
electrical  contact.  The  contact,  counterbalanced  by 
springs,  moves  in  proportion  to  the  amount  of  relative  dis¬ 
placement  between  the  axle  and  the  car  body.  The  statistic 
derived  is  a  function  of  the  weighted  sums  of  displacements 
(broken  into  1/8  in.  intervals).  That  is, 

road me  ter  statistic  =  CIA  +  4B  +  9C  +  1 6D  +  ....)/ 64 

where  A  =  no  of  car-axle  deviations  equal  to  1/8  in 
B  =  no  of  car -axle  deviations  equal  to  2/8  i n 
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C  =  no  of  car-axle  deviations  equal  to  3/8  in 
*  etc. 


Another  contact  method*  providing  roughness  measure¬ 
ments*  follows  from  a  physical  survey  of  a  pavement.  This 
method*  though  slow/  is  accurate  and  is  the  method  that  most 
other  methods  are  checked  against  It  provides  accurate 
locations  of  trouble  spots  and  data  that  are  easy  to  evalu¬ 
ate. 

Houbolt  (1961)  suggests  the  use  of  vertical  vehicle 
acceleration  as  a  measure  of  roughness.  The  greater  the 
accelerations*  the  more  the  passenger  discomfort  will  be. 

He  recommended  O.  3  g  as  the  maximum  allowable  acc e 1  era t i on. 
Hall  and  Kapelson  (1962)*  in  a  similar  study,  recommend  0.  5g 
as  their  maximum.  This  method*  while  sound  conceptually* 
fails  to  pinpoint  the  location  of  rough  spots.  Moreover* 
since  the  accelerations  airplanes  experience  depend  on  the 
history  of  accelerations  during  a  take-off  or  landing*  it  is 
difficult  to  evaluate  any  particular  section  of  a  runway. 

The  Highway  Research  Board's  Special  Report  No.  95 
(1968)  describes  the  philosophy  and  mechanics  of  contact 
roughness  measurement  and  evaluation.  Besides  contact 
methods  of  measuring  pavement  roughness,  there  are  also 
non-contact  methods 
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Dickerson  and  Mace  (1976)  describe  a  device  they  made 
that  uses  lasers  to  measure  distance  to  the  pavement  By 
arranging  three  of  these  on  a  rigid  beam  towed  behind  a  car# 
they  were  able  to  measure  the  profile  of  the  pavement  over 
long  distances.  Their  device  had  the  advantage  of  being 
able  to  measure  the  profile  with  respect  to  a  datum,  making 
absolute  variations  of  the  surface  measurable. 

Joyce  (1975)  used  an  acoustic  device  to  measure  the 
profile  of  the  pavement.  He  measured  the  change  in 
wavelength  of  a  signal  bounced  off  the  pavement  from  a  mov¬ 
ing  vehicle.  This  provided  a  measure  of  the  distance  from  a 
point  on  the  vehicle  to  the  pavement  surface. 

Both  methods  have  the  advantage  of  being  fast  and  con¬ 
tinuous.  The  data  they  provide  permit  the  location  of  trou¬ 
ble  spots.  They  suffer  the  disadvantages  of  being  expensive 
and  comp  1 icated. 

4.  Modulus  of  Elasticity  of  the  Subgrade.  In  Part  1  of 
this  subsection,  where  deflections  were  discussed#  it  was 
mentioned  that  these  deflections  could  be  used  to  evaluate 
the  modulus  of  elasticity  of  the  subgrade.  This  has  been 
done  using  vibrators#  such  as  the  Shell  and  Road  Rater 
models.  Cunny  and  Fry  (1973)  gave  a  brief  description  of 
how  the  vibrator  methodology  works. 


When  the  vibrator  is  placed  on  the  pavementi  it  sends 
waves  out  along  the  surface.  The  velocity  of  these  waves 
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can  be  obtained  from 


v  =  f  X 


where  v  =  velocity 


f  =  frequency  of  vibration 


X  =  wavelength 


This  is  approximately  equal  to  the  shear  wave  velocity,  V 

If  the  density,  p,  of  the  soil  is  known,  and  the  compres¬ 
sion  wave  ••velocity,  ,  is  known,  one  can  use 


o 

=  V  p 


=  U  -V 
r  c  ■■ 


>:  -  o  i 

_ r  _ 

2  <  -  1  > 


to  get  the  modulus 


E  =  '  2  v  »  £  :•<;  1  ri'  > 


In  the  above, 


G 


S  a  r  m  o  lulu  ; 
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V  =  Poissons''  ratio 


E  =  modulus  of  subgrade  reaction 


Cunny  and  Fry  note  that  the  properties  thus  derived 
apply  best  to  soils  at  depths  of  ~. 

Weiss  (1977),  working  with  the  Waterways  Experiment 
Station  vibrator,  proposed  two  methods  of  determining  the 
modulus.  The  first,  called  the  "dynamic  frequency  response 
spectrum  method",  models  the  soil  with  a  mass,  spring  and 
dashpot.  By  using  the  amplitudes  of  the  waves  produced  by 
the  vibrator,  the  properties  of  the  model  are  determined. 

It  is  an  iterative  process.  Once  obtained,  these  properties 
are  used  in  the  Chevron  layered  elastic  computer  program  to 
calculate  the  modulus.  The  moduli  of  the  pavement  and  base 
courses  also  must  be  known  before  the  Chevron  program  can  be 
run.  The  results  of  this  method  do  not  compare  well  with 
the  E  =  1500  CBR  equation.  The  derived  moduli  are  con¬ 
sistently  larger. 

The  second  method  Weiss  proposed  is  the  dynamic  load- 
deflection  curve  method.  Briefly,  this  method  entails  using 
a  non-linear  dynamic  response  theory  that  predicts  surface 
deflections  given  a  subgrade  modulus  and  pavement  moduli. 
When  the  measured  surface  deflections  agree  with  the 
predicted  ones,  the  correct  modulus  has  been  chosen. 
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The  results  of  this  method  generally  agree  with  the  E  =  1500 
CBR  formula. 

Weiss  noted  that  the  subgrade  modulus  is  a  function  of 
the  confining  pressure  of  the  soil  (and  hence,  the  overbur¬ 
den  pressure).  Consequently,  he  went  to  some  length  to 
remove  the  effect  of  the  weight  of  the  vibrator.  Wiseman 
(1973)  noted  this  same  effect  when  comparing  the  coeffi¬ 
cients  of  subgrade  reaction  gained  from  tests  using  the 
Benkelman  beam  and  the  Road  Rater.  The  Road  Rater  produced 
consistently  higher  values  of  the  coefficient,  which  Wiseman 
attributed  to  the  added  confining  pressure  on  the  subgrade. 
He  had  some  success  using  the  Hertz  theory  of  predicting  the 
coefficients  of  the  surface  course  and  the  underlying  soils. 

Witczak  (1980)  ran  vibratory  tests  on  frozen  and  thaw¬ 
ing  soils.  He  noted  that  the  subgrade  modulus  is  heavily 
dependent  on  the  temperature  of  the  subgrade.  However,  h e 
found  that  the  subgrade  modulus  did  not  have  much  effect  on 
the  maximum  deflection  measured  by  a  vibrator. 

Vaswani  (1971)  used  surface  deflections  measured  by  a 
Dynaflect  to  predict  the  modulus  of  subgrade  reaction.  He 
measured  the  maximum  deflection  and  the  deflection  profile 
in  front  of  the  maximum  deflection,  and  then  used  these  to 
compute  the  “spreadab i 1 i ty "  -  an  average  deflection  of 
sorts.  From  his  study,  he  drew  up  graphs  that  could  be  used 
to  approximate  the  modulus 
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All  of  the  methods  reviewed  require  contact  with  the 
pavement  in  order  to  ascertain  the  desired  property.  No 
non-contact  methods  were  reported. 


C .  Variables  Affecting  Pavement  P  er  f or mane  e 

This  section  is  intended  to  examine  the  variables  that 
affect  the  measurements  described  in  the  last  section.  The 
list  of  possible  factors  is  a  long  one  and  is  as  varied  as 
the  human  imagination  that  produces  them.  Consequently,  to 
be  expedient,  this  review  will  at  temp  t  to  cover  only  the 
more  obvious  factors. 

Probably  the  most  obvious  factor  affecting  pavement 
performance  is  the  asphalt  concrete  itself.  Thicker  pave¬ 
ments  deflect  less  under  a  given  load  than  do  thinner  pave¬ 
ments  -C<H'ihter  and  Harr  <1975>>.  One  would  expect  that 
less  deflection  would  mean  a  longer  life  for  the  pavement. 

If  all  other  variables  were  held  constant,  this  would  prob¬ 
ably  be  the  case  {(Nichols  (  19635  3- 

Temperature  is  an  important  factor  in  determining  when 
asphalt  will  show  distress  under  loading.  Trailer  and  Lay¬ 
man  (1975)  report  that  the  glass  transition  temperature  for 
some  asphalts  is  not  necessarily  below  32°  F.  This  is  the 
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temperature  at  which  asphalt  becomes  a  brittle  material. 
Kandahl  and  Wenger  (1975)  found  that  the  temperature  at 
which  asphalt  pavement  was  placed  affected  its  performance. 
They  also  noted  that  high  viscosity  asphalt  pavements 
(mixed  at  low  temperatures)  retained  entrained  air  longer 
than  low  viscosity  asphalt  pavements/  since  they  were 
compacted  more  slowly. 

Kondner  and  K.rizek  (1966)  noted  that  increased  tempera¬ 
ture  leads  to  increased  flexibility  of  the  pavement. 
Southgate  and  Deen  (1969)  also  noted  that  the  deflection  of 
an  asphalt  pavement  increases  with  temperature.  They  dev¬ 
ised  a  method  to  adjust  the  deflection  to  a  standard  tem¬ 
perature  by  noting  that  the  relationship  between  deflection 
and  the  logarithm  of  temperature  was  linear.  They  then 
scaled  the  measured  deflection  by  a  derived  adjustment  fac¬ 
tor  to  get  the  corrected  deflection. 

Lister  (1972)  tested  asphalt  pavements  at  different 
temperatures  and  found  a  linear  variation  of  deflection  for 
temperatures  between  10°  C  and  30°  C,  although  the  rate  of 
change  varied  with  pavement  type  (Figure  1 1  —  1 ) . 

Peterson  and  Shepherd  (1972)#  in  their  study  with  a 
Dynaflectz  were  unable  to  find  any  correlation  between  tem¬ 
perature  and  deflection.  They  do>  however/  provide  a  good 
review  of  the  temperature  corrections  proposed  by  other 
organizations  They  are  listed  below  - 


34 


X)  Canadian  Goad  Roads  Association  —  they  general¬ 
ized.  saying  that  the  change  in  deflection  was 
0.0002  inches  for  every  10°  F, 

2)  District  of  Northern  Vancouver  -  uses  a  linear 
relation  between  deflection  and  temperature. 

3>  AASH0  -  above  30'J  F  -  no  effect  of  temperature  on 
deflection  j  below  30°  F  -  curvilinear  effect  of 
temperature  on  deflection, 

4>  Transport  and  Road  Research  Laboratory  —  above  86° 
F  -  a  decrease  is  noted  above  this  temperature,  - 
below  86°  F  —  deflection  increases  approximately 
7.2*:  per  °F, 

5)  Colorado  -  linear  relation  between  temperature  and 
deflection,  but  different  at  every  site. 

The  U.  S.  Army  Corps  of  Engineers  (1975)  conducted  a 
study  using  their  16-kip  vibrator  on  asphalt  pavements.  By 
testing  at  different  temperatures  on  different  thicknesses 
of  pavements,  they  were  able  to  establish  the  correlation 
shown  in  their  Figure  1  <p75).  While  they  felt  this  was 
good  for  their  particular  site,  they  expressed  reservations 
about  using  it  at  other  sites. 


Dorman  and  Metcalf  (1965)  studied  the  effect  of  tem¬ 
perature  on  the  modulus  of  the  pavement  system.  They  found 
that  the  modulus  decreased  with  increasing  temperature. 

Meyer  and  Kummer  (1969)  noted  the  effect  of  temperature 
on  the  skid  resistance  of  pavements.  They  broke  skid  resis¬ 
tance  into  two  components  -  adhesion  and  hysteresis.  The 
adhesion  component  may  increase  or  decrease  with  tempera¬ 
ture.  while  the  hysteresis  component  only  decreases  with 
temperature  increases. 

Several  investigators  have  found  that  the  season  of  the 
year  influences  pavement  properties.  Using  deflection  as  a 
criterion,  Kondner  and  Krizek  (1966)  found  a  sinusoidal 
variation  with  time.  Poehl  and  Scrivner  (1971)  also  found 
that  deflections  varied  sinusoidally  during  the  year,  with 
the  largest  deflections  occurring  in  the  springtime.  They 
also  noted,  as  did  Peterson  and  Shepherd  (1972),  that  rain¬ 
fall  increases  deflections.  Yoder  (1962)  reported  that, 
among  the  many  factors  affecting  deflection,  climate  and 
temperature  are  important 

The  magnitude  of  the  load  applied  to  a  pavement  is  an 
important  factor  to  consider.  This  is  why  highway  depart¬ 
ments  specify  load  limits,  usually  18  kip?  per  axle  (static 
load).  Whittemore  (1969)  conducted  an  interesting  study  on 
the  dynamic  load  of  tires  on  pavements  Using  accelerome¬ 
ters  mounted  on  the  axle  of  a  truck,  he  measured  the  force 
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applied  by  the  tires  to  the  pavement  as  the  truck  moved  down 
the  road.  He  found  there  was  a  large  variation  in  the  force 
-  up  to  5000  pounds.  Compared  to  the  design  load  of  18.000 
pounds:  this  is  a  large  variation! 

The  last  factor  to  be  considered  is  site  dependency. 
This  variable  takes  into  account  all  the  variables  that  can¬ 
not  be  measured  or  accounted  for  directly;  particularly* 
weather  history*  fluctuations  of  the  groundwater  table*  con¬ 
struction  control  (eg.  variation  in  the  materials  used* 
mode  of  placement,  method  and  degree  of  compaction,  initial 
water  content  of  the  subsoils,  locked-in  stresses  in  the 
pavement  and  subsoils),  animal  activity  and  history  of  load¬ 
ing.  These  variables,  most  of  which  are  difficult  to  or  are 
impossible  to  quantify,  are  very  likely  to  be  influential  in 
determining  whether  a  pavement  fails  or  not. 

Kondner  and  Krizek  (1966)  acknowledged  the  importance 
of  site  dependency.  After  an  exhaustive  study  of  the  AASHO 
Road  Test  data,  in  which  correlations  between  PSI  and  other 
factors  were  derived,  they  concluded  that  derived  relation¬ 
ships  are  good  only  at  the  site  where  the  data  were  gath¬ 
ered.  Extrapolation  to  other  sites  would  not  be  likely  to  be 
val  id. 

Yoder  (1962)  showed  that  the  amount  of  deflection  under 
a  given  load  varies  with  the  lateral  placement  of  the  load 
on  the  pavement  Peterson  and  Shepherd  (1972)  noted  that 
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Colorado's  deflection  cr  iter  id  are  different  at  every  site. 

Poehl  and  Scrivner  (1971)  presented  a  very  good  example 
of  site  dependency  in  their  study  of  pavements  in  Texas. 

They  used  a  Dynaflect  to  monitor  1000  foot  sections  of  road 
in  five  areas  of  Texas.  Their  conclusions  u>ere  - 


1)  the  odds  are  2: 1  that  the  annual  mean  deflection 
at  the  ends  of  the  1000  foot  strip  are  very  dif¬ 
ferent, 

2)  the  odds  are  2: 1  that  the  annual  mean  change  in 
deflection  at  the  ends  of  the  1000  foot  strip  are 


also  very  different, 

3)  in  a  mile  of  "uniform"  pavement,  the  odds  are  high 
that  the  deflection  measured  at  different  points 
on  the  same  day  will  vary  more  than  the  deflection 
measured  at  any  one  point  during  a  year,  and 

4)  within  a  1000  foot  strip,  the  annual  percentage 
change  in  deflection  varied  s i gn i f i can 1 1  y . 


Finally,  they  stated  that  while  the  Dynaflect  does  not  meas¬ 
ure  the  same  deflection  as  a  static  or  moving  wheel  load,  it 
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does  provide  a  means  of  comparing  pavements. 

Vedros  and  Barker  (1977)  presented  the  results  of  their 
study  done  in  Kentucky.  They  compared  two  sections  of  pave¬ 
ment  within  51  feet  of  each  other.  One  had  cracked  and  the 
other  had  not.  Extensive  field  and  laboratory  tests  were 
made  to  determine  the  cause.  Moisture.  CBR.  density,  cyclic 
plate  tests,  deflections,  vibratory  measurements,  and  sur¬ 
face  profiles  were  all  checked.  Gradation,  resilient 
modulus,  creep  behavior  and  rutting  susceptibility  were  also 
determined.  They  concluded  - 


"Data  analysis  indicates  that  there  were  no  apparent 
differences  between  test  pits  to  explain  why  the  pave¬ 
ment  would  be  cracked  at  one  test  pit  but  not  at  the 
other.  " 


They  did  find,  however,  that  the  resilient  deformation  of 
the  pavement  could  be  predicted  by  a  finite  element  computer 


program. 
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D.  Interpretation  o f  Da ta  Qb ta i ned  from  Non— Destructive 
Tests 


1.  Texture 

Quantitative  measures  of  texture  are  usually  given  in 
terms  of  skid  numbers  (SN),  rather  than  absolute  measures  of 
the  texture.  Kummer  and  Meyer  (1967)  provide  a  comprehen¬ 
sive  review  of  skid  number  determination,  detailing  the 
several  devices  in  use  and  giving  correlations  between  them. 
Their  recommended  skid  number,  37,  as  determined  by  the  ASTM 
E-274  test,  falls  in  the  range  of  values  in  use  in  the 
United  States.  In  the  same  paper,  a  correlation  (with  skid 
numbers)  is  given  for  the  aluminum  foil  texture  test  men¬ 
tioned  in  the  previous  section.  They  noted  that  about  10 
holes  per  square  inch  are  required  to  gain  sufficient  skid 
resistance  for  medium  speed  roads. 

McCullough  and  Hankins  (1966)  studied  pavements  in 
Texas  with  a  skid  trailer.  They  recommended  that  two  skid 
numbers  (coefficients  of  friction)  be  specified,  each  meas— 
ui  id  at  a  particular  speed.  For  20  mph  and  50  mph,  they 
recommended  0.31  and  0.24.  respectively,  as  minimum  values. 
Pavements  that  fall  below  these  should  be  repaired. 

Beaton  (1971)  examined  the  texture  of  concrete  pave¬ 
ments  He  noted  that,  of  the  several  methods  available  for 

"te a  curing”  fresh  c  ncrete.  the  California  Highway 
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Department  prefers  the  broom  method.  However,  no  quantita¬ 
tive  evaluation  was  given.  Grooving  of  concrete 
pavements,  implemented  to  increase  drainage  and  road  tire 
friction,  is  done  in  California  by  1/8  in.  slots  on  3/4  in. 
centers  parallel  to  traffic  flow 

Rose  and  Galloway  (1977)  studied  the  effect  of  depth  of 


water  on  skid  resistance.  After  testing  various  combina¬ 
tions  of  speed,  tires,  tire  pressures,  pavements  and  water 
depths,  they  recommended  that  the  macrotexture  of  the  pave¬ 
ment  be  greater  than  0  04  in.  for  high  speed  roads.  Moore  (1969) 
suggested  that  the  microtexture  of  the  pavement  should  be  in 
the  range  0.0004  in  to  0.004  in.  to  insure  adequate  skid 
resistance  He  noted  (in  1969)  that  there  were  not  any 
methods  to  evaluate  this. 

2.  Roughness 

The  measurement  of  pavement  roughness,  as  was  discussed 
in  section  II,  S,  2,  c,  2.  takes  many  forms.  The  data 
obtained  from  the  different  devices  have  their  own  criteria 
of  acc ep tab  i  l  i ty .  Because  changes  in  PSI  reflect  changes  in 
roughness,  it  is  only  natural  that  developed  devices  should 
try  to  correlate  with  it.  Baker  (1975)  affirmed  that  this 
is  the  case.  He  noted  that  regression  equations  relating 
the  measured  variables  are  used  by  highway  departments  to 


derive  the  PSI. 


Because  the  requirements  of  each  department 
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and  the  peculiarities  of  their  measuring  techniques  vary, 
each  department  derives  its  own  equation.  The  Roughometer 
and  the  roadmeter  are  commonly  used  with  regression  equa¬ 
tions.  Roughometer  values  of  90  in. /mile  for  rigid  pave¬ 
ment  and  75  in.  /mile  for  flexible  pavements  are  sometimes 
used  as  construction  criteria.  The  CHLOE  pr of i 1 ome ter  ' s 
output,  slope  variance,  is  also  used  in  a  regression 
analysis  to  derive  the  PSI. 

The  GRM  prof i lometer,  described  in  Highway  Research 
Board  Special  Report  No.  95  <1968)  gives  the  pavement  pro¬ 
file,  but  not  a  measure  of  user  ac c ep tab i 1 i ty .  That  must  be 
determined  subjectively  fay  the  persons  conducting  the  study. 
The  non— contact  methods  of  pavement  profile  measurement  all 
require  a  subjective  correlation  with  user  acceptab  i  1  ity. 
Neither  Joyce  (1975)  nor  Dickerson  and  Mace  (1976)  suggested 
any  criteria  for  the  acceptability  of  pavement  roughness. 

Houbolt  (1961)  suggested  that  a  vertical  change  in  pro¬ 
file  of  more  than  0.  03  ft  per  250  ft  is  too  rough  for  air¬ 
planes.  When  using  a  straight-edge  to. measure  roughness,  he 
used 
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s  <  ■: 
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wh  ere 


6  —  deviation  from  straight-edge  <ft) 
L  =  length  of  straight-edge  (ft) 


Lastly*  he  proposed  that  acceleration  experienced  by  an  air¬ 
plane  can  be  used  as  a  criterion.  He  suggested  that  an 
upper  limit  of  0.  3  g  be  used.  Hall  and  Kapelson  (1962)  used 
0.  5  g.  They  also  noted  that*  although  a  particular  bump  may 
not  cause  the  maximum  ac c e 1 er a t i on.  it  may  do  so  when  run 
across  in  conjunction  with  previous  bumps,  since  the  air¬ 
plane  continues  to  vibrate  after  the  bump  has  been  passed. 
Eremin  (1962),  in  a  discussion  of  Houbolt,  noted  that  0  3  g 
is  less  than  the  design  acceleration  used  in  airport  pave¬ 
ment  construction. 


Sonnenberg  (1978)  used  a  statistical  analysis  of  air¬ 
port  profiles,  filtering  out  the  long  wavelengths#  to  get 
the  standard  deviation  of  the  profile.  Although  he  found  a 
range  of  standard  deviations  of  3  to  13  mm  for  the  13  air¬ 
ports  examined*  he  failed  to  declare  what  would  be  con¬ 
sidered  acceptable 


Walker  and  Hudson  (1973)  related  roadmeter  readings  to 
a  servi ceab i  1  i ty  index  thus 


SI 


5  « 


lr. 


5 
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where 

SI  —  serviceability  index 

M  =  roadmeter  roughness  reading  < in. /mile) 

B  =  roadmeter  instrument  coefficient  (calibration 
factor  for  test  vehicle) 

Dann  and  Schultz  (1978)  in  their  study  of  Wisconsin 
roads  with  a  roadmeter,  developed  the  relationship  between 
roadmeter  readings  and  PSI  shown  in  Figure  II-2  .  Chong  and 
Phang  (1978)  carried  their  roadmeter  study  further,  relating 
the  roadmeter,  pro f i 1 ometer ,  and  roughometer  to  either  the 
PSI  or  the  PPR  (Present  Performance  Rating  used  by  the  Cana¬ 
dian  Good  Roads  Association). 

3.  Deflection 

Finding  an  absolute  measurement  that  will  predict  when 
a  pavement  will  fail  is  difficult.  Deflection  studies, 
nevertheless,  are  often  employed  to  this  end.  This  section 
reviews  some  of  the  criteria  put  forth  in  the  literature 
that  attempt  to  define  when  a  pavement  will  fail. 

Highter  and  Harr  (1975)  studied  distress  in  airport 
pavements.  By  noting  the  number  of  airplanes  that  used  a 
runway,  the  maximum  deflection  caused  by  each  and  the  condi¬ 
tion  of  the  pavement,  they  wc-re  able  to  relate  the  sum  of 


Roadmeter  statistic 
(counts  per  mile) 


Figure  II-2  Pavement  evaluation  using  the  Roadmeter 

statistic  (after  Dunn  and  Schultz  (1972)) 
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the  deflections  the  pavement  experienced  to  the  PSI  Their 
hypothesis  -  that  the  condition  of  the  pavement  could  be 
related  to  the  cumulative  peak  deflections  -  was  shown  to  be 
valid.  Typical  results  are  shown  in  Figure  II-3.  The  three 
curves  represent  various  combinations  of  thicknesses  of  sur¬ 
face  courses,  bases  and  subbases. 

Huang  (1971)  proposed  that  the  magnitude  of  the  strain 
be  used  as  the  criterion  to  determine  when  asphalt  pavements 
will  crack.  By  measuring  the  maximum  curvature  and  the 
deflection  caused  by  a  wheel  load/  he  developed  a  method  to 
calculate  the  strain.  From  laboratory  investigations,  he 
determined  the  allowable  strain.  By  comparing  the  two.  it 
can  be  seen  if  the  pavement  will  crack.  He  noted  that  the 
time  of  year  that  the  measurements  are  made  also  affects  the 
results.  To  predict  rutting,  the  measurements  should  be 
made  in  the  summer.  To  predict  cracking,  winter  measurements 
should  be  used. 

Lister  (1972)  found  that 


life. 


where  d  —  def  lect”  ion  of  pavement 
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Figure  II-3  Cumulative  Deflection  vs.  PSI  (after 
Highter  and  Harr  (1975)) 


early  life  deflection 


4  7 


critical  life  in  cumulative  standard  axles  (I06) 

Figure  II-4  Summary  of  relation  between  life  and 
deflection  (after  Lister  (1972)). 
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Peterson  and  Shepherd  '1972)  used  a  Dynaflect  to  deter¬ 
mine  which  layer  beneath  a  pavement  needed  repair.  They 
defined  three  parameters  —  the  Dynaflect  Maximum  Deflection 
( DMD ) »  the  Surface  Curvature  Index  (SCI)  and  the  Base  Curva¬ 
ture  Index  (BCI).  The  DMD  is  the  maximum  deflection  meas¬ 
ured  by  the  transducer  nearest  the  vibrating  wheels  (the 
first  sensor).  The  SCI  is  the  difference  between  the 
deflections  measured  by  first  and  second  sensors,  and  the 
BCI  is  the  difference  between  the  fourth  and  fifth  sensor 
deflection.  The  SCI.  a  measure  of  curvature,  is  used  to 
define  the  strength  of  the  surface  layers.  The  BCI  is  used 
to  define  the  strength  of  the  lower  layers.  Some  examples 
of  how  the  BCI  and  SCI  are  interpreted  are  shown  in  Figure 
II-5.  Figure  II-6  shows  how  they  evaluate  pavements  using 
their  parameters.  They  then  developed  nomographs  to  compute 
overlay  thicknesses. 

A  comprehensive  review  of  deflection  analysis  methods 
is  presented  by  Beca.  et  al.  (1974).  They  give  the  maximum 
deflection  criterion  used  by  different  organizations.  These 
a^e  listed  in  Table  1 1  —  1  All  measurements  were  made  in  the 


springtime. 


DMD 

SCI 

BCI 

CONDITION  OF  PAVEMENT  STRUCTURE 

PAVEMENT  AND  SUBGRADE  WEAK 

SUBGRADE  STRONG,  PAVEMENT  WEAK 

SUBGRADE  WEAK,  PAVEMENT  MARGINAL 

LE.a*a 

LCO.U 

DMD  HIGH,  STRUCTURE  OK 

STRUCTURE  MARGINAL ,  DMD  OK 

% 

PAVEMENT  WEAK ,  DMD  OK 

N 

SUBGRADE  WEAK,  DMD  OK 

V-'/ 

PAVEMENT  AND  SUBGRADE  STRONG 

- - - 1 

KEY 

G.T.  -  greacer  than 
G.E.  »  greater  than 
or  equal  to 
L.T.  ■  less  than 
L.E.  -  less  than 

or  equal  to 


Figure  II-5  Evaluation  using  DMD,  SCI 

and  BCI  (after  Peterson  and 
Shepherd  (  1972) ) 


2 


Subgrade  Strong  ,  Pavement  Weak 


Subgrade  Weak,  Pavement  Marginal 


Figure  II-6  Evaluation  of  pavement  (after 
Peterson  and  Shepherd  (1972)) 


51 


TABLE  I I— 1 


Deflection  Criteria 


INCHES  OF  DEFLECTION 


AUTHOR 


0.05  for  0.  5  in.  asphalt 
to  0.  01S  for  4  in.  asphalt 

0.  05 

0.  036 

O.  020 


Hveem 


Canadian  Good  Roads 
Assoc iat  ion 

Virginia  Highway  Research 
Counc i 1 

Transport  and  Road 
Research  Laboratory 


They  noted  that  curvature  is  also  used  as  a  criterion, 
although  they  could  not  decide  if  this  was  a  better  measure 
than  deflection.  They  contended  that  a  given  deflection 
cannot  be  related  to  a  particular  degree  of  distress  or  loss 
of  serviceability. 

Green  and  Hall  (1975)  claimed  that  deflection  is  not  a 
good  indicator  of  the  remaining  life  of  a  pavement.  They 
claimed  that  there  is  no  discernible  change  in  deflection  up 
until  failure.  They  did  say,  however,  that  deflection  was  a 
good  measure  of  performance.  Moore,  Hanson  and  Hall  (1978) 
cited  a  study  by  Vaswani  in  Virginia  that  showed  that  the 
deflection  varies  over  the  life  of  the  pavement  as  shown  in 
F.gur*  1 1— 7  Region  I  of  Figure  II-7  represents  the  densif- 

lcation  that  takes  place  after  the  pavement  is  placed.  The 


deflection  then  remains  constant  (region  II)  until  it 
increases  to  failure  (region  III)  This  shape  of  figure  is 
called  a  “bathtub  function". 

Anderson  (1976)  used  the  results  of  studies  by 
Southgate  and  Vaswani  et  al.  to  evaluate  pavements  in  Utah. 
Southgate's  study  developed  the  method  for  temperature 
correction  discussed  in  Section  II.  C.  After  correcting  the 
deflections  for  temperature.  Anderson  used  Vaswani's  sprea— 
dability  concept  to  get  the  modulus  of  the  subgrade  and  the 
average  modulus  of  the  system.  They  are  determined  from  the 
chart  given  in  the  paper. 

Croney  (1772)  established  a  rough  correlation  between 
PSI  and  rut  depth  (permanent  deflection).  Using  a  six  foot 
straight-edge  laid  across  the  pavement,  he  found  that  a  rut 
depth  of  10  mm  was  the  most  common  one  in  which  cracks  were 
found . 

4.  Modulus  of  Elasticity  of  the  Subgrade 

The  modulus  of  elasticity  of  the  subgrade  is  used  for 
the  design  and  evaluation  of  pavements.  Methods  of  obtain¬ 
ing  E  were  reviewed  in  Section  II.  B.  2.  c.  4.  This  subsec¬ 


tion  examines  bow  different  investigators  use  the  modulus. 
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Dorman  and  Metcalf  (1965)  used  layered  system  theory  to 
develop  design  curves  for  asphalt  pavements.  They  were  con¬ 
cerned  with  the  tensile  strains  on  the  bottom  of  an  asphalt 
concrete  layer  subjected  to  a  wheel  load,  and  the  effect  of 
temperature,  thickness  of  the  asphalt  concrete,  and  modulus 
of  elasticity  of  the  suijrade  on  those  strains.  The  design 
curves  they  proposed  use  the  expected  strains  and  the  thick¬ 
ness  of  the  base  to  compute  the  thickness  of  asphalt  con¬ 
crete  required. 

Huang  (1971)  was  also  interested  in  the  tensile  strains 
at  the  bottom  of  asphalt  concrete  layers.  He  modeled  the 
pavement  in  three  layers,  each  having  its  own  modulus.  The 
lowest  layer  was  the  subgrade.  If  the  ratios  of  the  moduli, 
the  layer  thicknesses,  and  the  '‘curvature"  are  known,  the 
tensile  strains  can  be  computed.  These  are  then  checked 
against  allowable  strains  to  see  if  the  pavement  is  accept¬ 
able,  or  is  in  need  of  repair. 

Methods  of  defining  the  modulus  vary  from  study  to 
study.  Crawford  and  Katana  (1975)  discussed  different  ways 
to  get  E  from  stress-strain  plots.  Then,  having  chosen  a 
modulus,  they  show  how  it  is  used  in  finite  element  studies 
of  pavement  response.  A  variety  of  different  element  types 
and  configurations  are  considered 

Layered  elastic  theory  is  sometimes  used  to  analyze 
pavements,  despite  the  fact  that  pavements  and  soils  do  not 
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behave  elastically  in  all  stress  ranges.  Use  of  this  theory 
requires  a  knowledge  of  the  thickness  of  the  layers  and 
their  material  properties.  Weiss  (1980)  used  a  pavement 
vibrator  to  determine  the  subgrade  modulus  and  then  used  it 
with  layered  elastic  theory  to  determine  the  bearing  capa¬ 
city  of  the  pavement.  Vaswani  (1971)  used  a  Dynaflect  to 
get  a  measure  of  the  deflection  basin  he  calls  "spr eadab  i  1- 
ity".  With  this  and  the  various  layer  modulii  he  developed 
charts  to  evaluate  the  pavement.  Yoder  and  Witczak  (1975) 
note  that  the  modulus  of  elasticity  is  used  in  both  the 
Asphalt  Institute  and  FAA  methods  of  design  of  airport  pave¬ 
ments. 

The  modulus  of  subgrade  reaction  is  sometimes  used  in 
the  evaluation  of  concrete  pavement  Vesic  and  Saxena 
(1970)  used  the  subgrade  modulus  of  elasticity  in  conjunc¬ 
tion  with  the  concrete  slab  properties  to  derive  a  formula 
for  the  modulus  of  subgrade  reaction.  This.  in  turn,  was 
used  to  predict  the  deflections  of  the  slab  under  load. 


This  review  has  shown  that  there  are  a  great  number  of 
ways  to  test  pavements  available  to  the  pavement  engineer. 
The  divergence  of  opinions  on  the  subject  reflects  not  only 
the  uncertainty  in  present  methods,  but  also  the  need  to 
account  for  regional  factors  This  study  hopes  to  provide  a 


method  whereby  many  of  these  regional  factors  can  be 
accounted  for.  through  freguent  testing. 
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III.  DEVELOPMENT  OF  THE  SYSTEM 


A.  Introduction 


The  rapid  non-contact  measurement  of  pavement  deflec¬ 
tions  under  actual  wheel  loadings  requires  a  complex  state- 
of-the-art  system  of  electronic  components.  The  system  has 
to  be  capable  of  making-  recording,  and  processing  large 
amounts  of  data,  and  doing  so  very  quickly.  Such  a  device 
would  not  have  been  possible  ten  years  ago,  but,  because  of 
recent  advances  in  electronics,  it  is  possible  today. 

The  requirements  are  that  the  system  make  non-contact 
measurements  of  the  deflected  surface  of  asphalt  pavements 
caused  by  moving  wheel  loads  Furthermore,  the  device  is  to 
be  transported  by  the  vehicle  causing  the  deflection.  Impli¬ 
cit  in  these  requirements  are  the  following  details  - 

1.  The  device  must  be  sufficiently  sensitive  to  detect 
deflections  on  the  order  of  0.001  inch.  This  require¬ 
ment  was  determined  from  the  magnitude  of  the  deflec¬ 
tions  recorded  in  previous  studies  -C(Baladi  (1976); 


Ng-A-Qui  <1?76)>. 


56 


2  The  device  must  be  insensitive  to  changes  in  the 
color  of  the  surface. 

3.  The  device  must  be  insensitive  to  ambient  light  and 
sound . 

4.  The  system  must  be  able  to  provide  not  only  a  measure 
of  the  pavement  condition,  but  also  the  location  and 
speed  of  the  vehicle  at  the  time  of  measurement. 

5.  The  system  must  be  able  to  record  and  process  the 
data  gathered. 

6.  The  system  should  be  relatively  easy  to  set  up  and 
use. 

7.  The  range  of  the  device  must  be  large  enough  to 
accommodate  the  expected  vehicle  movement  plus  the 
pavement  deformation. 

Most  of  these  criteria  were  satisfied  by  the  selection 
of  the  non-contact  distance  measuring  gages,  to  be 
described  in  section  III,  C,  1  The  main  components  of  the 
system  are  the  electro- optical  distance  measuring  gages,  a 
rigid  beim,  a  fifth  wheel,  and  the  attendant  electronics. 
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el^ct’o-oprical  distance  iiieasur  ing  gages  are  mounted 
on  a  rigid  beam  which  in  tui-n  is  mounted  longitudinally  on  a 
load  vehicle  (see  Figure  1-1)  The  gage  heads  measure  the 
distance  from  the  beam  to  the  pavement  in  a  non-contact  way 
(using  lasers)  They  are  the  very  heart  of  this  system. 

In  order  jr  the  data  analysis  system  to  work,  the 
gages  cannot  move  relative  to  their  original  alignment.  For 
this  reason,  they  are  mounted  on  a  rigid  beam. 

The  data  analysis  system  requires  each  gage  to  read  the 
distance  from  the  rigid  beam  to  the  pavement  surface  at  the 
same*  horizontal  locations.  Each  succeeding  gage  must  read 
where  the  preceding  gages  have  read  A  distance  measuring 
device  called  a  "fifth  wheel"  is  used  to  provide  this  meas¬ 
ure.  The  speed  of  the  vehicle  is  determined,  indirectly, 
from  the  fifth  wheel 

The  attendant  electronics  are  the  power  supply,  com¬ 
puter,  and  the  electronics  for  the  fifth  wheel. 


*  In  the  statistical  sense 
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2  Method oloou 

1.  Theory  for  Profile,  Deflection  and  Texture  Measurements 

It  was  desired  to  measure  the  deflected  profile  of  the 
pavement  ad jacent  to  the  wheel  load  as  well  as  the  unde¬ 
flected  profile.  The  solutions  to  these  problems  go  hand  in 
hand.  Dickerson  and  Mace  <1976)  described  an  algorithm  that 
was  adapted  for  the  present  study.  Its  description  follows. 
The  computer  program  that  uses  the  solution  is  described  and 
listed  in  Appendix  C. 

The  algorithm  that  determines  the  undeflected  profile  of 
the  pavement  uses  five  values  -  three  gage  measurements  and 
two  arbitrary  elevations  above  a  datum.  These  five  inputs 
are  then  used  to  calculate  a  new  elevation  above  the  datum 
while  the  vehicle  moves  to  a  new  location.  In  Figure  III-l. 
the  two  arbitrary  elevations  are  designated  as  BB  and  CC, 
and  the  three  gage  head  measurements  are  A,  B>  and  C.  Given 
these  five  inputs,  the  algorithm  calculates  the  elevation, 

AA,  in  front  of  the  known  elevations.  All  five  inputs  are 
obtained  from  the  undeflected  portion  of  the  pavement.  AA 
defines  a  point  on  the  undeflected  profile.  For  this  reason, 
the  front  three  gages  are  called  the  "undeflected  profilome- 
ter".  When  elevations  are  determined  in  the  area  influenced 
by  the  wheel  load,  "deflected  profile"  points  are  calcu¬ 
lated  A s  any  five  inputs  may  be  used,  (if  chosen  in  a 
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direction  of  motion 


Figure  III-  I  Schematic  of  the  system.  Load  wheel 
not  shown  (normally  at  D) . 
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manner  consistent  with  that  just  described),  the  elevation 
can  be  calculated  at  gage  D>  near  the  load  wheel. 

The  derivation  of  the  algorithm  follows  from  Figure 
III— 1.  Here,  as  above.  BB  and  CC  are  arbitrarily  assigned, 
and  A.  B.  and  C  are  actual  gage  readings. 

By  constructing  00'  parallel  to  the  datum,  it  follows 

that 

C+CC=B+BB+b=A+  AA  +  a.  (1) 

If  the  gages  are  equally  spaced  on  the  beam 


a  =  2  b  (2) 

Assuming  this,  and  noting  that 

b  =  C  +  CC  -  B  -  BB  (3) 

gives  (by  substitution) 

AA  «  C  +  CC  -  A  -  2  (  C  +  CC  -  B  -  BB  ).  (4) 

which  is  the  desired  quantity,  the  new  undeflected  eleva¬ 
tion.  AA  can  be  calculated  from  Eq.  (4)  as  all  the  quanti¬ 
ties  an  the  right  side  of  the  equation  are  known  (either 
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measured  or  assumed).  This  provides  the  initial  step  of  the 
procedure.  The  process  is  then  continued  in  the  following 
manner:  the  gages  move  forward  until  gage  B  is  over  the  pre¬ 
vious  position  of  gage  A*,  and  the  algorithm  is  repeated# 
using  the  five  designated  quantities  (BB  is  the  previous 
AA>  and  CC  is  the  previous  BB ) .  A  new  undeflected  profile 
point  is  then  calculated.  Repetition  eventually  yields  the 
entire  undeflected  profile  at  increments  corresponding  to 
the  gage  spacing. 

The  spacing  of  the  gages  on  the  beam  is  very  important. 
By  mounting  all  the  gage  heads  at  equal  intervals  along  the 
beam#  the  gage  head  readings  can  be  timed  so  that  each 
reads  at  the  same  location  as  the  previous  one.  Thus,  any 
three  gages  can  be  used  as  a  p ro f i 1 ome t er .  By  placing  gage  D 
an  integral  number  of  gage  spacings  from  gages  A,  B,  and  C, 
gages  B,  C,  and  D  can  be  used  as  a  prof i lometer.  Whereas  AA 
was  calculated  in  a  "forward"  manner,  DD  is  calculated 
in  a  "backward"  manner.  The  five  inputs  for  the  latter  are 
B,  C,  D,  BB,  and  CC.  DD  is  calculated  from 

DD  =  C  +  CC  -  D  +  (  r  )  (  C+CC-B-BB  )  (5) 

where  r  is  the  ratio  of  the  distance  between  gages  C  and  D 
to  the  distance  between  gages  B  and  C.  When  the  B-C-D  profi- 
lometer  is  used,  a  deflected  profile  is  calculated, 

*  In  the  statistical  sense 
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gage  D  being  near  the  load  wheel. 

The  deflections  (relative  to  the  undeflected  profile) 
caused  by  the  load  wheel  can  now  be  calculated.  The  unde¬ 
flected  profile  was  calculated  by  the  A-B-C  prof i lometer 
The  deflected  profile  was  calculated  by  the  B-C-D  profilome- 
ter.  The  fifth  wheel  assures  that  all  gages  read  at  the 
same*  points  on  the  pavement.  So<  two  profiles  are  obtained. 
The  difference  between  the  two  provides  the  deflections 
caused  by  the  load  wheel. 

Two  important  pavement  parameters  con  be  ascertained 
from  the  measurements  described  above  -  the  profile,  and  the 
deflection  adjacent  to  the  wheel.  Because  the  gages  can  make 
measurements  at  a  very  rapid  rate  (16  kHz),  another  impoi — 
tant  pavement  characteristic  -  the  pavement  texture  -  can 
also  be  determined.  The  texture  can  be  obtained  from  the 
variations  in  readings  in  the  near  vicinity  of  a  point  on 
the  pavement.  This  variation  is  possible  because  the  laser 
gages  that  were  used  had  a  light  spot  size  much  smaller  than 
the  variation  in  the  surface  texture,  and  because  the  gage 
resolution  is  smaller  than  the  surface  texture. 

Many  readings  are  required  over  a  short  distance  not 
only  so  that  the  texture  can  be  measured,  but  also  to  keep 
null  readings  due  to  "dropout"  from  entering  the  profile 
algorithm.  Dropouts  occur  when  the  light  spot  from  the  gage 

*  In  the  statistical  sense. 
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is  hidden  from  the  view  of  the  photodetector  (Figure  III-2). 
In  o^der  to  keep  from  losing  the  data  point.  the  gage  head 
makes  many  readings  over  a  "short  distance")  averages  them, 
and  considers  the  average  to  be  the  distance  from  the  gage 
to  the  pavement.  Four  inches  of  the  pavement  in  the  line  of 
travel  was  chosen  as  the  short  distance  in  this  study. 

Averaging  many  readings  provides  a  very  stable  statisti¬ 
cal  measure  of  the  distance  from  the  gage  to  the  pavement. 
The  large  number  of  readings  allows  an  extra  order  of  magni¬ 
tude  to  be  added  to  that  accuracy. 

Because  the  resolution  of  the  gages  is  much  smaller  than 
the  variation  in  the  surface  texture,  the  variations  in  the 
readings  over  each  of  the  four  inch  segments  of  pavement 
reflect  the  coarseness  of  the  surface  of  the  pavement  at 
that  location.  A  large  variation  indicates  a  coarse  surface, 
while  a  small  variation  indicates  a  smoother  one.  A  useful 
measure  of  the  variation  in  the  readings  is  the  statistical 
variance  <V>  - 


V 


<Xx  -  X)a 

—  y71  - * - -  — 

<T>  -  l> 
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wh  ere 


X.  —  measurement 


X  —  mean  of  n  measurements 


photodetector 


Figure  III-2  The  dropout  case 


n  =  nuffitior  of  measurimorito 

2.  2  p  e  e  d  and  Distarat'  Calculation 

The  speed  and  location  of  the  vehicle  during  a  test  are 
recorded  with  each  pvof i le/def lection  measurement.  By  know¬ 
ing  the  frequency  of  the  readings  <16  kHz)/  the  number  of 
readings  (n)>  and  the  distance  over  which  the  readings  were 
taken  (four  inches)#  the  speed  can  be  calculated  from 

speed  =  (4)(160C0)/n  inches/second  (7) 

An  alternate  method  also  was  used  to  calculate  the 
speed.  The  time  to  traverse  four  inches  was  measured  by  a 
clock  in  the  computer  which  ran  the  data  collection  program. 
By  dividing  the  this  time  into  four  inches#  the  speed  was 
6  e v e r mined. 

Gene  readings  were  maue  uver  the  four  inch  distance  for 
every  ••pot  of  pavement  That  is#  the  gages  would  read  for 
four  inches*  skip  eight  inches#  read  for  another  four 
inches,  etc  As  the  readings  are  one  foot  apart#  the  number 
of  readings  gives  the  distance  travelled  in  feet.  Thus# 


locations  arc  secured. 
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C.  The  Eouiomen 


1.  The  Gage  Heads 


A  review  of  the  literature  revealed  that  there  were 
several  non-contact  prof i 1 ometer s  currently  in  use.  Most  of 
the  systems  used  gages  mounted  on  a  rigid  beam  towed  by  a 
vehicle.  The  gages  were  electronically  connected  to  and  read 
out  on  a  piece  of  recording  equipment.  In  order  to  know 
where  the  measurements  occurred,  a  distance  meter  was  con¬ 
nected  to  the  vehicle. 


Sound,  visible  and  infrared  light  gages  (using  lasers 
and  polarized  light  sources)  have  all  been  tried  with  some 
modicum  of  success.  Having  seen  some  prior  success  with 
light  emitting  diodes,  for  static  pavement  deflection  meas¬ 
urement  {(Harr  and  Ng-A-Qui  (1977)>,  the  use  of  light  sen¬ 


sors  was  the  first  course  pursued  in  this  study.  Several 
arrangements  of  lenses,  mirrors,  and  optical  detectors  were 
considered,  including  one  particularly  promising  one  that 
utilized  the  Scheimflug  condition  {(Hallert  (1979)>.  The 
basic  arrangement  of  the  optics  for  all  such  schemes  is 
shown  of  Figure  III-3.  A  light  source  is  reflected  from  a 
pavement  onto  a  detector  that  is  calibrated  to  provide  its 
distance  above  the  point  of  reflection.  Unfortunately,  the 
detectors  used  (either  CCDs  or  linear  photodiodes)  were 


found  to  be  sensitive  to  the  intensity  of  light.  This  meant 


Figure  III-3 


Arrangement  of  optics 


70 


that  readings  were  sensitive  to  the  color  of  the  pavement 
and  to  the  intensity  of  sunlight. 

Several  commercial  companies  produce  non-contact  dis¬ 
tance  measuring  gages.  The  gage  produced  by  Selective  Elec¬ 
tronic  was  chosen  for  this  study.  Their  gage,  the  "Optoca— 
tor",  had  the  required  resolution,  range  and  standoff*.  In 
addition,  its  outputs  were  compatible  with  modern  computers. 
Thus,  the  most  important  component  of  the  system  was 
secured. 

2.  The  Fifth  WJheel 

As  was  noted  in  the  introduction  to  section  III,  it  was 
necessary  to  have  the  gages  read**  where  the  other  gages  had 
read  The  gages  read  con t inuous 1 y ;  however,  their  readings 
were  sampled  only  occasionally,  as  noted  above.  The  sam¬ 
pling  rate  was  controlled  by  the  fifth  wheel. 

Several  methods  of  timing  were  considered.  First,  the 
use  of  a  shaft  encoder  on  the  axle  of  the  load  vehicle  was 
contemplated  to  tell  when  the  required  distance  had  been 
traversed.  This  method  was  discarded  because  of  its  inflexi¬ 
bility  The  encoder  would  have  to  be  recalibrated  for  each 
vehicle  used,  causing 

*  Standoff  is  the  distance  from  the  gage  to  the  center  of 

its  measuring  range. 


> e  [ r.  t li e  statistical  sense 


a  substantial  increase  in  set-up  time  when  changing  vehi¬ 


cles. 

Next,  long  distance  laser  measuring  devices  were  con¬ 
sidered.  These  devices,  commonly  used  in  surveying,  are 
extremely  accurate,  and  have  ranges  on  the  order  of  a  mile. 
These  devices  were  eliminated  from  consideration  because 
their  response  time  u/as  too  slow  (0.6  Hz  vs.  about  7.0  Hz 
required).  Moreover,  this  method  would  have  required  the 
laser  to  be  coupled  with  a  tracker  so  it  could  follow  the 
load  vehicle.  The  tracker  would  add  expense  and  complication 
to  the  system. 

Finally,  it  was  decided  to  use  a  calibrated  wheel  of  the 
kind  used  to  measure  the  performance  of  four-wheeled  vehicles 
(hence  the  term  "fifth  wheel").  "Labeco"  manufactured  the 
fifth  wheel.  A  shaft  encoder,  manufactured  by  Madison  Elec¬ 
tric.  was  attached  to  the  wheel  to  give  a  pulse  every  four 
inches  (for  this  wheel.  21  pulses  per  revolution).  The  p ower 
supply  for  the  encoder  was  made  at  Purdue  University. 

An  error  analysis  for  the  performance  fifth  wheel  is 
given  in  Appendix  B,  1 

3.  The  Rigid  Beam 

For  the  profile  algorithm  described  in  the  second  part 
of  this  section  to  work,  the  gages  must  remain  fixed  in 
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space  relative  to  one  another.  As  it  was  the  purpose  of  this 
study  to  develop  a  system  that  could  be  used  with  almost  any 
load  vehicle<  fixing  the  gages  rigidly  onto  the  load  vehicle 
was  ruled  out.  Attaching  the  gages  to  a  rigid)  portable  beam 
was  thought  to  provide  the  necessary  flexibility. 

Considerations  in  choosing  the  beam  included 

1.  weight) 

2.  stiffness) 

3.  avai lab  i  1 i ty )  and 

4  cost 

Steel  and  aluminum  I-sections  were  considered.  Steel 
sections  are  too  heavy  and  difficult  to  work  with.  Aluminum 
sections  are  available  only  in  one  ton  lots.  Consequently 
it  was  necessary  to  fabricate  a  beam.  The  design  was  accom¬ 
plished  using  the  computer  program  SAP  IV  -C(Bathe.  Wilson, 
and  Peterson  (1974)>  at  Purdue  University.  A  truss  afforded 
the  best  stiffness  to  weight  ratio.  A  rectangular  truss  was 
designed,  as  it  permitted  easier  gage  mounting  and  afforded 
some  protection  for  the  gages.  The  truss  was  fabricated  at 
Purdue  University.  Figure  JII-4  shows  a  picture  of  the 
beam. 


The  lengtn  of  the  beam  is  of  considerable  impor tantance. 


The  beam  had  to  be  long  enough  to  mount  one  gage  near  the 


73 


Figure  III-4  Rigid  beam  with  laser  gages 


load  wheel,  and  to  hold  the  undeflected  profilometer  (gages 

A,  B  and  C)  outside  the  zone  of  influence  of  the  load  wheel 

Highter  and  Harr  (1973)  and  Baladi  and  Harr(1976)  give  data 

that  show  the  deflection  basin  for  flexible  pavements  rarely 

extends  beyond  five  feet  in  front  of  the  load  wheel,  in  the  direction 

of  motion.  A  ten  foot  beam  (nine  feet  between  gages)  was  chosen 

The  beam  had  to  be  stiff  enough  to  prevent  the  gages 
from  moving  out  of  their  original  relative  alignment.  The 
allowable  movement  was  taken  to  be  the  resolution  of  the 
gages.  Any  greater  movement  would  not  allow  one  to  determine 
if  the  pavement  had  deflected  or  whether  the  beam  had  moved 
vertically. 

4.  Data  Collection, Storage  and  Reduction 

a  Overview.  The  data  acquisition  system  collected  the  data 
generated  by  the  lasers  and  the  fifth  wheel.  It  consisted  of 
a  microcomputer  and  various  peripherals  The  computer  was 
triggered  by  the  pulses  from  the  fifth  wheel.  When  the  fifth 
wheel  gave  the  proper  pulse,  the  computer  clock  started,  and 
the  computer  began  'looking  at'  the  readings  from  the 
lasers  It  then  averaged  the  readings  from  each  gage  and 
computed  standard  deviations  for  all  of  them  When  the  fifth 
wheel  sent  the  next  pulse  (at  four  inches  of  travel),  the 
computer  stopped  looking  at  the  readings,  stopped  looking  at 
the  clack,  and  the  data  were  put  in  core  memory  When  the 
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entire  desired  pavement  section  had  been  traversedi  the 
operator  flipped  a  logic  switch  which  told  the  computer  to 
stop  the  data  gathering  program.  The  program  stopped  and  the 
data  were  written  onto  a  floppy  disc.  The  profile/deflection 
program  later  uses  the  data  from  this  disc 


b.  Th  e  Eou  i  ament. 

1.  Hardware.  The  electronics  for  the  data  collection, 
storage  and  reduction  system  consisted  of  - 

1.  a  Heath  H-ll  microcomputer, 

2.  a  Heath  H-19  video  terminal 

3.  a  Heath  H-27  dual  floppy  disc  drive 

4.  a  Texas  Instruments  TI-810  line  printer 

5.  a  Digital  Equipment  Corporation  BA11-NE  bus 
extender,  and 

6.  a  'logic  box'  with  four  switches. 

The  H-ll  microcomputer  is  based  on  the  DEC  LSI-il 
microprocessor.  Consequently,  the  H-ll  uses  DEC  assembly 
language.  The  H-ll  is  a  lo-bit  machine.  This  was  required 
because  the  Optocators  have  16-bit  parallel  output. 


76 


The  maximum  memory  capacity  was  56k  bytes. 

The  H-19  video  terminal  is  used  to  'talk'  to  the  com¬ 
puter.  It  uses  the  RS-232c  serial  interface,  and  can 
display  24  80-character  lines. 

The  H-27  dual  floppy  disc  drive  accepts  two  eight  inch 
discs.  It  uses  a  single-side,  single  density  IBM  3740  for¬ 
mat.  The  two  discs  provide  about  500k  bytes  of  storage. 

The  TI-810  line  printer  is  made  by  Texas  Instruments.  It 
accepts  14-inch  paper  <132  characters),  and  uses  the  RS— 232 c 
serial  interface.  The  printing  rate  is  150  cps. 

The  BA11— ME  bus  extender  is  used  to  extend  the  number  of 
parallel  input/output  <I/Q)  boards  the  computer  will  accept. 
The  parallel  I/O  boards  for  the  lasers  were  used  in  the 
BA  1 1 -NE. 

The  logic  box  was  used  to  send  one  bit  signals  to  the 
computer.  It  consisted  of  a  five  volt  DC  transformer  and 
four  no-bounce*  switches.  The  switch  output  was  either  zero 
or  five  volts,  depending  on  tne  position  of  the  switch  They 
were  used  to  control  the  data  gathering  program. 

All  of  the  above  equipment  used  110-volt  60  Hz  AC  power. 
►4 ,.  t  h  the  exception  cf  the  terminal  and  the  logit  box,  all 
wore  approximately  the  size  and  weight  of  a  modern  electric 

*  A  no— bounce  switch  is  one  whose  contacts  do  not  bounce 

upon  closure  Bouncing  switches  can  be  made, 

electronically,  to  act  like  no-bounce  switches 


typewriter.  The  bus  extender  was  purchased  from  Hamilton— 
Avnet  electronics.  The  logic  box  was  fabricated  at  Purdue. 
The  remaining  equipment  was  purchased  from  Heath  Corpora¬ 
tion. 

2.  Software.  Software  is  the  set  of  programs  that  is 
used  to  operate  a  computer.  For  this  study,  the  term  system 
software  will  be  used  to  describe  that  which  came  with  the 
computer.  User  software  refers  to  that  which  was  written 
specifically  for  this  study. 

System  software  included  a  FORTRAN  IV  compiler,  editor, 
linker,  library  functions,  assembly  language  compiler,  and  a 
file  manager.  User  software  included  a  data  gathering  pro¬ 
gram  (GATHER),  a  data  file  previewer  (LIST),  and  the  data 
reduction  program  (CALC).  The  usefulness  of  the  data  file 
p-r'Pviewer  is  described  in  Appendix  G. 

The  system  software  functions  will  be  described  by 
tracing  typical  steps  in  the  program  creation  and  execution 
process.  A  program  is  written  and  modified  in  the  editor. 
When  it  is  correct,  it  is  compiled  by  the  computer.  Then  it 
is  linked  to  the  library  files  by  the  linker.  If  no  errors 
are  found,  the  program  can  be  run.  Should  the  program  be 
written  in  assembly  language,  the  assembly  language  compiler 
is  substituted  for  the  FORTRAN  IV  compiler  The  file  manager 
is  used  to  view,  compress,  and  copy  files. 


The  data  gathering  program.  GATHER,  was  written  by  Mr. 
Cary  Cox  of  the  U  S.  Army  Corps  of  Engineers  Waterways 
Experiment  Station,  Vicksburg,  Mississippi.  It  was  written 
in  assembly  language  so  it  would  be  fast  enough  to  keep  up 
with  the  lasers.  When  GATHER  is  run,  it  asks  the  operator 
to  name  the  data  to  be  gathered.  After  accepting  this, 

GATHER  prompts  the  operator  to  flip  the  start  switch  on  the 
logic  box.  This  being  cone,  GATHER  monitors  the  fifth  wheel 
shaft  encoder,  the  logic  box,  the  line  time  clock  on  the 
computer,  and  the  four  laser  gages  The  fifth  wheel  encoder 
gives  a  pulse  every  four  inches  of  travel  over  the  pavement. 
GATHER  checks  the  pulse  frequently,  and  counts  the  pulses 
The  beginning  of  the  pulse  representing  the  first  four 
inches  of  every  foot  tells  GATHER  to  begin  reading,  counting 
and  averaging  the  readings  from  the  lasers.  At  the  same 
time,  it  begins  counting  1/60-ths  of  a  second  on  the  line 
time  clock.  It  also  begins  calculating  standard  deviations 
of  all  four  of  the  lasers.  Later,  the  time  is  used  to  calcu¬ 
late  speed  While  doing  this,  the  computer  checks  the  pulse 
to  see  if  four  inches  of  pavement  have  been  traversed.  When 
the  pulse  changes,  the  four  inches  have  passed,  and  GATHER 
puts  the  four  average  readings,  the  number  of  1/60-ths  of  a 
second,  and  the  four  standard  deviations  into  memory.  When 
the  operator  flips  the  start  logic  switch  to  stop.  GATHER 
puts  values  of  negative  one  in  for  all  the  readings,  stores 
them  in  memory,  and  writes  the  entire  data  set  onto  the  disc 
under  the  name  assigned  by  the  operator.  The  negative  ones 


signal  the  end  of  the  data  to  the  data  reduction  program. 


After  the  data  are  gathered,  the  operator  may  wish  to 
preview  the  data  before  reducing  them  Of  interest  are  the 
number  of  readings,  the  number  of  zero  readings,  the  stan¬ 
dard  deviations,  and  the  length  of  the  file.  The  user  pro¬ 
gram  LIST  enables  the  operator  to  do  this  e *  ami na 1 1  on. 

The  averaged  laser  readings  are  converted  to  profile  and 
deflection  measurements  by  the  user  program  CALC.  CALC  reads 
the  data  from  the  file  created  Dy  GATHER.  It  then  uses  all 
the  data  to  compute  two  constants  -  the  two  out— of— line 
measurements  for  gages  A  anfl  D  (Recall  that  for  the  algo¬ 
rithm  to  work,  the  gages  must  be  in  a  straight  line.  Gages  B 
and  C  are  used  to  define  a  straight  line  CALC  calculates 
the  deviation  of  gage  A  is  any  time,  and  that  for  gage  D  if 
the  pavement  is  undeflected  (see  Appendix  E.  ))  These  two 
constants  are  then  added  to  their  respective  calibration 
equations  to  correct  the  "out— of— lineness".  CALC  then 
rereads  the  data,  using  the  algorithm  described  in  section 
III,  B,  2,  to  obtain  the  profile  and  deflections.  Then,  if 
the  operator  desires,  plots  of  the  profile  vs.  distance  and 
deflection  vs.  distance  are  produced 

Occasionally,  null  readings  occur.  The  effect  of  a  null 
reading  on  the  algorithm  depends  on  when  and  at  which  gage 


it  occurs. 
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There  are  two  ways  to  create  a  null  reading.  First,  the 
pavement  may  move  out  of  range  of  the  gage  The  range  of 
each  gage  is  +  2.  5  inches  about  the  standoff,  which  is  ten 
inches.  The  gage  then  sends  a  null  reading  to  the  computer. 
Second,  the  gage  may  make  less  than  thirty  readings  in  the 
travelled  four  inch  span.  This  is  interpreted  by  the  program 
as  the  gage  having  gone  out  of  range  part  way  through  the 
four  inches.  CALC  also  interprets  this  as  a  null.  The  CALC 
program  is  sensitive  to  both  cases. 

A  special  condition  occurs  when  a  test  is  started  The 
first  seven  feet  of  travel  (the  distance  between  gages  C  and 
D)  is  a  special  zone  at  this  time  No  deflections  are  meas¬ 
ured  there,  as  the  undeflected  profile  has  not  been  measured 
in  that  region.  Consequently,  if  the  average  reading  of  gage 
D  is  null  there,  it  has  no  consequence.  If  gages  A,  B,  or  C 
produce  a  null  (or  any  combination  of  nulls)  in  this  rone, 
CALC  ignores  past  data  and  Degins  the  algorithm  afresh. 

This  happens  because  the  algorithm  cannot  tolerate  gaps  in 
the  continuity  of  the  incoming  data. 

After  the  first  seven  feet  are  traversed,  deflections 
are  calculated  and  the  profile  continues  to  be  calculated. 

I?  gage  A  produces  an  average  reading  that  is  a  null,  a 
deflection  is  calculated,  and  the  profile  is  ended.  At  the 
next  reading,  a  new  profile  is  begun  and  six  deflections  are 
lost  (those  between  gages  C  and  D>  The  algorithm  begins 
anew  as  described  in  the  pre-'icus  paragraph  If  gages  3  or 
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C.  or  both,  produce  an  average  reading  that  is  a  null,  no 
deflection  is  calculated,  the  profile  ends,  and  then  begins 
anew  at  the  next  reading.  For  this  case,  seven  deflections 
are  lost.  The  CALC  program  performs  all  these  operations 
automat i ca 1 1 y  . 

The  CALC  program  puts  out  the  following  information  - 

1.  the  distance  travelled, 

2.  the  number  of  null  readings  at  each  gage, 

3.  the  two  out-of-line  constants, 

4.  the  undeflected  profile, 

5.  the  deflections,  their  mean  and  standard  deviation, 

6.  a  measure  of  the  texture,  and 

7.  plots  of  the  deflected  profile  or  deflections  versus 
d  istance. 


All  of  the  user  programs  are  listed  in  Appendix  C 
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IV.  RESULTS  AND  DISCUSSION 


A.  Profile  and  Deflection 

The  developed  apparatus,  described  in  Section  III/  was 
tested  at  the  Waterways  Experiment  Station  (WES)  of  the  U.  S. 

Army  Corps  of  Engineers  in  Vicksburg,  Mississippi.  This  was 
done  in  December  1981  The  first  tests  were  run  with  the  beam 
mounted  on  a  load  cart,  simulating  one  main  gear  on  an  F-4  air¬ 
craft.  The  computer  was  carried  in  a  station  wagon  which  towed 
an  electric  generator.  These  tests  uncovered  many  hardware 
problems  —  particularly  the  need  to  use  no-bounce  switches  in 
the  logic  box.  Several  bugs  in  the  software  were  also  discovered. 

The  next  tests  were  run  with  the  beam  mounted  on  the 
side  of  a  semi— trailer  at  the  Waterways  Experiment  Station. 

This  was  done  in  January  1982  Gage  D  (Figure  III— 1)  was 
adjacent  to  the  rear  wheel.  The  trailer  carried  the  computer 
and  the  electric  generator.  The  profile  calculated  from 
these  tests  plunged  with  distarce  (although  the  actual  pro¬ 
file  was  nearly  level).  The  magnitude  of  the  plunge  was  as 
great  as  fifty  feet  over  a  one  hundred  foot  horizontal  dis¬ 
tance.  It  was  exponential  -  meaning  that  the  amount  of 
plunge  increased  rapidly  with  distance  traveled  It  was 
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obvious  that  there  was  a  serious  problem  with  the  system. 

The  calculated  deflections  from  these  tests  were  also 
wholly  unrealistic.  The  weight  of  the  semi-trailer,  coupled 
with  the  pavement  structure,  should  not  have  produced 
deflections  greater  than  about  one  half  inch. 

After  more  testing,  and  a  review  of  the  algorithm,  a 
solution  to  the  plunging  profile  problem  was  found.  It  was 
discovered  that  the  gages  had  to  be  aligned  extremely  care¬ 
fully  on  the  rigid  beam.  Their  relative  positions  on  the 
beam  must  be  known  with  great  precision  A  slight  error  in 
their  alignment  causes  a  large  error  in  the  calculated  pro¬ 
file.  This  error,  perhaps  small  at  first,  accumulates 
rapidly  -  each  error  being  added  to  the  previous  error.  This 
error  manifested  itself  in  causing  the  measured  profile  to 
p 1 ung e  rap i d  1  y . 

A  source  of  the  noted  error  will  be  illustrated  by 
example.  Consider  the  following  case  -  the  laser  gages  lie 
on  a  straight  line  on  the  rigid  beam,  the  beam  moves  paral¬ 
lel  to  a  perfectly  flat  pavement,  and  the  datum  is  parallel 
to  the  pavement.  The  equation  that  calculates  the  next  pro¬ 
file  point  is  - 


AA  =C+CC-A-2  (  C+CC-B-  BB  )  (1) 


This  equation  was  derived  in  section  III,  equation  (4)  For 
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the  conditions  just  stated,  equation  (1)  will  calculate  a 
perfectly  flat  surface.  If  gage  A  is  out  of  vertical  align¬ 
ment  with  respect  to  gages  B  and  C  (Figure  1 1 1  —  1 )  by#  say 
one  inch  at  the  first  measurement,  an  error  in  A  of  one  inch 
will  occur.  AA.  then,  will  be  too  small  by  one  inch.  Recall 
that,  in  the  next  iteration,  AA  will  become  BB.  and  the  pre¬ 
vious  BB  will  became  the  present  CC.  So.  in  the  next  itera¬ 
tion.  A  is  again  one  inch  too  large,  and  BB  is  one  inch  too 
small.  The  net  error  is  that  after  the  second  measurement. 

AA  will  be  too  small  by  three  inches.  The  next  time.  AA  will 
be  too  small  by  six  inches.  The  relation  between  any  two 
consecutive  errors  is  — 


iteration  numbe 


Figure  IV- 
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the  relative  vertical  positions  of  the  gages  on  the  rigid 
beam.  It  did  this  by  using  the  measurements  from  the  gages 
themselves,  thus  having  the  same  accuracy  as  the  entire 
system.  The  result  of  this  scheme  was  to  produce  the  "out- 
of-line"  constants  -  the  deviation  from  a  straight  line  of 
any  three  gages. 

In  order  to  get  the  out- of -line  constant  for  a  gage, 
the  beam  had  to  be  used  where  no  measurable  deflection  would 
occur  At  first,  the  beam  was  hung  from  cantilever  beams  on 
the  F-4  load  cart,  far  away  from  the  zone  of  influence  of 
the  loaded  wheel  Unfortunately,  the  spring  mounts  used  to 
attach  the  rigid  beam  to  the  vehicle  permitted  too  much 
lateral  movement.  Moreover,  it  was  suspected  that  the  gages 
moved  out  of  their  range  of  reliable  measurement  during 
the  test.  A  special  lightweight  iricycle  was  then 
constructed.  The  beam  was  attached  to  its  undercarriage 
in  such  a  way  that  lateral  movements  were  greatly 
reduced. 

Using  the  tricycle,  new  out-of-line  constants  were  cal¬ 
culated  for  the  gages.  Once  these  were  obtained,  the  beam 
was  removed  and  transferred  to  the  semi-trailer  and  addi¬ 
tional  measurements  were  made.  The  results  from  these  tests 
indicated  that  the  calculated  profile  wandered  only  slightly 
from  the  actual  (surveyed)  profile. 
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Figure  IV -2  provides  the  comparison  between  the  transit 
surveyed  profile  and  three  laser  beam  calculations  of  the 
same  profile.  This  and  all  succeeding  transit  surveys  were 
made  at  one  foot  intervals  The  trend  of  che  calculated 
profiles  is  the  same  as  that  of  the  surveyed  profile.  The 
variance  of  the  endpoints  is  seen  to  be  only  a  few  inches. 
Some  of  the  variance  can  be  removed  by  forcing  the  endpoints 
of  the  laser  surveyed  profile  through  the  endpoints  of  the 
transit  survey  (recall  that  BB  and  CC  were  arbitrarily 
chosen  initially).  This  same  technique  is  practiced  at  the 
Transport  and ’ Road  Research  Laboratory  The  result  of  this 
technique  is  shown  in  Figure  IV-3.  The  ADJUST  program  (in 
Appendix  C>  was  used  to  do  this  The  deviations  were 
reduced  to  less  than  an  inch.  However,  the  deviations  them¬ 
selves  still  varied.  Ower  a  distance  of  100  feet,  the 
deviations  ranged  fi  ..  i  positive  three  tenths  inch  to  nega¬ 
tive  three  tenths  inch,  which  wa s  still  considered  to  be 
excessive  Testing  with  a  Benkelman  beam  showed  this  to  be 
the  case  The  standard  deviation  of  the  deflections  was  on 
the  order  of  one  tenth  inch. 

Several  extended  profile  tests  were  run.  Three  five 
hundred  foot  sections  were  surveyed  with  both  the  present 
la-er  system  and  a  transit  All  the  transit  surveys  were 
mule  on  points  one  foot  apart  The  results  of  some  of  the 
tests  are  shown  in  Figures  IV  -1  through  Figure  IV-6.  The 
profiles  in  these  have  been  adjusted  so  that  the  endpoints 
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pass  through  the  same  elevation.  The  maximum  deviation  of 
the  laser  system  profile  from  the  transit  surveyed  profile 
varied  from  0.290  feet  to  1.205  feet,  over  the  several  pro¬ 
files. 

The  search  for  the  reason  for  the  lack  of  agreement 
between  the  transit  and  laser  surveyed  profiles  continued 
The  shape  of  the  laser  system  profile  changes  radically  with 
small  changes  in  the  out-of-line  constant.  The  profiles  in 
Figure  IV— 4  through  Figure  IV-6  use  an  out-of-line  constant 
chosen  to  minimize  the  maximum  difference  between  the  tran¬ 
sit  and  laser  surveys.  The  constants  were  chosen  to  four 
decimal  places.  This  is  one  place  beyond  the  accuracy  of 
the  laser  gages.  Closer  agreement  could  have  been  accom¬ 
plished  if  more  decimal  places  had  been  taken,  however,  this 
could  not  be  justified.  Figure  IV-6  and  IV-7  represent  the 
same  transit  survey.  The  out-of-line  constant  in  figure 
IV-7  is  0.0277  inches,  while  that  used  in  Figure  IV-6  is 
0.0267  inches.  The  maximum  deviation  from  the  transit  pro¬ 
file  in  Figure  IV-6  is  1.205  feet,  while  the  maximum  devia¬ 
tion  in  Figure  IV-7  is  5.869  feet. 

The  raw  data  from  the  long  profile  tests  contained 
inherent  errors.  Over  each  five  hundred  foot  length,  the 
laser  system  recorded  about  five  hundred  and  ten  readings, 
instead  of  the  expected  five  hundred  and  one  readings.  The 
itandarl  deviations  of  the  gj.;>.i.  usually  between  zero  and 
fifty,  sometimes  jumped  an  o’-der  of  magnitude  Both  the 


extra  distance  and  the  large  standard  deviations  were  con¬ 
sidered  indications  that  the  system  was  still  not  operating 
properly 

Testing  with  the  developed  system  has  only  just  begun 
and  improvements  are  to  be  made.  Suggestions  For  these  will 
be  addressed  below.  At  this  point  it  should  be  noted  that 
the  main  advantage  of  the  system  (its  ability  to  do  frequent 
and  rapid  testing)  is  worthy  of  the  effort  necessary  to 
achieve  these  ends.  As  more  testing  is  done,  more  refine¬ 
ments  will  come  to  light.  The  ability  to  test  pavements  fre¬ 
quently  and  rapidly  will  provide  the  pavement  engineer  with 
a  quantitative  measure  of  the  performance  of  his  pavement 
system. 

The  system  described  herein  has  no  built-in  redundancy. 
Some  form  of  redundancy  may  make  it  easier  to  debug  the  sys¬ 
tem.  The  addition  of  more  laser  gages  or  of  a  device  that 
provides  some  measure  of  the  tilt  of  the  beam  (such  as  a 
tiltmeter)  will  aid  in  accomplishing  this. 

The  addition  of  more  laser  gages  to  the  present  sytsem 
will  provide  a  redundancy  If  gages  are  added  between  gages 
C  and  D  (Figure  1 1 1  —  1 ) .  and  these  gages  are  out  of  the  lone 
of  influence  of  the  load  wheel,  they  can  be  used  to  calcu¬ 
late  the  undeflected  profile  in  combination  with  any  other 
two  gages.  These  additional  calculations  of  the  profile 
could  be  compared  with  the  first  calculation  (and  with  each 
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other)  to  help  determine  if  the  system  was  working  properly. 
The  same  method  that  was  used  to  align  gages  A.  B  and  C 
(Figure  1 1 1  —  1 )  could  be  used  to  align  any  additional  gages 
(see  Appendix  E).  If  many  more  gages  are  added,  it  may  be 
necessary  to  lengthen  the  rigid  beam  in  order  to  keep  the 
gages  out  of  the  zone  of  influence  of  the  load  wheel.  Cal¬ 
culating  several  profiles  s imu 1 taneous 1 y  allows  the  algo¬ 
rithm  to  "bridge"  over  gaps  in  the  data.  With  the  present 
system,  if  one  of  the  gages  goes  out  of  range,  the  profile 
skips  a  point,  and  then  begins  anew,  relative  to  a  new 
datum.  It  should  not  be  expected  that  the  new  and  old  data 
will  be  coinc  ident. 

If  oroperly  placed,  the  additional  gages  would  also 
provide  better  definition  of  the  deflection  basin.  The 
deflections  at  these  additional  gages  could  be  calculated  in 
a  way  similiar  to  that  of  gage  D  (Figure  1 1 1  —  1  ) .  The  equa¬ 
tion  used  to  calculate  the  deflection  at  gage  D  (equation 
(5)  in  Section  III)  can  be  changed  by  modifying  "the  ratio 
of  the  distance  between  gages  B  and  C"  to  "the  distance 
between  gages  C  and  the  additional  gage" 

A  tiltmeter  (or  any  device  that  measures  the  tilt  of 
the  beam'  could  also  be  used  to  provide  a  redundancy.  The 
tiltmeter,  in  conjunction  with  two  laser  gages,  can  be  used 
to  calculate  the  undeflected  or  deflected  profiles  The 
methodology  is  detailed  in  Appendix  A.  As  was  the  case  with 
additional  gages,  another  calculation  of  the  profile  is 
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made.  These  profiles  could  then  be  compared.  Adding  a 
tiltmeter  to  the  present  system  will  enable  the  system  to 
"bridge"  lost  readings,  as  does  the  addition  of  laser  gages. 
The  addition  of  a  tiltmeter  mould  not  require  the  lengthen¬ 
ing  of  the  rigid  beam.  An  extremely  accurate  and  responsive 
tiltmeter  mould  be  required. 

At  present,  the  laser  system  uses  the  first  four  inches 
of  travel  per  foot  to  record  data.  This  mas  chosen  because 
it  mas  thought  that  the  time  required  to  travel  the  other 
eight  inches  per  foot  mould  be  needed  to  do  calculations  and 
store  the  data.  However,  the  present  data  collection  program 
is  fast  enough  to  double  the  rate  previously  anticipated. 
This  means  that  the  present  system  can  be  used  at  double  its 
present  rate  of  output.  The  system  can  be  used  to  record 
readings  over  the  first  tmo 

-  four  inch  sections  of  every  foot  of  travel.  The  last 
four  inches  of  travel  time  (so  to  speak)  are  sufficient  to 
do  the  necessary  calculations  and  to  store  the  data  The 
result  is  that  tmo  profile  and  tmo  deflection  points  can  be 
obtained  for  every  foot  of  travel,  instead  of  the  present 
single  profile  and  one  deflection  point.  Clearly,  this  is  a 
function  of  horn  fast  the  rigid  beam  is  moving.  At  the 
present  "creep"  speed,  the  output  can  be  doubled  The  upper 
bound  for  the  speed  to  do  this  remains  to  be  found. 

Being  able  to  gather  pavement  deflection  data  fre¬ 
quently  and  rapidly  yields  the  advantage  of  being  able  to 
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utilize  transfer  function  theory  to  a  high  degree  to  predict 
pavement  perf  or man  c  e.  Given  the  deflection  caused  by  one 

vehicle  and  the  transfer  function  of  the  pavement  system/ 
the  deflections  for  other  anticipated  vehicles  can  be 
predicted.  The  transfer  functions  change  with  time  and  sea¬ 
son.  The  laser  system  allows  frequent  updating  of  the 
transfer  function  and  hence  renders  it  as  a  valuable  tool 
(Boyer  and  Harr-  1973). 

It  should  be  emphasized  that  the  present  system  offers 
the  ability  to  survey  a  pavement  in  a  global  sense.  No  other 
device-  at  the  present  time,  does  so.  In  addition-  its  speed 
enables  it  to  be  deployed  in  such  a  manner  as  to  minimize 
the  interruption  of  traffic 

B .  Texture 


As  noted  in  Section  III,  S,  1,  texture  measurements 
were  made  by  recording  the  standard  deviations  of  the  read¬ 
ings  each  gage  made  over  the  short  distance 

Table  IV-1  shows  partial  sets  of  uncalibrated  standard 
deviation  data  for  the  different  pavements.  As  expected, 
the  standard  deviations  for  the  concrete  pavement  are 
smaller  than  those  for  the  asphalt  pavements  This  reflects 


the  fact  that  the  concrete  pavement  is  smoother 


TABLE  IV-1 


riple  VolutJ 


of  un calibrated  data 


gage 

numb  er 

of 

standard 

time 

T*  £■  c 

dire 

C. 

de  vi  a 

t  i  on 

1/60 

sec 

A 

B 

C 

A 

B 

c 

A 

B 

C 

c  or  ere 

te 

i 

23  5E 

2853 

r>  t 

V.  i  t_ 

•;.<  1 2 

”1  n 

X  o 

8 

72 

2 

12 

1 

L  ^  c. 

25:6 

291 2. 

c.  WC. 

28  c.’ 

232 

14 

17. 

54 

11. 

pr  t 

2  C  32 

2908 

265 

265 

266 

B 

1 1 . 

B. 

10. 

2943 

n  rt  i 
lOOc. 

2995. 

246 

24* 

247 

17 

11. 

6. 

10. 

1 

2930. 

2322. 

2912. 

248 

243 

249 

17 

19. 

7. 

10. 

. 

Loot 

2  a  i  0 . 

2 G  i'  cl. 

ao5 

265 

266 

e 

7. 

11. 

10. 

2870. 

2798. 

2980 

24° 

249 

24  9 

B 

6. 

12 

10. 

2375. 

2767. 

2899. 

241 

241 

242 

11 

7 

o 

9. 

escr.alt 


/>  ^  /  v_J 

1826. 

r'  r*  »-> 

a  j. 

1986. 

1922. 

2032. 

241 

1  8°0 

1901 

1985 

263 

1947. 

1  888. 

1977 

239 

1 Q97. 

1905 

2006 

234 

2041 

1947. 

2026. 

230 

14-64 

1962. 

2C 1 6. 

22? 

■?  ~  Z~.  C3 

1  'Jl'o.'c- 

i  CC-j 

A  -  W 

CO 

r  c  i;  s 

frier 

ion  surface 

4  r>r- 1 

>  w  ■  *  . 

«  n  n  ■» 

d  1  '*  G  . 

wl  W  < 

1909 

1 3BS 

2052. 

282 

1655 

1 651 

2013. 

236 

1 861 

2030 

217 

a  n  j  . 

1  C  '7 

>  •r-J  *"*  ' 

r£  \i  2  4 

206 

4  c  2:n 

1  a~'V 

2036 

139 

lAj- 

1  904 

2032 

212 

1  “22 

1  9  1 0 

2030 

134 

252 

52 

13. 

16. 

10. 

12 

241 

241 

13. 

17. 

13. 

12 

263 

268 

15. 

17. 

13. 

13 

229 

238 

14. 

20. 

20. 

1 2 

234 

234 

11. 

13. 

15. 

12 

231 

231 

136. 

10. 

13. 

12 

o"»c 

C-  7 

eLc.  7 

17. 

19. 

134 

12 

249 

1  2. 

24 

12 

w  /  w 

379 

46. 

54. 

20. 

15 

279 

281 

31. 

47. 

23. 

11 

244 

245 

34. 

36. 

40. 

9 

£26 

226 

30 

36 

31 

9 

•  nr 

y  '  ... 

208 

40. 

38. 

200. 

a 

212 

23 

44. 

34. 

a 

^  i  A 

C  -  A  w‘ 

21  1 

26. 

36. 

143 

9 

i  —  r> 

>  *  — . 

:  ?8 

25, 

29 

30. 

a 

IOC 


No  correlation  with  standard  texture  measurements  were 
made  Table  IV-1  shows  the  standard  deviations  measured 
over  three  pavements  exhibiting  different  textures.  The  con¬ 
crete  pavement  was  smoother  than  either  of  the  asphalt  pave¬ 
ments.  This  is  demonstrated  by  the  magnitude  of  the  stan¬ 
dard  deviations  listed  in  Table  IV-1.  The  concrete 
pavement's  standard  deviation  is  less  than  the  asphalt 
pavement's  standard  deviation.  It  is  expected  that  this 
measure  of  pavement  texture,  standard  deviation,  might  be 
correlated  with  other  standard  texture  measurements. 

C.  System  Ref i nement 

As  was  noted  above,  the  following  suggestions  are 
offered  to  refine  and  improve  the  developed  system. 

Some  advantage  can  be  had  by  using  a  different  computer 
system.  The  computer,  a  Heath  H-ll,  is  no  longer  offered  by 
the  manufacturer.  Buying  a  "military  specification"  com¬ 
puter  from  a  well  established  manufacturer  would  yield  two 
benefits  -  reliability  and  service.  A  DEC  PDP-11  is  recom¬ 
mended.  The  Heath  H-27  dual  disc  drive  used  a  single  side, 
single  density  format.  A  Winchester  hard/floppy  disc  drive 
could  collect  a  lot  more  information.  The  Heath  operating 
system  editor  was  difficult  to  use.  A  better  editor  would 
make  the  programmers  job  easier  The  TI-BIO  printer  was 
sufficient,  but  a  printer-plotter  would  produce  better  look¬ 


ing  plots.  A  Printronix 


printer-plotter  is  recommended 


r 
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Replacing  the  laser  gages  with  Polaroid  sonic  dis¬ 
tance  measuring  gages  warrants  further  examination.  Prelim¬ 
inary  tests  showed  the  sonic  gages  are  not  affected  by  winds 
up  to  fifty  miles  per  hour.  Their  resolution  was  determined 
to  be  0  01  inch.  The  gag'es  can  be  made  to  read  as  fast  as 
44  Hz  with  the  addition  of  a  precision  clock.  It  has  not 
been  determined  if  the  gages  are  sensitive  to  ambient  noise, 
the  temperature  of  the  pavement,  humidity  or  air  tempera¬ 
ture.  The  replacement  of  the  laser  gages  with  these  sonic 
gages  would  realize  a  considerable  savings  in  money  (1982). 
These  sonic  gages  would  always  measure  the  shortest  distance 
to  the  ground,  reducing  the  error  described  in  Appendix  B 
(relating  to  rotating  gages). 

The  gage  readings  taken  over  the  four  inch  distance  are 
averaged  and  recorded.  At  present  operating  speeds,  over 
1600  readings  are  averaged.  This  average  is  truncated  to  a 
four  digit  integer,  the  same  number  of  digits  that  the 
lasers  output.  Because  a  large  number  of  samples  is  taken, 
another  digit  can  be  accepted.  It  is  felt  that  this  extra 
digit  would  improve  the  accuracy  of  the  system.  In  lieu  of 
this,  the  average  could  be  rounded  off  instead  of  truncated. 

Presently,  the  system  records  more  readings  than  it 
should  That  is,  in'  five  hundred  feet  of  travel,  it  records 
about  five  hundred  and  ten  readings  (instead  of  five  hundred 
and  one).  It  is  felt  that  slack  in  the  connection  between 


the  shaft  encoder  and  the  fifth  wheel  is  causing  this 
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problem.  This  slack  should  be  removed.  As  it  is.  extra 
readings  might  result  if  the  encoder  rotates  relative  to 
the  frame  of  the  fifth  wheel 

The  large  standard  deviations,  noted  in  section  IV.  B, 
might  be  an  indication  that  the  laser  gages  are  defective 
At  one  point  in  the  study,  one  of  the  gages  failed  com¬ 
pletely.  Before  it  did  so.  it  gave  large  standard  deviations 
intermittently  and  then  constantly.  It  is  suggested  that 
the  resistor  that  failed  on  the  failed  gage  be  replaced  on 
the  other  gages. 


V.  SUMMARY  AND  CONCLUSIONS 


This  study  has  presented  a  method  whereby  the  texture  and 
profile  of  a  pavement  may  be  determined  from  a  moving  vehicle. 
The  system  is  relatively  portable  and  requires  little  contact 
(a  transporting  vehicle  and  a  fifth  wheel)  with  the  pavement 

The  profile  algorithm  has  been  derived  c,nd  refined  so 
that  a  fair  agreement  with  surveyed  profiles  has  been 
obtained  Important  advances  in  the  algorithm  pertaining  to 
the  gage  alignment  have  been  discovered  and  reviewed.  The 
profile,  which  is  used  in  roughness  studies,  can  now  be 
accomplished  reasonably  quickly,  not  only  in  the  sense  that 
the  test  is  quick,  but  also  inasmuch  as  the  data  reduction 
is  also  rapid. 

The  laser  system  presented  can  also  measure  texture. 

While  the  values  presented  have  not  been  correlated  with 
other  methods,  the  present  system  allows  highway  and  airport 
personnel  to  easily  and  rapidly  detect  c  h  a  n  q  e  s  in  their  tex¬ 
ture  of  the  pavements. 

A  method  to  measure  the  deflection  caused  by  a  moving 
vehicle  from  a  moving  vehicle  has  been  developed.  When 
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perfected,  this  system  mill  allow  the  pavement  engineer  not 
only  to  measure  deflections,  but  also  to  determine  the  remain¬ 
ing  life  of  a  pavement  system  (based  on  the  deflection). 

The  present  study  has  developed  a  working  prototype.  Some 


refining  of  the  system  mill  result  in  its  improvement. 
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VI.  SUGGESTIONS  FOR  FURTHER  RESEARCH 


A  method  to  obtain  rapid  non-contact  measurement  of 
deflections  from  a  moving  vehicle  was  developed  in  the 
present  study  This  section  lists  various  directions  that 
future  research  could  take. 

1  Concrete  pavements  could  be  tested  for  pumping 

potential  and  groove  depth  Pumping  requires  the 
movement  of  the  edge  of  a  slab  By  starting  the 
undeflected  profile  just  before  the  edge  of  the 
slab,  the  deflection  caused  when  the  load  reaches 
that  point  can  be  measured  The  procedure  might  be 

a  move  the  load  vehicle  into  position 
b  slowly  turn  the  fifth  wheel  by  hand  until 
the  first  reading  is  detected,  and 
c  move  the  load  vehicle  across  the  slab  edge 

It  would  be  helpful  to  have  a  visible  indicator 
of  where  the  first  gage  was  reading.  A  visible 
light  spot  from  a  HeNe  laser  could  be  used.  The 
above  described  method  method  would  allow  the  use 
of  the  existing  program  GATHER 
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To  measure  groove  depth,  a  new  program  would  have 
to  be  written.  This  one  would  not  calculate  a 
profile  It  would  monitor  the  fifth  wheel  and 
take  readings  from  the  first  gage  at  the  signal 
from  the  operator.  As  the  gage  approached  a 
groove,  the  operator  signals  the  computer  to  take 
readings  until  the  gage  has  passed  the  groove. 

If  the  beam  does  not  move  much  vertically  during 
the  test,  a  profile  of  the  groove  can  be  measured. 
Again,  a  visible  light,  indicating  the  position 
of  the  first  gage,  would  be  of  assistance  to  the 
operator. 

This  test  would  have  to  be  conducted  slowly,  to 
keep  the  rigid  beam  from  moving  vertically.  The 
operator  could  use  a  remote  switch  to  trigger  the 
the  computer,  so  he  could  walk  alongside  the 
gage 

2.  Running  the  system  transverse  to  the  runways  can 
evaluate  rutting  The  system  can  also  be  adapted 
to  railroad  cars  to  find  soft  spots  in  the  roadbed. 
The  steel  wheels  of  railroad  cars  could  be  used 
to  replace  the  fifth  wheel. 

3  Replacing  the  laser  gages  with  the  Polaroid  sonic 
distance  measuring  gages  might  be  explored.  This 
replacement  would  realize  a  considerable  savings 


in  money  (1982).  In  addition.  the  sonic  gages 
measure  the  shortest  distance  to  the  ground, 
reducing  the  error  described  in  Appendix  B 
(relating  to  rotating  gages).  With  their  suc¬ 
cessful  development,  a  measuring  system  would  be 
affordable  to  all  engineering  organizations 
interested  in  pavement  performance. 
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APPENDIX  A 


Ai terns  t e 
and 


Theory  for  Profile-  Texture 
Deflection  Measurement 


As  an  alternative  to  using  three  gages  for  the  profilom- 
eter,  tmo  gaqes  and  a  tiltm?ter  can  be  used  (presuming  one 
with  sufficient  accuracy  can  be  located)  A  til ‘meter  is  an 
electro-mechanical  device  that  measures  the  angle  of  a  sur¬ 
face  with  respect  to  a  level  datum 

The  algorithm  for  finding  the  profile  is  similar  to  that 
used  with  three  gages.  Consider  Figure  1 1  I  —  1  again.  This 
time-  gages  B  and  C  w  i  i  1  be  used.  The  inclination  of  the 
beam  is  8.  Initial  ly,  B ,  C,  and  CC  are  known  CC  mas  chosen 
to  define  the  datum  -  the  datum  being  level  and  passing 
through  the  end  of  CC,  with  x  known. 

b  =  x  sin  8  C 1  > 


B  +  b  +  BB  =  C  t  CC 


(2) 
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so> 


ED  =  C  +  CC  -  Q  -  b  (3) 


substituting  (1)  into  (3)  yields 


BE  =  C  +  CC  —  B  —  Cx  sin  9  )  <4> 

Again#  note  that  the  terms  on  the  right  side  of  the  equation 
are  known.  Thus*  the  next  pavement  profile  point  is 
obtained.  The  beam  moves  forward-  as  before,  until  gage  C 
is  over  where  gage  B  was,  and  the  process  is  repeated.  The 
method  can  be  used  with  gage  C  and  D  to  get  the  deflected 
profile. 

The  deflection  and  texture  measurements  can  be  calcu¬ 
lated  in  the  same  manner  described  an  Section  III,  B,  1. 
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APPENDIX  B 

Analysis  of  the  Theory 


1.  Fifth  wheel 

An  analysis  was  carried  out  to  see  the  effect  of  the 
fifth  wheel  slipping  on  the  pavement  If  the  wheel  slips- 
the  succeeding  gages  would  not  read  where  the  preceding 
gages  had  read.  If  the  vehicle  was  travelling  on  a  slope,  a 
false  elevation  would  be  computed.  This  was  a  function  of 
the  amount  of  slippage  and  the  angle  of  the  pavement.  Furth¬ 
ermore,  it  was  a  cumulative  effect,  that  is,  once  the  wheel 
slipped,  it  shifted  the  elevations  at  ell  subsequent  points. 
If  it  slipped  again,  this  error  was  added  to  the  previous 
error  for  all  the  points  subsequent  to  that' si ippage,  and  so 
forth  Figure  A-l  is  a  plot  of  this  effect. 

Figure  A-l  was  obtained  from  the  geometry  of  the  system, 
shown  in  the  upper  part  of  the  figure.  Here  are  shown  the 
front  three  gages  of  the  beam  (without  the  load  vehicle) 
after  one  iteration  The  gages  have  moved  forward  one  incre¬ 
ment  (dashed  lines)  and  are  making  a  reading.  Normally,  they 
weld  reed  directly  over  whore  the  previous  gage  had  read. 
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Since  the  fifth  wheel  slipped,  however,  the  beam  moved  more 
tnan  the  standard  increment  before  the  computer  received  the 
signal  from  the  fifth  wheel  to  take  a  reading.  The  horizon¬ 
tal  component  of  the  distance  slipped  is  6 .  The  vertical 
component  is  S  tan  £  As  8  increases  the  vertical  (eleva- 
ticn)  error  increases.  As  the  amount  of  slippage  increases, 
the  elevation  error  also  increases.  This  is  shown  on  the 
plot  in  Figure  A— 1.  the  errors  in  elevation  are  cumulative 
because,  in  each  succeeding  iteration,  the  (incorrect) 
elevations  are  used  from  the  previous  iteration.  When  these 
elevations  are  larger  than  they  should  be,  the  resulting 
elevations  will  be  calculated  as  being  larger  than  they  are. 
Also,  the  deflections  will  be  calculated  at  the  wrong  loca¬ 
tions.  They  will  be  in  error  by  the  amount  C  tan  6  until  the 
deflection  gage  reaches  the  position  where  the  first  gage 
occupied  at  the  time  of  the  slippage. 

The  fifth  wheel  may  go  faster  than  the  load  vehicle 
(the  opposite  of  slippage).  When  this  occurs,  the  exact 
opposite  effect  as  that  of  slippage  is  introduced.  Instead 
of  an  elevation  error  being  added  to  the  profile,  it  is  sub¬ 
tracted 

Consider  the  fifth  wheel  again.  Recall  that  the  fifth 
wheel  tells  the  computer  when  the  succeeding  gages  are 
located  at  a  position  where  the  preceding  gages  had  been. 

If  the  pavement  is  level,  every  foot  traversed  on  the  pave¬ 
ment  is  one  foot  traversed  on  the  datum  Thus,  by  counting 
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the  number  of  readings  made<  one  can  locate  any  particular 
reading  in  space.  If/  however,  the  pavement  slooes  at  an 
angle  8  -from  the  datum  ^positive  or  negative},  then  the 
fifth  wheel  does  not  measure  one  foot  intervals  along  the 
datum.  Thus/  for  the  horizontal  distance  to  be  properly 
determined/  they  must-  be  multiplied  by  cos  8.  Theta  can 
easily  be  determined  from  any  two  successive  readings.  This 
also  assumes  the  slope  of  the  pavement  between  the  fifth 
wheel  and  the  other  end  of  the  beam  is  constant  -  not  an 
unreasonable  assumption.  If  this  error  is  not  corrected/  the 
calculated  profile  will  be  flatter  and  longer  than  it  actu¬ 
ally  is. 


Yet  another  problem  can  occur  with  the  fifth  wheel.  This 
concerns  timing  the  gages  so  that  each  one  reads  where  the 
preceding  gage  had  read*.  Figure  A-2 a  shows  the  beam  and 
the  fifth  wheel  moving  on  a  surface  that  is  parallel  to  the 
datum.  Note  that  the  beam  and  the  surface  are  parallel.  As 
the  fifth  wheel  moves  forward  one  increment/  each  gage  moves 
up  to  where  the  previous  gage  had  been.  In  Figure  A-2b/  how¬ 
ever/  when  the  beam  is  not  parallel  to  the  surface,  one 
increment  moved  by  the  fifth  wheel  does  not  place  the 
succeeding  gages  over  where  the  preceding  gage  was.  How 
important  is  this  error  ? 


The  size  of  the  error  is  a  function  of  the  slope  of  the 
surface,  the  slope  of  the  beam,  and  the  length  of  the  beam 


♦  In  thr  stat/stical  sen;,  e 


Assume  the  slope  is  5  degrees  (1  :  11)  and  the  slope  of  the 

beam  is  5°  +  2.7°  =  77°  The  angle  2.7°  was  chosen  as  it 
is  the  upper  limit  the  beam  can  pivot  about  the  end  gage 
before  the  other  end  gage  goes  out  of  range  (that  is,  2.  5 
inches  upward  for  the  Optocators).  The  error  that  occurs  is 
that  the  elevation  of  the  pavement  (above  the  datum)  is  cal¬ 
culated  incorrectly.  Figure  A-3  shows  the  details  of  the 
problem.  In  this  figure*  the  solid  beam  is  parallel  to  the 
slope  and  measures  the  pavement  profile  at  B  when  the  fifth 
wheel  has  moved  one  foot  up  the  slope.  If  the  beam  had 
tilted  just  before  the  the  reading  was  made,  the  reading 
would  have  occurred  at  A.  the  elevation  at  A  would  then  be 
computed  correctly,  but  attributed  to  the  horizontal  dis¬ 
tance  E.  The  difference  bet ween  these  two  points,  S  y,  is  a 
function  of  the  slope  of  the  surface  and  the  amount  the  beam 
has  tipped.  For  the  severe  conditions  just  stated,  6  u  is 
significant.  For  a  slope  of  0.7  degrees  (1  115)  and  a  beam 

tilt  of  3  4  degrees,  S  y  is  0.002  inches.  This  is  signifi¬ 
cant,  when  one  considers  that  deflections  are  measured  in 
0.001  inch  increments. 

The  solution  to  this  problem  is  to  keep  the  b earn  as 
parallel  as  possible  to  the  pavement. 


m 


direction  of  motion 


datum 
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2.  Rotating  Gages 

In  the  foregoing  figures/  whenever  the  beam  tipped/  the 
gages  continued  to  make  readings  perpendicular  to  the  datum. 
This  could  only  occur  if  the  gages  were  -free  to  rotate  on 
the  beam.  If  the  gages  were  fixed  rigidly  tc  the  beam/ 
incorrect  data  could  be  generated.  A  numerical  example  is 
given  in  Figure  A-4.  This  figure  shows  that  a  small  tilt  of 
the  beam  over  a  level  pavement  produces  an  error  in  the  hor¬ 
izontal  distance  measurement. 

This  is  not  to  say  that  this  error  would  not  occur  if 
the  gages  were  free  to  rotate  on  the  beam  -  it  would,  only 
the  magnitude  would  be  about  twenty  percent  of  that  of  the 
fixed  gage  case  (see  Figure  A-5).  Thus,  from  this  stand¬ 
point.  it  is  more  desirable  to  use  pivoting  gages. 

It  is  important  tc  note  that  the  gages  must  pivot  in 
unison  .  if  they  do  not.  their  random  movements  could  gen¬ 
erate  an  error  larger  than  that  caused  by  gages  fixed  to  the 
beam. 

Even  for  readings  over  a  four  inch  distance  on  the  pave¬ 
ment  (instead  of  single  readings),  the  pivoting  gage  system 
is  still  more  desirable  The  advantage  comes  from  the  fact 
that  one  must  measure  at  the  same  location  as  the  previous 
gage 
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The  horizontal  distance  error  for  the  -loating  gage  system 
was  absolutely  less  than  that  for  the  fixed  gage  system.  As 
a  percentage  of  the  four  inch  distance/  the  floating  gage 
system  is  closer  to  the  theoretical  spacing  than  the  fixed 
gage  system.  The  closer  the  gages  read  where  the  previous 
gages  had  read/  the  smaller  will  be  the  error  in  the  deflec¬ 
tion  calculation. 

If  the  fixed  gage  system  is  used  on  a  slope,  with  the 
beam  parallel  to  the  slope,  the  error  in  the  horizontal 
reading  will  be  the  same  as  in  the  pivoting  gage  case.  That 
is,  the  fifth  wheel  will  not  measure  distance  along  the 
datum  uni  *:=■?  the  cos  6  correction  is  made.  Should  the  beam 
with  fixed  gages  become  oblique  to  the  slope,  the  distance 
that  the  gage  laser  light  spot  will  be  off  target  will  be 
the  same  as  in  the  zero  slope  case.  The  same  problem  occurs 
with  approximately  the  same  amplitudes  for  the  pivoting  gage 
case.  This  problem  can  be  alleviated  (though  not  eliminated) 
by  using  the  average  of  readings  taken  over  a  short  dis¬ 
tance,  rather  than  single  readings.  This  averaging  was  done 
in  this  study  (section  III,  B,  1). 

The  complexity  of  motion  that  the  beam  can  go  through 
has  only  been  touched  upon  in  this  Appendix.  It  has  examined 
only  one  of  the  many  gyrations  that  the  beam  can  exhibit 
(the  pitch  mode).  The  errors  due  to  roll  and  yaw  remain  for 
future  investi gator s. 
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APPENDIX  C 

COMPUTER  PROGRAMS 


This  Appendix  lists  the  dets  gathering  program 
(GATHER),  the  data  reduction  program  fCALC),  the  data  pre¬ 
view  program  (LIST),  and  the  profile  adjusting  program 
ADJUST. 

GATHER  is  written  in  the  DEC  assembly  language.  This 
was  done  because  it  had  to  operate  very  quickly  It  gathered 
the  data  from  the  gap^-s.  monitored  the  fifth  wheel,  and 
counted  pulses  from  the  line  time  clock.  GATHER  writes  the 
data  onto  the  disk  in  binary. 

CALC  is  written  in  FORTRAN  IV.  It  reads  the  data  (con¬ 
verted  from  binary  to  ASCII  by  program  LIST)  and  computes 
the  undeflected  ar.d  deflected  profiles  and  the  deflections. 

P aa  data  points  are  noted,  as  are  the  standard  deviations  of 
the  readings  made  over  the  four  inches  of  pavements. 

LIST  is  a  data  previewing  program.  After  the  data  are 
gathered,  the  operator  may  wish  to  see  if  it  is  good  before 
reducing  it.  LIST  tells  the  operator  how  many  out— of— range 
points  there  are  If  there  are  too  many,  the  operator  may 
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wish  to  drive  ever  that  particular  stretch  of  pavement 
again.  LIST  is  also  used  to  translate  the  data  file  into 
ASCII.  In  this  form/  the  data  is  compatible  with  most  any 
computer. 

ADJUST  adjusts  the  endpoints  of  the  profile  to  user 
specified  elevations.  The  effect  is  to  change  the  datum  to  a 
user  defined  datum.  In  particular/  it  is  useful  to  adjust 
the  datum  to  be  horizontal. 


d 


non 
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Prc-grajri  GATHER,  use?,  to  collect  the  data. 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

C  DATA  ACQUISITION  PROGRAM  FOR  DAVE  ELTON  OF  PURDUE 

c 

C  ••  GATHER  YE  DATA  WHILE  YE  MAY  " 

C 

C 

C  LOGICAL  VALUES  RECEIVED  FROM  LASER  SUBROUTINE 

C  (WHEEL,  THREE  SWITCHES,  ?■  TIME-CUT-ERROR) 

C 

r> 

L 

C  LOG IC=1=TRUE  LOG  I C=G=F ALSE 

C 

C 

BYTE  SPACE 
DATA  SPACE/ 1H  / 

BYTE  I YES 
BYTE  IANS 
DATA  IYES/1HY/ 

BYTE  NAMEF(BO) 

DIMENSION  JTIM2 ( 2 ) 

DIMENSION  IBUF2 ( 0192 ) ,  IBUF3<8191 ),  I3UF4(B190),  IBUF5(81B9>. 
*v.I3UF6<9183)  ,  IBUF7  (  SI  67  ) ,  1BUFS(S186),  IBUF9(8l95>,  IBUFl  0  ( 81 94  ) , 
& IBUF 11(81 33 ) <  IBUF 12 ( 8182  > ,  I BUF 13 ( 3181 ) ,  IBUFI4(B180) 

DIMENSION  IBUF (8193) 

DIMENSION  J7IM(2) 

INTEGER  TURN, SD,  SB,  SC, SD2 

EQUIVALENCE  <IBUF(14), IBUF14(1 }, IBUF13(2> ) 

EQUIVALENCE  ( IBUF( 12),  IBUF12C 1 ) ,  IBUFl 1(2),  IBUFl 0(3) , 

ScI  BUF9  ( 4  ) ,  IBUF8  <  5 ) ,  IBUF7(6),  IBUF6<7),  IBUF5(8),  IBUF4<9), 
fcIBUF3< 10),  IBUF2U0)  ) 

DT=1. /SAMP  RATE  OF  LASERS) 


DATA  DT/.  0005/ 


DT  IS  SET  FOR  2000  HZ 


C 

c 
c 
c 

cccccccccccccccccccccccccccccc; 

r 


cccccccccccccccccccccccccccccccccc 


c  GET  FILE  NAME 

C  OPEN  FILE  h 

C  WRITE  REC 

C 
C 

1  TYPE  8377 

SS77  FORMAT  ('  ENTER  6  CHARACTER  NAME:  ' , %) 

ACCEPT  8876,  NC,  (  NAMEF  (  I  ) ,  !--=!,  NC  ) 
FORMAT  (G.80A1) 

CALL  1RAD50(6, NAMEr  YNANS ) 


8S76 


■0  CD  -r~i  h  (J  (j 
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CALL  NAME I T ( XNAME  > 

CALL  OPEN 

N=  1 

IR  =  1 

TURN=1 

IWHEEL=1 

IWH=1 

TYPE  7010 

7010  FORMAT  ('  SHALL  1  TAKE  THE  WHEEL  PULSES  FROM  SWITCH?  ',*) 
ACCEPT  7011, IANS 

7011  FORMAT  <A1) 

IF( IANS.  NE.  IYES)  GO  TO  7111 
I  WH=5 
IWHEEL=16 
7111  TYPE  1005 

C 

C  START  DATA  ACQUISITION 

C 

C 

1005  FORMAT  <'  TURN  ON  START  SWITCH  TO  START  DATA  ACQUISITION ' ) 

C 

C 

I F  (  I  SW  (  2  ) .  NE.  1  )  GO  TO  6b 

IF<  <  1SW<  IWH).  EQ.  1  ).  AND.  ( 1SW(2).  NE.  0)  )  GO  TO  8 
IF ( ISW<2) .  EQ.  O)  GO  TO  1720 
0  IF(  (  ISW(  IWH).  EQ.  0).  AND.  <  ISW(2).  NE.  O)  )  GOTO  10 

IF  < ISW  <  2 ) .  EQ.  O )  GO  TO  1720 
GO  TO  9 

720  TYPE  1721 

721  FORMAT  <'  YOU  HIT  STOP  BEFORE  YOU  STARTED???'-'???') 

CALL  EXIT 


CALL  GTIM  <  JTIM2 ) 

CALL  CLASS <LA2,  LAI,  KA. LB2, LB  1 > KB. 

?yLC2,  LC1,  KC,  LD2,  LD1,  KD, 

MSD3,  ISD2,  ISD1,  ISDB3,  ISDB2,  ISDB1, 

&ISDC3,  ISDC2.  ISDC1,  ISDD3,  ISDD2,  ISDD1,  I WHEEL) 
CALL  GTIM(JTIM) 

SUMA=CONV ' LA2,  LAI ) 

SUMB=CONY  <  LB2, LEI ) 

SUMC=CONV LC2,  LC 1 ) 

SUMD=CONV ( LD2, LD 1 ) 

I AVGA=0 
I AVGB=0 
I AVGC=0 
1AVGD=0 

IF ( KA.  NE  O)  I AVGA=SUMA/FL3AT  <KA) 

IF(KB  NE.  0)  IAVGB  =  SUMB/FLOAT (KB ) 

IF(KC.  NE.  O)  1 AVGC=SUMC /FLOAT  <  KC ) 

IF(KD  NE.  0)  IAVGD=SUMD/FLOAT(KD) 
SSA=FIXSD( ISD1 ,  ISD2,  ISD3) 

SSB=FIXSDi ISDB1,  ISDE2.  ISDD3) 
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SSC=FIXSD( ISDCli  ISDC2,  1SDC3) 

SSD=F I  X SD  < I SDD 1 ,  ISDD2.  ISDD3) 

SA=0 

SB=0 

SC=0 

SD=0 

•4777  FORMAT  (  TESTQ=  ' .  F  1  3.  4  /  ) 

IF  <  KD.  NE.  0)  SD=SQRT < < SSB- < SUMD*SUMD > /FLOAT < KD ) ) /FLOAT ( KD ) ) 

IF  ( KB.  NE.  G)  SB=SGRT(  (  SSB-  (  SUMB*SUMB  )  /FLOAT  KB  )  )  /FLOAT  ( KB  )  ) 

IF  <  KC.  NE  0)  SC -SORT (  ( SSC- <  SUMC*SUMC ) /FLOAT ( KC ) )/FLOAT<KC) ) 

IF  <  KA.  NE  0)  SA=SQRT  (  SSA-  <  SUM A* SUM A  >  /  FLOAT  (  KA )  >  /  FLOAT  (  KA )  ) 
9995  IF<  ISU<4).  NE.  1  >  CO  TO  333 
IBUF  <  N )  =  I A VGA 
I BUF2  <  N )  =  I AVGB 
I BUF3 ( N ) = I AVGC 
I BUF4 ( N ) = I AVGD 
IBUF5  <  N ) =KA 
I BUF6  <  N ) =KB 
IBUF7CN)=KC 
I  BUFf3 ( N  >  =KD 
I BUF9 ( N ) =SA 
I BUF 1 0  ( N ) =SB 
IBUF1 1 <N)=SC 
IBUF 12 ( N ) =SD 
I D UF 13 <N>= JT I M2 ( 2 ) 

IBUF 1 4 ( N  5  =UT IM <  2 ) 

N-N+36 

333  IF  <  I SW  ( 3 )  .  EG.  1  )  TYPE  1112,  TURN,  IAVGA,  I  AVGB,  I AVGC,  I  AVGD,  KA, 

&  KB, KC. KD 

33  IF<  C  ISW(IWH).  EQ.  1  >.  AND.  <1SW<2>.  NE.  0>  )  GO  TO  33 
IF(  ISW<2) .  EG.  0)  GO  TO  1717 

34  IFC (ISW(IWH). EG.  0).  AND  < ISW(2) .  NE.  0) )  GO  TO  34 
IF < ISW <  2 > .  EG.  O)  GO  TO  1717 

IF  < ISW  <  4 ) .  EQ.  0 )  GO  TO  43 
1112  FORMAT  ('  ' ,  14,  X,  BI 6 ) 

IF(N.  GE.  3192)  CALL  WR ITE ( IBUF,  JR > 

43  IF(  (ISW<  IWH).  EQ.  1  ).  AND  ( ISWC2).  NE.  0)  )  GO  TO  43 
I F  < I SW  <  2 )  EG.  0)  GO  TO  1717 

44  IF  <  <  ISW  (  IWH) .  EQ.  0 )  AND.  <  I  SW  \  2  <  NE.  0>)  GO  TO  44 
IF (ISW (2). EG  0)  GO  TO  1717 

C 

TURN=TURN+1 

IF  (  N.  LT.  81  92)  GO  TO  9 

N—  1 

IR=IR+32 
GO  TO  9 
C 
C 
C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

c 

C  NOT  START 


non  non  ►*  ~  n  o  r> 
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TYPE  EXIT 


717  TYPE  1007,  TURN 

007  FORMAT  !'  END  OF  TEST  AFTER  '.16,  '  FEET'//) 
IF! ISW!4).  EG.  0)  GO  TO  1 
DO  1562  1  =  1,  12 
562  IBUF'.N+I  — 1  )=-l 

CALL  WRITE! IEUF, IR) 

CALL  WAIT 
CALL  CLOSEX 
GO  TO  1 
END 


FUNCTION  FIXSD! ISD1,  ISD2,  ISD3) 

X  X  =  I SD 1 
YY=ISD2 
ZZ=ISD3 

IF  (  XX  LT.  0  )  XX=6553o.  +XX 
IF  !  YY.  LT.  O.  )  YY=65536.  t-YY 
I F <  Z Z .  LT.  0.  )  ZZ=65536.  +ZZ 
FIXSD=XX+65536  * < YY+65536.  *ZZ) 

IF  <  F I  XED  LT  0.)  TYPE  12  ,  ISD1,  ISD2,  ISD3,  FIXSD 
12  FORMAT (  2X ,  '*■*##«■**■«■*■**■»■>?■**#■*#■*■#■*■*#  /,  2X, 

*  '  ISD1  -3,  FIXSD  *  315,  X,  G12.  6> 

RETURN 

END 


FUNCTION  C0NVU2,  ID 
Y=I  1 

IF  (  Y.  LT.  0.  )  Y=65536.  +Y 
C0NV=65536.  *FLOAT < 1 2 >+Y 
IF(CONV.  LT.  0.  )  TYPE  12,  12,  Il.CONV 

12  FORMAT!/,  X,  &*&*&*&•*&■*&*&*', 

5*  /  2X,  '12,  12,  CONV  =  ',  217,  X,  G16.  3) 

RETURN 
END 
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Program  CLASSY#  a  subroutine  to  GATHER. 

»  *  »»  i  >  k  t  »  i  i  >  *  t  »  \  *  t  \  't  't  t  ‘.  i  t  t  i  i'l'ti'tii  >>  i  '  \  i  t  t  )  i  t 


CALL  CLASS < LAE,  LAI,  KA,  LB2,  LB1.  KB#  LC2.  LC1,  KC, 

It  LD2#  LD1#  DK,  LSD  A3.  LSDA2  •  LSDA1,  LSDB3.  LSDB2.  LSDB1. 
&  LSDC3# LSDC2. LSDC1,  LSDC3,  LSDD2# LSDD1,  IWH) 


I  =  I SW ( LSW ITCH) 


t  P  >**>>*>>*  l  >  i  *  »  t  »»  t  i  >  t  >  *  >  i  >>!>>  >  *  >  >  i  \  't  \  t  \  t 


.  ENABL  AMA 
.  GLOBL  CLASS,  ISW 


;  DEFINE  SYMBOLS  AND  REGISTERS 

i 

AR=  1 77550 

BR=3 77540 

CR=1 77230 

DR-177530 

ER=177510 

2=7.0 

Z  2=7.1 

MASK=7.2 

R=-.'5 

PC -7.7 

START=2 

i  START  LASERS  AT  ENTRY  POINT  “'CLASSY" 

CLASS:  INC  AR 

INC  BR 

INC  CR 

INC  DR 

MOV  #2,  Er? 

;  CLEAR  SAMPLE  COUNT 

i 

CLR  KA 

CLR  KB 

CLR  KC 

CLR  K  D 

t 

,  CLEAR  AVG  REGISTERS 

i 

CLR  SUM1A 

CLR  SUM2A 

CLR  SUM1B 

CLR  5UM2B 

CLR  S'.  W  3  C 
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CLR  SUM2C 

CLR  SUM ID 

CLR  SUM2D 

CLEAR  RMS  REGISTERX 

SSI  A 
SS2A 
SS3A 
SS1B 
SS2B 
SS35 
SS1C 
SS2C 
SS3C 
SS1D 
SS2D 
SS3D 

SET  WHEEL  SWITCH  BIT 
MOV  <*62<R),WH 

i  CLEAR  THE  HALF  CYCLE  DETECTOR  &  SET  MASK 

f 

CLR  HALF 

MOV  #100007,  MASK 

* 

;  LASER  "A"  ACQUISITION.  AVG.  Z-  RMS  CALCULATIONS 

i 

A:  TSTB  AR  ; IS  LASER  "A"  READY 

BPL  B  i  NO  -  SO  GO  GET  LASER  ,,B*' 

CMP  KA,  #32767. 

BEG  B 

MOV  AR+2, 2  ; YES-GET  DATA  FORM  PARALLEL  CARD 

INC  AR  ; RESTART  LASER 

BIC  MASK, Z  i CLEAR  UNWANTED  BITS 

ROR  Z  i  (THESS  CAN  BE  DELETED- 

ROR  Z  ;  IF  CABLE  IS  REWIRED  SO  THE 

FOR  Z  ;  DATA  IS  RIGHT  JUSTIFIED) 

ADD  Z , SUM 1 A  > ADD  DATA  TO  SUM 

ADC  SUM2A 

MOV  Z.TMP 

MUL  TMP, Z 

ADD  Z2, SSI  A  ; ADD  SQUARE  TO  RMS  SUM 

ADC  SS2A 

ADD  Z.SS2A 

ADC  SS3A 

INC  KA  ; INC  THE  SAMPLE  COUNT 

LASER  "B" 


CLR 

CLR 

CLR 

CLR 

CLR 

CLR 

CLR 

CLR 

CLR 

CLR 

CLR 

CLR 


TSTB 
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BPL 

C 

CMP 

KB,  £32767. 

BEQ 

C 

MOV 

BR  +  2,  Z 

INC 

BR 

BIC 

MASK, Z 

ROR 

Z 

ROR 

Z 

ROR 

z 

ADD 

Z,  SUM IB 

ADC 

SUM2B 

MOV 

Z,  TMP 

MUL 

Trip,  z 

ADD 

12,  SS1B 

ADC 

SS2B 

ADD 

Z,  SS2B 

ADC 

SS3B 

INC 

KB 

LASER  "C" 

TSTB 

CR 

PPL 

D 

CMP 

KC- £32767. 

BEQ 

D 

MOV 

CR+2.  Z 

INC 

CR 

BIC 

MASK.  Z 

ROR 

Z 

ROR 

Z 

ROR 

Z 

ADD 

Z,  SUM  1 C 

ADC 

SUM2C 

MOV 

Z,  TMP 

MUL 

TMP,  Z 

ADD 

12,  SS1C 

ADC 

SS2C 

ADD 

Z, SS2C 

ADC 

SS3C 

INC 

KC 

LASER  "D" 

TSTB 

DR 

BPL 

E 

CMP 

KD, #32767 

EEQ 

n 

MOV 

DR+2.  Z 

INC 

DR 

BIC 

MASK.  Z 

ROR 

Z 

ROP 

T 

A. 

ROP 


7 
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ADD  Z,  SUM ID 

ADC  SUM2D 

MOV  Z,  TMP 

MUL  TMP,  Z 

ADD  Z2, SS1D 

ADC  SS2D 

ADD  Z.  SS2D 

ADC  SS3D 

INC  KD 


READ  AND  PROCESS  THE  SWITCHES 


E:  MOV 

BIT 
BEQ 
TST 
BNE 
BIT 
BNE 
MOV 

GOTQA:  JMP 

SNDHLF:  BIT 
BEQ 


ER+2,  Z 
#START,  Z 
EXIT 
HALF 
SNDHLF 
WH,  Z 
GOTOA 
#1,  HALF 
A 

WH,  Z 
GOTOA 


;  IS  START  SWITCH  STILL  SET? 

;NO  -  GO  EXIT 

> YES-ARE  WE  IN  FIRST  HALF  OF  CYCLE 
;  NO-GO  TO  SECOND  HALF 

;  YES  -  ARE  WE  NOW  IN  FIRST  HALF  OF  CYCLE? 
;  NO-GO  ACQUIRE  MORE  DATA 
;  YES-SET  SECOND  HALF  FLAG 
;  AND  GO  ACQUIRE  MORE  DATA. 

; SECOND  HALF  -  ARE  WE  NOW  IN  FIRST  HALF? 

;  NO-GO  ACQUIRE  MORE  DATA 


EXIT  WITH  ARGUMENTS 


EXIT:  TST 

MOV 
MOV 
MOV 
MOV 
MOV 
MOV 
MOV 
MOV 
MOV 
MOV 
MOV 
MOV 
MOV 
MOV 
MOV 
MOV 
MOV 
MOV 
MOV 
MOV 
MOV 
MOV 
MOV 
MOV 

* 

RTS 


(R  )  + 

SUM2A, ©<R)+ 
SUM1A, e<R)+ 

KA,  e<R )+ 
SUM2B, @<R>-<- 
SUM1B,  e(R  )-f 

KB,  <5  <  R  )-*- 
SUM2C, £<R>+ 
SUM1C, @(R)+ 

KC ,  ®  <  R  )  + 
SUM2D,  ©<R)  + 
SUM1D,  @  ( R  >  •+• 

KD,  @<R>  + 
SS3A, £<R  >  + 
SS2A,  *2HR  ; 
SSI  A,  0(R  *■*- 
SS3B, @(R )+ 
SS2B, ®<R  >  + 

SSIB,  e<R)+ 
SS3C, @<R  )  + 
SS2C,  £<R  ')* 

SSIC,  ©<R )  + 
SS3D, <*<R  '  + 
SS2D, G<Ri+ 

SSID,  @  ( R  )  + 

PC 
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DATA 


WH 

.  WORD 

0 

KA 

.  WORD 

0 

KB: 

.  WORD 

0 

KC: 

.  WORD 

0 

KC 

.  WORD 

0 

SUM  1  A 

WORD 

0 

SUM2A 

WORD 

0 

SUM  IB 

.  WORD 

0 

SUM2B 

.  WORD 

0 

SUM  1C 

.  WORD 

0 

SUM2C 

.  WORD 

0 

S'JM  ID 

.  WORD 

0 

SUM  2D 

.  WORD 

0 

SSI  A 

.  WORD 

0 

SS2A: 

.  WORD 

0 

SS3A : 

.  WORD 

0 

SS 1 B : 

.  WORD 

0 

SS2D . 

.  WORD 

0 

SS3D 

.  WORD 

0 

S?1  c 

WORD 

o 

V.* 

SE2C  • 

.  WORD 

w 

cnr  • 

.  WORD 

0 

SS1D: 

.  WORD 

0 

SS2D: 

.  WORD 

0 

mn  • 

wuw ’t-f  . 

.  WORD 

0 

HA_F 

.  WORD 

0 

TMF  : 

.  WORD 

0 

FUNCTION  SUBROUTINE  TO  READ 
ANY  BIT  OF  SWITCH  CARD 


MOV 

#2.  ER 

; START  SAMPLING 

MOV 

#1  •  Z 

; ROTATE  ONE  BIT  TO  TEST  LOCATION 

MOV 

£2  <  R  ) ,  R 

Ci_R 

•1  CL 

DEC 

k 

ASHC 

R.  Z 

EIT 

I,  EFN2 

i TEST  IF  THIS  BIT  IS  SET  IN  CARD 

BcG 

NOTSET 

; NO  -  GO  EXIT  WITH  ZERO 

INC 

12 

; YES  -  GO  EXIT  WITH  ONE 

RTS 

PC 

;E>::T  -  SUBROUTINE  "ISW" 

END 

ENABL 

AM  A 

.  GLGBL 

ISW 

A 


N0T3FT 
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ER=177510 
2=7.1 
r=>:5 
22=7.0 
PC  =7.7 


Program  DISK.  MAC/  a  subroutine  to  GATHER. 
;  DISK  READ  WRITE  SUBROUTINES 
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.  MCALL  .  LOOKUP,  .  FETCH 

.  MCALL  . READW.  .  WAIT,  .  WRITE,  .  LOCK,  .  ENTER,  .  CLOSE 
GLOBL  OP END 
.  GLOBL  NAME IT 

.  GLOBL  OPEN, CLDSEX,  READ, WRITE, WAIT,  IER 
ENABL  ANA 

R0=7.0 
R  1=7.1 
R2=7.2 
R  3=7.3 
R  4=7.4 
R  5=7.5 
FC-7.7 


OPEN. 

.  LOCK 

. FETCH 

aCORSPC, #NAME 

CLR 

FIRST 

MOV 

#AREA,  7.5 

MOV 

#1, 7.4 

.  ENTER 

R5,  R4,  #NAME,  #■ 

BCS 

ERR 

MOV 

#1,  E 

RTS 

PC 

ERR 

CLR 

RTS  PC 

E 

NAME I T : 

MOV 

2<R5),  R 1 

MOV 

( R 1 >+, NAME+2 

MOV 

( R 1 > , NAME+4 

i 

RTS 

PC 

i 

NAME: 

.  RAD50 

/DK  / 

.  RAD50 

/DAT/ 

.  RAD 50 

/A  / 

.  RAD50 

/DAT/ 

AREA: 

» 

.  BLKW 

10 

i 

IER: 

MOV 

E,  RO 

/ 

RTS 

PC 

/ 

CLOSEX 

MOV  #1,R4 

CLOSE 

R4 

RTS 

PC 

CLOSEX 


#1,  R4 
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OPENO: 

LOCK 

FETCH 

CLR 

MOV 

MOV 

. LOOKUP 
BCS 

MOV 

RTS 

ERG: 

CLR 

RTS 

E . 

.  WORD 

READ: 

MOV 

MOV 

MOV 

MOV 

. READW 
BCS 

MOV 

RTS 

RDERR: 

CLR 

RTS 

FIRST: 

.  WORD 

WRITE. 

MOV 

MOV 

MOV 

MOV 

FIR: 

WRITE 

BCS 

MOV 

MOV 

RTS 

WTERR : 

CLP. 

FTP 

WAIT. 

MOV 

WAIT 

ECS 

NOV 

RTS 

W4ER  CLR 
RTS 

CORSPC  .  BLKW 
.  END 


#CORSPC. #NAME 
FIRST 
*AREA,  R5 
#1 ,  R4 

R5i  R4<  #NAME 
ERO 
#1,  E 
PC 
E 
PC 
1 


2 <  R5 ) » R4 

@4 ( R5 ) • R 1 

#1,  R3 

#AREAj  R2 

R2,  R3»  R4,  #8192. 

RDERR 

#1 ,  E 

PC 

E 

PC 

0 


2(R5)i  R4 

£4 ( R5 ) »  R 1 

#1 ,  R3 

#AREA,  R2 

R2.  R3<  R4,  #8192. 

WTERR 

#1 . FIRST 

#1,  E 

PC 

E 

FC 


#1-  R3 
R3 

WAER 
#1  .  E 
PC 
c 

PC 

400 


R  1 


R  1 


ru  n 
N  rj 
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Program  LIST,  used  to  preview  the  data. 

C  THIS  IS  LIST  FDR  IT  TAKES 

C  A  DATA  RILE  CREATED  BY  GATHER  FOR  AND  RE-WRITES 
C  IT  AS  AN  ASCII  CODE  FILE.  EITHER  ON  DISK, 

C  PRINTER  OR  SCREEN.  SEE  PROGRAM  FOR  THE 
C  ORDER  THE  DATA  IS  PRESENTED 
C  N=NUMBER  READINGS 

C  S=STD  DEV  OF  ONE  OF  THE  GAGES  OVER  FOUR  INCHES 

C 

DIMENSION  I  BUFF (81 ^3) 

BYTE  NAMEF  <  SO )  ,  NAMES ( BO  > 

BYTE  I  TT.  I  AA,  ILL,  IDD 

DATA  ITT,  It  L.  IDD/1HT,  1HL,  1HD/ 

I  BOZ0=0 
TYPE  5555 

5555  FORMAT  <  / ,  X ,  '  LISTING  ON  TT,  LP  OR  DISK''-  (T.L  OR  D)  ',*) 
ACCEPT  55,  I  AA 
55  FORMAT  (AD 

I F  <  I  AA.  NE.  IDD)  GO  TO  22 
TYPE  20 

20  FORMAT  (  /,  2X>  '  OUTPUT  FILENAME  ?  ' ,  S> ) 

ACCEPT  21  NOC,  (NAMED  .>),  JN,  NOC  > 

FORMA  I (G, 80A1 ) 

NAME2  <  NGC+ 1 )=0 

CALL  ASEIGNdl,  NAME2,  NOC,  'NEW') 

CONTINUE 

FORMAT  (/,2X,  '  DATA  FILE  '.20A1,/) 

1  =  17 

DO  10  11=1,32767 
IF( II  GT. 1 )  BOZO= 1 

CALL  READC  ( A,  B,  C,  D,  NA.  NB,  NO,  ND,  SD1,  SD2,  SD3,  SD4.  IS,  IE,  I  ) 
DT=FLOAt < IE-IS) 

I DT  = I E- 1 S 
ISD1  =  IFI X  <  £D1 ) 

I SD2= I F I  X ( SD2 , 

I SD3= I F I  X ( SD3 ) 

I  SD4=  I  FI  >.'  (  SD4  ) 

I A= I F I  X . A ) 

I B- J  F I  X 1 D ) 

I  C  =  I F I  X  '.  (  i 
I D= I F I  X  t  D , 

IF (A  LT  G)  GO  TO  1 1 

1F(IAA,  EG  IDD)  WRITE  (11.  cjO  )  I  A,  IB,  IC  ID,  NA,  NB,  NC,  ND.  ISD1, 

I SD2,  ISD2,  ISD4,  IDT 

I  F  (  I  AA  EG  ITT)  TYPE  1000.  I  T,  IA,  IB,  IC,  ID,  NA,  NB,  NC,  ND.  ISD1, 
ISD2,  ISD3-  ISD4,  IDT 

IF(  IAA.  EG.  ILL  )  PRINT  1  000,  1  1 ,  1  A,  I B ,  I C ,  I D.  NA,  NB ,  NC,  ND.  I SD1 , 
ISD2.  ISD3.  I SD4,  IDT 
10  CONTINUE 

II  IF  (  IAA  EG  IDD)  CALL  CLOSE(U) 

1000  FORMAT  (  X,  1315,  X,  14) 

fjj  FQRMAT(  X,  1315) 
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GO  TO  1 
END 

C  SUBROUTINE  READC  <  A,  B.  C.  1).  NA,  NB.  NC,  NC.  SD1 .  SD2.  SD3.  SD4,  IS. 

C  IE,  I  ) 

C 

C  A,  B,C,D  =  LASERS  A,  B>  C»  D 

C  NA, NB. NC, ND  =  NO  OF  SAMPLES 

C  SD1 , SD2, SD3, SD4  =  STANDARD  DEVIATIONS 

C  IS  =  TIME  TICKS  READ  BEFORE  CLASS 

C  IE  =  TIME  TICKS  READ  AFTER  CLASS 

SUBROUTINE  READC  <  A,  B,  C.  D,  Nr,,  NB,  NC,  ND,  SDl ,  SD2,  SD3,  SD4, 

I  START  -  IEMD,  I  ) 

COMMON  /RMECOM/ IBUF ( 16, 512) 

DATA  I OPEN/O/ 

DATA  LR/-1/ 

BYTE  NAME (20) 
r 

I F ( I  OPEN  EQ.  1 )  GO  TO  1 
I OPEN= 1 
WRITE'S, 2) 

2  FORMAT  (  / ,  3X ,  'RAW  DATA  FILE  NAME  ...  ',*> 

ACCEPT  5, NOC ,  < NAME ( I Z Z ) ,  IZ7-1.N0C) 

NAME ( NOC+ 1 ) =0 
5  FORMAT <Q, 2CA1 ) 

CALL  ASSIGN ( 1, NAME, NOC,  'OLD') 

DEFINE  FILE  1  < 1000, 0192, U.  IV) 

WRITE (6. 4)  (NAME (ITT), 1TT=1, NOC) 

4  FORMAT ( / , 3X , 'DATA  FILE  NAME  ',80A1) 

WR I TE ( 5, 3) 

3  FORMAT  (  /  ,  3X,  'ROARRRRR  •  !  '  CHOMP  CHOMP  CHOMP',*) 

1  =  17 

1  1R=< < 1-1 )/512>+l 

IF  < IR.  EQ  LR)  GO  TO  30 

READ ( 1  ' IR, END= 1717)  <  ( IBUF'MZ, LZ), MZ=1,  16). LZ  =  1,  512) 

LR  =  I R 

30  INDX=MOD( I,  512) 

JF<  JNDX.  EG.  0)  I NDX  =  5 1  2 
A= I BUF ( 1,  I ND  X ) 

D=IRUF(2, I NDX ) 

CV;  BUF  (  3.  I  NDX  ) 

D= IBUF <c,  I NDX  ? 

NA= IBUF ( 5,  INDX) 

NB= I BUF ( 6, INDX ) 

NC=IBUF ( 7 , INDX) 

ND= I BUF (0,  INDX > 

SDi=IBUF (9, INDX ) 

BD2= IBUF ( 10. INDX ) 

SD3= I BUF ( 11, INDX) 

£D4= I BUF ( 12,  JNDX ) 

I START= I BUF (13-  INDX) 

I END= IBUF (  14,  INDX ) 

1  =  1  +  1 
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IF(A.  NE. -1)  RETURN 

CONTINUE 

RETURN 

END 
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Program  CALC,  used  to  reduce  the  data. 


c  this  is  CALC.  FOR.  It  uses  another  way.  The  terror  of  the  bwy. 

C  data  ab,  be.  cd/1 1.  974.  12.  149.  83.  75/ 

data  ab.  be.  cd/12.  04.  12.  .  34.  / 
data  ascl.  dscl.  dflt/1.  0.  1.  0.  5h*****/ 
byte  name(30).  cue2.  YES.  NO.  preue.  PP.  PD 
data  YES,  ND,  key/lhY,  IhN,  0/ 

dimension  bb( 20 ),  c  c (20 ) .  dd( 20 ) ,  b ( 20 ) ,  c ( 20 ) . f c  tra  <  20 ) . 

&  d (20 ) .  aa  s cron (20 ) ,  aascmx (20) ,  dm in (20) ,  d max (20) 

S'.  .  gr ad  ( 20 ) .  i num  ( 20 ) ,  aa  s  tr  t  ( 20) ,  e  1  evl  ( 20 ) ,  el  ev2  (20) , 

S'  aascal  (20) »  pr  ( 20) .  oef  d  i  f  ( 20) » 

Sc  saend  (20) ,  pltct  (20) 

integer  fa,  fb. fc. fd,  to, consec. deflct,  reset, cue3 

common  /p 1 1 / i p 1 1 ,  pmi n,  dmi n 

common  a (20) . aa ( 20) ,  def 1 ( 20) 

c  ue2=N0 

c  ue3=0 

uir  i  te  (  5,  92) 

92  f ormat ( ///, 3x. '10  -  4  good  buddy  !!  Ready  to  run  CALC  ?', 

&  3x,$) 

accept  73,  :  d  z 
write  (  5,  93) 

93  f ormat (/, 3x .  'Will  you  want  an  extended  pr intout?  ' , 2x, * ) 
accept  73,  preue 

do  122  i z i=l . 20 
p 1 tc t ( i z i  ) =0.  0 
aas  tr  t  <  i  z  i  )  =0. 
aaend  (  i  z  i  )  =0. 

122  continue 

if  (cue2.  eq.  NO)  go  to  85 
78  urr  x  te  (  5,  68  ) 

JJ  =  17 

68  format(/,3x,  'Type  in  full  name  of  profile  plot  file...  ' 

S-  .  $) 

call  ass i gn ( 8, . -1 . 'new', 'CC') 
wr  i  te  (  5.  69 ) 

69  f  ormat  </,  3x .  'Type  in  full  name  of  def  1 .  plot  file... 

$  S  5 

call  assign  (9,, -1,  'new'.  'CC') 

85  continue 

xn=  <ab+bc)/bc 

if  <  c  ue2.  eq.  NO )  write(6,43) 

43  format  ( 1  h  1 ,  ///,  45x,  '*#■***  PURDUE  -  WATERWAYS  LASER  SYSTEM'. 

&  '  ■•HHMUt-  '  ,  //  ) 

c  i n i t i a  1 i z a t i on s 
ictt=0 
incre=0 

if  <cue2.  eq.  NO)  afctra=0  0 
if  (cue2  eq.  NO)  afctrd=0.  O 
sa=0. 
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sc=0.  0 
sd=0.  0 
d  sb  =0. 
dsc=0. 
d  sd=0. 
nf  a=0 
nf  b=0 
nf  c=0 
nf  d=0 
nnfa=0 
nnf b=0 
nnf c=0 
nnf d=0 
inn=00 
inni=0 
JJ  =  17 
aval  i  d=0. 
dva 1 i  d=0 

if  (cue2  eq.  YES)  go  to  63 

1  format(/»2x,  '  Total  number  of  invalid  readings  is  *  '. 

*  5<  i4,  x  ),  /) 

S  forraat(/,  2*.  '  No  of  data  pts.  is  =  ',  i7) 

22  formate/, 3x,  'start/pnd  -  ',£) 

21  formate  f  20.  5  ) 

20  formate/, 2x,  '  start/end  =  ' , 2( i 10,  x >, / > 

J  J  ~  1 7 

call  redcal  (i,  z ,  z ,  z,  z,  z,  z,  z,  z,  z,  z,  z,n,n,n,n.  z,  z,  z,  z,  jj) 
scheme  for  getting  out-of-line  constants  (delta  a, d ) 

JJ  =  17 

uir  i  t-  e  (6,  94  ) 

94  formate/,  3x,  'Bad  A,  B,  and  C  gage  readings  as  follows  —  '. 

Sc  /,5x,  '  1  o  c  at  i  on  ' .  6x  ,  'gages',/, 

&  8x»  'ft.  ',8x>  'A  B  C  ' .  /  ) 

do  23  innn=l, 32575 

call  redcel(ax,  b  x .  cx,  dx,  ascl,  dscl.afctra,  afetrd, 

&  az,  bz.  cz,  dz,  na,  nb,  nc,  nd,  stda,  stdb,  stdc,  stdd,  jj) 
if ( a z .  1 1 . 0.  )  go  to  34 

call  check  <  a  z ,  bz.  cz.  dz,  na,  nb,  nc,  nd,  fa,  fb,  fc,  fd) 
nnfa=nnf6+fa 
nnf b=nnf o  +  f b 
nnf c=nn? c+f c 
nnf d=nnf d+f d 
inn=inn+l 

ifffaeqlor.  fb  eo  1  orfeeqi)  write(6,95)  innn,  fa.  f  b.  f  c 

95  f ormat ( 7 x , i 3. t 1 9 , 2i2 > 

if(fa  eq.  1.  or.  fb.  eq  1  or  fc.  eq.  1)  go  to  90 

sa-sa+ax 

sb  =  sb+b  x 

SC*5C+C X 

avalid=avalid+l 

of  c  t’a=  (  sa  'a  /a  1  i  d  )  -  '■  (  s<  /aval  m))-(  (bc*ib  )/bc  '*(  (sc-*b)/a  valid)) 
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90  if(fb.  eq  1.  or.  fc.  eq.  1.  or  fc.  eq.  1)  go  to  91 
d sb=dsb-*-b  x 

dsc=dsc+c  x 
dsd  =  dsd-t-dx 
dval i d=d va 1 i  d  +  1 . 

ofctrd=(dsd/d valid  )-<  <dsc/dvalid)+<cd/bc>*(  (dsc-dsb)/d valid)  ) 

91  continue 
23  continue 

34  incre=nnfsTnnfb+nnfc+nnfd 
J  J  =  17 

wr  i  te  (  5,  1)  nnfa,  nnf  b ,  rsnf  c ,  nnf  d ,  incre 

wr  i  te  (6,  1  )  nnfa,  nnfb.  nnfc.  nnf  d,  incre 

wr  i  te  <  6,  8 )  inn 

wr  i  te  <  5.  8 )  inn 

wr i te  <  5, 22) 

accept  21,  rist 

wr  i  t  e  <  5,  22 ) 

accept  21 ,  r i end 

l s t-  i f i x l r i st ) 

iend=ifix (riend ) 

wr  i  te ( fc,  20 )  ist.  iend 

write(5,20)  ist,  iend 

u>r  i  t  e  ( 6,  ^'6  ) 

96  formate/, 3x. 'Bad  readings  between  start  and  end  are  -  ',/• 

$<  5  x  ,  'location',  6x,  'gages',/, 

&  8x,  'ft  ',  8x,  'A  B  C  '.  /) 

i  f  (  f  a  eq  1.  or  fb.  eq  1  or.  fc.  eq.  1)  wr i te ( 6.  95)izz,fa.  fb.fc 
do  11  i z z  =  ist, iend 

call  read c ( a  z ,  fcz,  cz,  dz,  na.  nb, nc, nd,  stda,  stdb,  stdc, stdd, 

&  iistar, iiend,  jj) 

call  check (az, bz, cz, dz, na, nb, nc, nd, fa, fb, fc, fd> 
if  (fa.  eq  1.  or  fb  eq  l  or  fc.eq  1)  write<6,95)izz»fa,  fb,fc 
if  (az.  It.  0.  )  go  to  9 
nfa=nfa+fa 
nf b=nf b+f b 
nfc=nfc+fc 
nf d=nf d+f d 
11  continue 

9  ictt=nfa+nf  b+nf  c-enf  d 

J  J  =  l" 

25  forma  t  (  / »  3x ,  'Gage  spacing  ab,bc  am  cd  =  ',3<fl0.  5,  x),/) 

50  f  orma  t  <  2x ,  i  5,  x ,  8  (  f  10  5,  x  >  ) 

uirite(5,39)  ofetra,  ofctTd,nfa,nfb,nfc» 

&  n  f  d ,  i  c  1 1 

uir i t e (6,  39)  ofetra,  ofctrd,  nfa,  nfb,  nfc, 

Sr  n  f  d ,  i  c  1 1 

39  format ! /, 2*,  '  Calculated  delta  a,d  =  ' , 2 ( f 1 0.  6.  x  ) , 

S-,  //,3x,  'There  were  ',x,5(i6,  x).'  inva  lid  readings',/, 

&  4x,  'in  between  the  start  and  end.  './) 

wr  i  t  e  <  5,  42 ) 
accept  41.  afetra 
wr i t  e ( 5.  401 
accept  41,  afetrd 
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42  format(/,2x.  '  input  fctra  please.  .  .  '.*) 

41  f  ormat  (  f  10.  6) 

40  format(/.2x»'  input  f  c  tr  d  please..  ' .  * ) 

write(6.3B)  afctra,  afctrd 
38  format(3x<  'Input  delta  a.  d  —  2(f  10.  6.  x  )  ) 

write(6.25)  ab.bc.cd 

c  establish  the  arbitrary  distances  from  gages 

63  cc(l)=l 

b  b ( 2 )  =  1 . 
d ef I c  t=0 
s  umd  e  f =0.  0 
ssd e f  =0.  0 
i p 1 t— 0 
1  =  1 

c  AHEM!  now  do  the  first  seven  points,  for  which  the 
c  deflections  at  gage  d  cannot  be  calculated. 

JJ  =  17 

i  f  <  c  ue2.  eq.  YES)  ibbi=ist 
if  (cue2.  eq.  NO)  write<6.  54) 

54  format!  ///,  13(  '=»=========  '  >  ) 

if  <cue2.  eq  MO)  write (6, 53) 

53  format  (//.  t3,  'FT.  tl2,  'AA'.  t20,  'BB  '.  t29,  'CC  ', 

l:  t38>  'DD',t47.  'DEFL  ' .  t56.  'THETA',  t65,  'STD  A', 

Sc  t74,  'STD  B  ',  t83,  'STD  C  t92,  'STD  D',  tl02,  'FLAGS'.  //) 

call  advnce  <  ri st ) 

31  iplt=iplt+l 

i  f  (  c  ue2  eq  YES  and.  cue3.  cq.  0)  wr i te ( 5. 81 > 

81  format!//. 3x.  'Input  starting  el evati on  in  feet.  ..  '*> 
if ( cue2  eq.  YES.  and  c ue3  eq  0)  accept  82.  elevl(iplt) 

82  format(flO.  5) 

i f ( c  ue2  eq.  YES.  and.cue3  eq  0)  wr  i  te  (  5.  83) 

83  format  <//.3x.  'Input  ending  elevation  in  feet. 

if  (cue2.  eq  YES.  and.  cue3  eq  0)  accept  82.  elev2(iplt) 

if(cue2.  eq.  YES.  and.  cue3.  eq.  0)  el  evl  ( l  p  1 1 )  =e  1  evl  (  i  p  1 1 )  *12 

if(cue2.  eq.  YES.  and.  cue3.  eq.  0)  e 1 ev2 ( i p 1 t ) =e 1 ev2 ( i p 1 t >*1 2 
if  <  c  ue3.  eq.  1  )  call  ex  i  s.  <  aascmx  (  ip  1 1 ) .  aascmn  (  i  p  1 1 ) . 

S<  iplt.  pr(  iplt)  .  8.  'PROFILE  -  ENLARGED 

S-:  elevl(iplt).  e 1 e  v2 ( i p 1 1 ) .  1  ) 

if(cue3.  eq.  1)  call  asisldmax ( ip  It ). dmin (iplt).  iplt. 
defdif  (-pit).  9,  'DEFLECTIONS  -  ENLARGED 
S<  e levl  (  i  p  1 1 ) ,  el ev2 (  i  p  1 1 ) .  2) 
if(cue3.  eq.  1)  go  to  71 
dma  x ( i p 1 1  >  =-l e36 
dmin( i p 1 t ) = 1 e36 
a  a  s  c  m  x ( i p 1 1 ) =- 1 e  36 
aascmn! iplt)=le36 
71  do  3  i b b i  =  i st.  32573 

if(ibbi.gt. iend)  go  to  44 

if(cue2.  eq  NO)  p 1 tc t < i p 1 1 ) =p 1 1 c t ( i p 1 1 >  +  l . 
iri=(ab+tic  >/bc 

call  ~ed  ca  2 ( a  < i  +2  ) ,  1  i  ■  c(i),  d  x  x  ,  asc  1 .  dscl.  a  fctra. 
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&  afctrd.az.bz,  cz<  dz,  na.nb.nc.nd.  stda,  stdb.stdc.  stdd.  . 

if  <az.  It.  0.  )  go  to  44 

call  check  <  a  z ,  bz> cz.  dz> na»  nb ,  nc. nd»  fa.  f  b ,  f  c »  fd) 
if(fa.  eq  0.  and.  fb  eq.  0.  and.  fc  eq.  0)  go  to  32 
i  f  <  c  ue2.  eq.  NO)  write(6,26)  ibb  i.  fa.  f  b»  f  c.  f  d 
26  format  <  2x ,  i4,  x ,  'reset'.  1 100,  4i3) 

if(cue2.  eq.  NO  >  uir  i  t  e  (  5 ,  lOllibbi,  fa.  fb,  fc 
101  format!/.  3x,  'Bad  data  point,  at  location  '.i4,3x.4i2) 
i st= i bb i +1 
bb (2)=1 . 
cc  ( 1  >=1. 

x  f  (  i  st.  gt.  iend.  and.  i.  gt  1)  aaend(iplt)=aa(i+2> 
i  =  l 

if  (ist.  gt.  iend)  go  to  46 
if(cue2.  eq.  YES.  and.  cue3.  eq.  1)  key  =  l 
if(fcey.  eq  1)  go  to  100 
go  to  31 
32  continue 

c  the  order  the  algorithm  requires  the  data  is  a  b  c  d. 

c  a,  b  and  c  are  the  prof ilometer.  and  d  is  the  deflection  pt. 
c  First,  note  that  the  gage  spacing  will  affect 
c  the  algorithm.  Gage  spacing  factor  is  xn. 

c  Now  to  calculate  the  profile  point. 

theta=atan( <c(i )+cci i >-b ( 1+1 )-bb( i  +  1  ) )/bc ) 
c  theta=0. 

aa(i+2)=c(i  )*cos(theta)+-tc(i  )  -a  <  i  +2 )*cos(theta)-xn* 
(c(i)*c05(theta)+cc!  i)— b(i+l)*cos(theta)-bb(i+l)) 
i  f  <  c  ui?2.  eq.  NO.  and.  i.  eq.  1)  aastrt(iplt)  =aa  <  i+2> 
i f ( c ue2.  eq  YES )  call  sea  1 e ( e 1 evl ( i p 1 1 ) .  e  1  ev2< i p 1 1 ) , 

&  aastrt<iplt),aaend<iplt). pltct(iplt)»aa!i+2>, 

&  aa seal  l  i+2 ) ,  i  ) 

if  ;cue2.  eq  YES.  and  cue3.  eq.  O  >  aascmx  ( ipl  t  )=amax  l(aascmx(iplt) 
S<  aa  s  ca  1  <  i  +2 )  ) 

if(cue2.  eq.  YES.  and .  cue3.  eq.  O)  aascmn(iplt)  =aminl  ( aasemn  (  i  p  1 1 ) 
$<  aascal  (  i+2)  ) 

64  if (cue3.  eq.  1)  call  ppt<aascal(i  +2 ),  ibbi,aascmn(iplt),pr(iplt) 
&  8,  1.  > 

17  continue 

if  (preue.  eq  YES.  and. 

S<  c ue2.  eq  NO )  wr  i  t  o  <  6,  35 )  ibbi.  aa(i-*-2),bb(i  +  l),cc(i).  theta, 
stda.  stdb,  s  t  d  c ,  stdd,  fa,  fb,  fc,  fd 
25  fnrmat(2x,  i4,  x,  3(fB  3,  x  ) ,  9x ,  9x ,  5  <  f  8.  3,  x  ) ,  1 100.  4i2 ) 
c  reset 

cc  <  i  +  1  )=bb  <  i-el  ) 
bb ( i+2)=aa ( i+2) 

if < ibbi  eq  iend.  and  c ue2  eq. NO )  aaend(iplt) =aa < i+2) 
if(i.eq.  7)  go  to  30 
i  =  i  +  l 

3  continue 

c  now  to  do  points  8  to  infinity .  and  beyond  !!. 

C  +•  r  e  •». -f- *■#*#*■  *+*■■**■#■*+*+ 

30  i=B 
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ibbi=ibbi+l 
i i c  t=7 

ifiibbi.  gt.  iend)  go  to  46 
do  10  i i  =  i b  b i ,  i end 
i i c  t=i ic  t+1 

if(cue2.  eq.  NO )  pltct(iplt)=pltct(iplt)+l. 
xn=iab+bc ) /be 

c  note  that  ii  indicates  inhere  gage  c  is. 

call  r ed  ca 1 ( a (  i+2 ) .  b ( i  +  1  ) ,  c  (  i  > .  d(  i— 7  ) .  asc  1 »  dscl»  afetra. 

&  afetrd.  a  z  >  b  z .  cz,  dz,  na .  nb ,  nc,  nd ,  stda»  stdb.  stdc.  s td d.  j  j  ) 
i  f  (  a  z .  1 1.  0.  )  go  to  46 

call  check  (az.  bz.  cz.  dz.  na.  nb.  nc.  nd.  fa.  fb.  fc.  fd) 
if  (fa.  eq.  0.  and.  fb.  eq  0.  and.  fc.  eq.  O)  go  to  28 
if(fa.  eq.  1.  or.  fb.  eq.  1.  or.  fc.  eq.  1)  uiri  te  (6.  104)  ii,  fa,  fb.  fc,  fd 
if  (fa.  eq.  1.  or.  fb.  eq.  1.  or.  fc.  eq.  1)  go  to  27 
104  format(2x,  i4,  x,  /reset2'>  tlOO,  4i3) 
c  Emergency  reset  (for  bad  date  points 

27  ist=ii+l 

i  f  <  c  ue.  eq.  NO )  aaend(iplt)=aa(i  +2 ) 
if  (ist.  gt.  iend)  go  to  46 
i  =  l 

cc ( i >  =  1 . 
bb ( i  +  1 )=1 . 

if!cue2.  eq.  YES.  and.  cue3  eq.  1)  k  ey  =  l 
if(key.eq.  I)  go  to  1 02 
go  to  31 

28  continue 

th  e ta-atan ( ( c  < i ) +c  c ( i ) -b ( i+1 ) -b  b ( i  +  1 ) )/bc) 
c  theta=0. 

aa  (  i  +2)  =c  (  i  )  *c  os  <  th  eta  )+c  c  (  i  )  -a  (  i+2)  *c  os  (theta)-xn* 

&  (c(i)*cos(theta)+cc(i)-b<i+l)*cos(theta)— bb(i+l)) 

l f ( c ue2.  eq.  YES )  call  seal e ( e 1 evl ( i p 1 t ) > e 1 ev2( i p 1 t ) « 
aastrt(iplt).aaend'.iplt).  p 1 1  c  t ( i p 1 1 ) , a  a ( i +2 ) » 

S(  aasca  1  (  i+2  ) ,  i  i c  t ) 

i  f  (  cue2.  eq.  NO.  and  ii.  eq.  iend)  aaend  (  i  p  1 1 )  =aa  (  i+2) 

if  ( f  a.  eq.  0.  and.  fb.  eq  0  and.  fc.  eq.  0.  and.  f d.  eq.  1  )  go  to  84 

29  dc!<i-7)=c(i)*cos(t.heta)  +  cc(i)-d(  i-7)*c  os  (theta)+(cd/bc>* 

?<  (c(i)*cos(  theta)  +c c ( i )-b  < i+1 > *cos  <  theta )-bb  < i+1 >  > 

def 1 ( i-7)-c c  < l -7 ) -d ri ( i -7 ) 
d  e  f  1  c  t  =  d  e  f 1 c  t+ 1 
sumdef=sumdef+defi  (  i  -7 ) 
ssdef=ssdef+(defj  i  i-7,m-*2) 

84  if<cue2.  eq.NO.  and.  ii.eq.  iend)  aaend(  ip  It  )«=aa  (i+2) 

18  continue 

if(cue2.  eq.  YES)  go  to  13 

if  (fa  eq  O  and  fb  eq  O  and  fc.  eq  0.  and.  fd.  eq.  1. 

&  and.  preue.  eq.  YES) 

S>.  iur  ite(6>  36  )  i  i  ,  a  a  (  i  +2  ).  bb(i  +  l)»  cc(i),  dflt,  d  f  1 1 .  theta. 

S'  s  t  d  a  >  stdb.  stdc.  stdd.  fa.  fb.  fc.  fd 

36  forma t(2x.  i4,  x.3<f8  3.  x),2(4x.a5),5(f8.  3.  x).2x.4i2) 
if(fa  eq  O  and  fb  eq  0  and.  fc.  eq.  O.  and  fd.  eq.  1) 
r.  ci  o  Co  13 
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if (cue2.  eq.  YES)  go  to  13 
if  (prcue.  eq  YES) 

Sc  write<6.16)  i  i  .  aa  <  i+2)  .  bb  (  i  +  1  )  .  c  c  (  i  ) .  d  d  (  i-7 )  .  def  1  (  i-7  ) » 

&  theta,  stda.  stdb,  stdc,  stdd.  fa.  fb.  fc.  fd 
16  format  <2x»  i4.  x .  10<f8.  3.  x).  1 100.  4  i2 ) 

13  continue 

if  (cue2.  eq.  NO.  and.  cue3.  eq.  0)  go  to  66 
if(cue2.  eq.  YES.  and.  cue3.  eq.  1)  go  to  66 
aascmn (  iplt)=aminl (  aa  s  cmn ( ip  It). aascal  (  i  +2 )  ) 
aascmx  (  i p 1 1 ) =ama  x  1  (aascmx  (  iplt).  aascal  (  i+2)  ) 
dmin  <  ip  1 1 )  =amm  1  <  dmi n  <  i  p  1 1 ) .  def  1  <  i-7  > ) 
dmax  < i p 1 1 ) =ama  x 1 (dmax ( iplt).  defl ( i-7) ) 

66  if  ( cue3.  eq.  l)call  ppt(aascal(i+2).  ii. aasc  mn  (iplt).  pr  (  ip  1 1 ) .  8. 
&  1.  ) 

i  f  (  c  ue3.  eq.  1 )  call  ppt ( def 1 ( i-7) » i i < dmin < ip  It )• 

&  defdif  (iplt).  9.  0.  ) 

c  reset  cc(l)  -  cc(8)  for  next  iteration 
c  c  ( 1 )  — c  c  ( 2 ) 
cc  <2)=cc (3) 
c  c  <  3 ) =c  c  (  4  ) 
cc  <4)=cc  <  5) 
c  c (5)=cc <  6 ) 
c  c ( 6 ) =■ c c <  7 ) 
cc (7 )=cc (8) 
cc (8)=bb ( i+1 ) 
b  b  < i  +  1 ) =aa  < i +2 ) 


if  < i i .  eq.  iend)  go  to  46 
10  continue 

44  if  (cue2.  eq.  NO)  ur  i  t  e  (  5.  45  > 
if  (  c  ue2.  eq.  NO)  wr  i  te  <  6.  45  ) 

45  f ormat ( //. 3x . ' You  ran  out  of  points  before  7  were  read*!') 
go  to  43 

46  if(cue2.  eq  NO )  u;rite<5.47) 
if  (  c  ue2.  eq  NO)  write(6>47) 

47  format(//3x. 'You  ran  out  of  data  points.  End  of  data  set.  ') 

48  continue 

ifCdeflct.  It.  1)  go  to  89 

stddef =sqrt ( (1.  /  (deflct-1 ) )*(ssdef-( <  sumdef **2 )/<deflct) ) ) ) 
avodef=sumdef/def3ct 


99 

if  (:u<?2 

eq 

NO)  wr i te ( 6, 87 )  sumdef 

.  d  e  f  1  c  t 

if  (cuc-2. 

eq. 

NO)  write <5. 37)  sumdef 

.  def  1  c t 

37 

format ( / 

,  3  x 

.  'The  sum  of  def  1 . 

i  s 

= 

'.  gl6 

S  2x .  i  5.  2x. 

'  1  s 

=  ',  f  10.  3,  2x.  '  std. 

dev.  of 

JJ»17 
i  f  (  c  ue2. 

eq. 

YES.  and.  cue3.  eq.  0) 

go 

to 

87 

i  f  <  c  ue2. 

eq 

YES.  and  cue3  eq  1  ) 

go 

to 

74 

write(5. 

72) 

72 

forma  t ( / 

.  3  x 

, 'Do  you  want  plots? 

/ 

i 

% ) 

accept  73. 

cue2 

T  — 

/  vi 

format ( a 

1  ) 

i  f  <  c  ue2 

eq 

YES  and  cue 3  eq.  0) 

go 

to 

78 

j  f  <  c  u  e  2 

eq 

NO)  go  to  75 

'  C  0 

u>r  ;  \  e  (  5 . 

103)  (  l-l  )  •  .i..rcj!  <  i+2) 

3.  2x.  'd  ivided  by  '. 
defl.  =  flO.  3. /) 
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go  to  74 

102  write<5,103)  ( i i-1 ) .  aascal ( i+2) 

103  f ormat ( //, 3x ,  '***•*  You  hit  a  bad  data  point  while' 

$<  '  writing  on  the  disc.  The  plot  up  to  this  point  has  been' 

&  ,  /. 

&  '  saved.  To  get  the  rest  of  the  plot,  run  CALC  again' 

&  '  starting  at  the  next  point.  The  last  point  plotted' 

S<  ,  /  ,  '  uia  s  '.  i5.  x .  'The  last  scaled' 

&  '  elevation  was  ' .  f 1 0 .  5. / / ) 

74  wr  i  t  e  <  6, 88 ) 
wr  i  te  (  5,  88  1 

9B  format ( //, 33x. '  WHEW  ! ! ! !  Done  with  all  those  plots!  ') 

do  123  i=8. 9 
call  cl ose ( i ) 

123  continue 

87  if(cue3.  eq  1)  go  to  75 

cue3=l 
go  to  85 

75  call  c 1 ose i 7  ) 

c  bozo  printer  scheme  (empty  printer  buffer) 

62  i23=19 

wr i  t  e  ( 6. 59 )  i23 

5*?  format  (  a  1  > 

121=17 

wri te  <6.  59  >  i21 
stop 
end 
C 

SUBROUTINE  REDCALCA,  3,  C.  D.  ASCL,  DSCL.  AFCTRA.  AFCTRD. 
l>.  AZ,  BZ,  CZ.  DZ,  NA.  NB  •  NC.  ND,  STDA.  STDB.  STDC.  STDD,  JJ) 

C  READS  DATA  AND  CALIBRATES  IT.  VALUE=  (DIST.  TO  MEASURING 
C  RANGE  +  REGRESSION  EQ.  )  +  BEAM  ALIGNMENT  CONST. 

INTEGER  TO 

CALL  READC  ( AZ.  BZ.  CZ,  DZ,  NA,  NB,  NC,  ND,  STDA,  STDB,  STDC.  STDD, 

&  ISS,  IENND,  JJ) 

IF<  AZ.  LT.  0)  GO  TO  174 
CONST=.  000032/.  0254 
DX=DZ*CONST 
CX=CZ*CONST 
EX=BZ*-CGNST 
AX*AZ*CONST 
C  FOR  GAGE  203 

STDF=0.  001 26 1 *STDB 
B=13.  695-.  001261*BZ 
C  FOR  GAGE  269 

STDC=0  001261 *STDC 
C  =  1 3  652-0.  00 126 1  *C  Z 
C  FOP  GAGE  270 

STDA=STDA*0  001261 
A=1 3  935-0  00 1 26 1 *AZ  -  AFCTRA 
C  FOR  GAGE  271 

STDI>=STCD*0.  001261 
D=13  903-0  00 1 26 1  *  DZ  - 


AFCTRD 


on  non 


174 
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RETURN 

END 


SUBROUTINE  ADVNCE ( R I  ST ) 

ADVANCES  TO  THE  DESIRED  STARTING  POINT  IN  THE  DATA  SET 
I ST = I F I X ( R I ST ) 

IF<  1ST.  EQ.  1  )  GO  TO  2 
DO  1  1=1,  < I ST-1 ) 

CALL  READC  (A,  B,  C,  D,  NA,  MB,  NIC,  ND,  STDA,  STDE  STDC,  STDD,  ISS, 
S<  I  ND,  J  J ) 

1  CONTINUE 

2  RETURN 
END 


C 


C 

C 

C 

C 

C 

c 

c 

c 

c 


c 


SUBROUTINE  CHECK (AZ,  B2,  CZ,  DZ, NA, NB, NC, ND,  FA, FB, FC,  FD) 
CHECKS  FOR  INVALID  DAT  POINTS 
INTEGER  FA,FB,FC,FD 
FA=0 
FE=0 
FC=0 
FD=0 

IF  (AZ  EG.  0,  OR.  NA.  LT.  30)  FA=! 

IF  ( B  Z.  EG  0  OR.  NB.  LT  30)  FB=1 
IF(C2.EQ  0.  OR.NC.LT.  30)  FC=1 
IF(DZ  EQ  O  OR.ND.LT.  30)  FD=1 
RETURN 
END 

SUBROUTINE  READC  (A,  B,  C,  D,  NA,  NB,  NC,  NC,  SDl.  SD2,  SD3,  SD4,  IS, 

IE,  I  ) 


A.  B.C.D  =  LASERS  A,  B,  C,  D 
NA,  NB,  NC ,  ND  =  NO.  OF  SAMPLES 
SDl , SD2,  SD3,  SD4  =  STANDARD  DEVIATIONS 
IS  =  TIME  TICKS  READ  BEFORE  CLASS 
IE  =  TIME  TICKS  READ  AFTER  CLASS 
SUBROUTINE  READC (A,  B,  C, D, NA, NB, NC, ND, SDl,  SD2, SD3,  SD4, 
I  start,  IEND,  I  ) 

BYTE  NAME (20), ZZ2 

COMMON  /RMECOM/ I EUF (16, 51 2) 

DATA  I OPEN/O/ 

DATA  LR/-1 / 

IF  < I  OPEN  EG  1 )  GO  TO  1 

I0PEN=1 

WRITE  (5,  2) 

FORMAT  (/,  3X.  'RAW  DATA  FILE  NAME  ...  ',*> 

ACCEPT  5.N0C,  (NAME(IZZ),  IZZ  =  1,N0C> 

NAME ( NOC+ 1 ) =0 
FORMAT (G, 20A 1 ) 

CALL  ASSIGN ( 7,  NAME-  MGC  •  'OLD ' > 


noon 
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DEFINE  FILE  7  <  1000,  8192,  l>,  IV) 

WRITE (6. 4)  ( NAME < ITT ) ,  ITT=1,N0C) 

4  FORMAT  (  /,  3X,  'DATA  FILE  NAME:  '.80A1) 

WR  I  TEC  5,  3) 

3  FORMAT (  /,  3X,  'ROARRRRR !  !  !  CHOMP  CHOMP  CHOMP') 

1  =  17 

1  IR=( ( 1-1 )/512)+l 

IF <  IR  EG.  LR  )  GO  TO  30 

READ ( 7 ' IR, END=17 17 )  ( ( IBUF(MZ,  LZ  >, MZ  =  i,  16),  LZ  =  1, 512) 

LR=IR 

30  INDX=MOD( I,  512) 

IF( INDX.  EG. O)  INDX=512 
A=IBUF< 1, INDX) 

D= IBUF < 2, INDX) 

C=IBUF ( 3, INDX) 

D=IBUF(4,  INDX) 

NA=  I BUF  (  5,  INDX) 

NB=IBUF<6,  INDX) 

NC=IBUF ( 7,  INDX) 

ND= I BUF ( 8, INDX) 

SDl = IBUF ( 9, INDX) 

SD2= I BUF  <10,  INDX) 

SD3= IBUF ( 1 1, INDX) 

SD4=ISUF( 12, INDX ) 

ISTART=IBUF (13. INDX) 

I END* I BUF ( 1 4, INDX) 

1  =  1+1 

IF  (A.  NE.  -1)  RETURN 
1717  A=— 1 

RETURN 

END 


HERE  ARE  THE  PLOT  ROUTINES 

SUBROUTINE  AXIS(PMAX,  PMIN,  IPLT, PR,  IU,  ID,  ELEV1.  ELEV2.  IX) 
BYTE  ID  <  1  > 

EELEV1=ELEV1 /12. 

EELEV2=ELEV2/12. 

RANGE=PMAX— PMIN 
PR  =  1 . 

IF  <  RANGE.  GT.  .  005)  PR  =  .  01 
IF  (RANGE  GT.  .01)  PR=.  02 
IFIRANGE.  GT.  .  02)  PR  =  .  05 
IF  (  RANGE.  GT  .  05)  PR  =  .  1 
IF  (RANGE  GT  .  1  )  PR=.  5 
I  F  (  RANGE.  GT.  .  5)  PR=1 
IF  (  RANGE.  GT.  1 )  PR=2. 

IF  ( RANGE  GT.  2)  PR  =  5. 

IF(RANGE.  GT.  5)  PR=10. 

IF(RANGE.  GT  10)  PR=50. 

IF  <  RANGE  GT  50)  PR=100. 

IF < RANGE.  GT  100)  PR  =  250 
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IF  (RANGE.  GT.  250)  PR  =  500. 

IF < RANGE.  GT.  500)  PR=70Q. 

IF  (RANGE.  GT.  700)  PR=1000. 

IF (RANGE  GT. 1000)  PR=1500. 

IF(RANGE.  GT.  1500)  PR  =  2000. 

IF  (RANGE.  GT.  2000)  PR=2500. 

P  P  M 1 N=P M I N / 1 2 . 

PPR=PR  / 12. 

IF( IX  EQ.  1 ) WRITE (8,  1 ) ( ID( IV),  IY=1, 25),  IPLT, EELEV1 , 

EELEV2, PPR, PPMIN, (PPP+PPMIN) 

1  FORMAT ( 1  HI ,  ///, 60X , 25A1 , 2X,  '*  ' ,  14,  /, 

&  5X,  'START  ELEV.  =  ',F10  3,  '  FT.  ',/,5X,  'END  ELEV.  =  \ 

?<  F10.  3,  '  FT.  ', 

&  70X,  'RANGE  =  ',  F8.  3,  X-  'FT.  ',  //, 

12X,  '  MIN.  ' ,  F8.  3,  S3X,  '  MAX.  ',FB.  3,//) 

IF  (IX.  EQ.  2) WRITE (9,  4)  (  ID(IY).  IY*1,  25),  IPLT,  PR,  PMIN.  (PR+PM1N) 

4  FORMAT  ( 1H1 ,  ///,  60X,  25A1.3X,  'fc',14,  /  ,  T102,  'RANGE  =  ', 

%>.  F8.  3,  X,  'IN.  ',  //,  14X,  'MIN.  ',  F8.  3,  '  IN  ', 

78X,  'MAX.  FS.  3,  X,  'IN  ',  //) 

WRITE! IU, 3) 

3  FORMAT  ( 8X,  'VALUE',  BX,  'DIST  ',/, 

4X,  'IN  ' ,  7X,  'FT.  ' ,  T23,  'FT.  ', 

•5-  T27 .  '  :  10(9X,  '  !  '  )  ) 

WRITE< IU, 2) 

2  FORMAT (T27, 10< '0123456789'),  '0') 

RETURN 

END 

C 

SUBROUTINE  PPT(C,  1 1 ,  XMIN, PR,  IU,  X ) 

BYTE  PCHAR( 100), STAR,  PD,  PP,  PE DT 
DATA  STAR/1H*/ 

factor  =  ioo.  /pr 
c 1 2=c / 1 2 

CC=C— XMIN 
NUM=FACTOR*CC-*-l 
IF(NUM.  GT.  100)  NUM= 1 00 
IF ( NUM. LT  1)  NUM= 1 
DO  1  M=l, 100 
1  PCHAR ( M ) =STAR 

WR I TE ( ! U, 2 )  C , C 1 2,  I  I ,  C “CHAR ( N) ,  N= 1 ,  NUM  > 

2  FORMAT  (  X  .  2  '  F9  3,  X),  14,  1 02A 1  ) 

RETURN 

END 

C  THIS  IS  THE  SCALE  ROUTINE  (SCALES 
C  THE  AA'S  TO  USER  IN“UT  ENDPOINTS) 

SUBROUTINE  SCALE  (ELEV1 ,  ELEV2,  AASTRT,  AAEND,  PLTCT,  AA, 
y  SCALAR, N ) 

FN^FLOAT ( N ) 

SC  ALAA=AA-*-  (  FN—  1  )*(  (  (  ELEV2-AAEND  )-<  ELEV1 -AASTRT )  >/(PLTCT-l  )  ) 

£,  t- ( ELEV  1 -AASTRT) 

RETURN 

END 
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Program  ADJUST,  used  to  adjust  the  endpoints  of  the  profile. 


35736. ouh . 

pfiles(get. hngrsur) 
pfiles<get.  hngr02c> 
mnf ( t ) 

pfiles<put. compare) 
fceor 

program  ad  just (input,  output. hngrsur.  hngr02c»  compare. 

Jk  tape5=input.  tap e6=output,  tape7=hngrsur,  tape8=hngr02c. 
&  tap  e9=c  ompar  e  ) 

real  las(600),  laser,  1  a ss tr t , newsur ( 600 ) ,  lasend.  lasct. 

?<  neuil  ( 600 ) 


c  converts  the  laser  file  to  feet,  and  adjusts  both  to  70' 
c  endpoints 
c 

elevi=70. 
e  1  ev2=70. 
is=0 
i  1=0 

c  survey  pndpts 

do  1  ii  =  l,2222 
i s=i s+1 


33 

2 

1 

3 

c  laser 


survey=z survey 

r eati ( 7. 33. end=3 )  i  »  zsurvey 
format  (  i  5.  f  10.  5  ) 
if(ii.  eq.  1)  surst=z survey 
format(i5.  g20.  5) 
continue 
surend=survey 
endp  oints 
do  4  i i i = 1  2222 
i l=i 1+1 


laser  =  z 1 aser 

read ( B, 2.  end=5)  i.zlaser 
iffiii.  eq.  1)  lasstrt=z  laser 

4  continue 

5  continue 
lasend=laser 
reii'ind  7 
remind  8 

c  adjust  and  store 

lasct=float< i 1  —  1  ) 
surct=float< is— 1 ) 


llasct=ifix ( lasct) 


isurct=ifix(surct) 

do  6  i=l»l lasct 

read < 8, 2. end=7 )  ii,.  laser 

laser=laser /12. 

x l=e 1 evl - ( 1 asstrt/ 12.  ) 

x2=r  1  ev2-  (  1  zisend  /  12  ) 


i 
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grad~ ( x2-x 1 ) / ( 1 asc t-1 .  > 
neu»l  (  i  )  =  laser+grad*(  i  —  1  )  +  xl 

6  continue 

7  do  8  i=l,  i surct 

read <7, 2, end=9 >  iii» survey 
x l=elevl— surst 
x2=e 1 ev2-s  ur end 
grad=(  x2-x  1  >  /  <  sure  t—  1 .  ) 
neuisur  (  i  )  =survey +gra d*  (  i- !  )  +  xl 
3  continue 

9  continue 

write(6i 12)  lasstrti lasend, surst, surend, isurct. llasct 
12  format(//,  5x,  'starts  end=  ' ,  4<  f  lO.  5,  x  )  >  x>  2i  10) 

ma x=ms xO ( i sur c t , llasct) 
do  10  ii=l, max 
d  i  f=neu)sur  (iz>-neu)l<iz> 

urr  i  te  <  6,  11)  i?,neu)sur(iz),r,eu)l(iz)>dif 
a>rite(9,  11)  iz,netitsur(i7>,newl(iz),dif 
11  formate  i  5.  3(f  10.  5,  x)  ) 

10  continue 
stop 

end 

#eor 
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APPENDIX  D 

System  Usage 


1.  Technical 

This  section  explains  briefly  how  to  set  op  and  how  to 
use  the  system. 

The  following  steps  describe  how  to  set  up  the  system  - 

1.  Mount  the  laser  gages  in  the  rigid  beam  so  that 
all  the  laser  beams  lie  in  a  plane*  and  are 
paral lei* 

2.  attach  the  rigid  beam  to  the  load  vehicle  with 
the  spring  loaded  supports* 

3.  attach  the  fifth  wheel  to  the  load  vehicle* 

A.  make  all  electrical  connections*  and 


4.  drive  over  the  test  pavement. 


Better  equipment  operation  will  result  if  - 


1.  "The  operator  reads  all  the  manuals  on  all  parts 
of  the  equipment,  and  observes  the  manufacturers 
directions". 

2.  the  disk  drive  is  not  subjected  to  shock.loads 
(this  may  cause  the  format  to  change,  making  all 
previously  made  disks  unreadable),  and 

3.  the  equipment  is  at  a  steady  state  temperature 
before  operation. 

2.  Philosophical 

This  system  was  intended  to  measure  longitudinal  pro¬ 
file.  deflection  at  a  point,  and  texture.  If  the  rigid  beam 
was  removed  from  the  load  vehicle  and  cantilevered  normal  to 
the  direction  of  travel,  a  definition  of  the  lateral  deflec¬ 
tion  basin  can  be  obtained  when  the  load  vehicle  is  driven 
past  the  beam. 

If  the  system  is  run  laterally  acr oss  a  pavement,  the 
profile  will  reflect  the  rut  depth. 
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APPENDIX  E 


Gage  Alignment  Scheme 


In  order  for  the  algorithm  described  in  Section  III,  B, 
1  to  work,  one  of  the  following  must  hold  - 


1.  The  gages  must  not  move  relative  to  one  another, 
and  their  relative  positions  must  be  known,  or 


2.  the  gages  can  move  relative  to  one  another,  and 
the  movement  is  known. 


The  first  condition  was  used  in  this  study.  The  following 
scheme  was  used  to  get  the  relative  positions  of  the  gages. 


Referring  to  Figure  III  —  1 ,  let  *'i"  be  one  step  in  the 
direction  of  motion.  Note  that  - 


AA .  ■  BB.  _  =  CC.  _ 

i  l+l  1+2 


Also  note 


.n— ^  _ , 
i  =  l  l 


it  * 

i  =<=:  i 


„n  __ 

2  CC  . 

i=J  l 
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For  large  numbers  of  readings  (large  n)»  equation  (2)  can  be 
approximated  by 

2*  .  AA.  =  2^  ,  BB.  =  2n_1  CC  C3> 

1=1  i  i=l  i  i=l  i 

If  n  is  not  large  enough,  an  error  may  result.  This  is  exam¬ 
ined  later. 

Note  that,  at  any  time,  the  endpoints  of  A  +  AA.  B  -r 
BB.  and  C  +  CC  define  tu>o  straight  lines  (the  datum  and  the 
rigid  beam).  Therefore,  the  endpoints  of 


„n 


S'.'  .  (A.  +  Aft.,  .  211  .  (E.  +  BB.S  and  2"  ,  <C, 
i=l  l  i>  1=1  l  l >  i=l  i 


^T' 


<C.  +  CC . .  <4> 


i> 


also  form  two  straight  lines. 


Because  of  (3),  it  follows  that 


•"  ,  <  A. 
i  =  l  l 


Ah  > 

l 


,.n 


.  .  < AA . > 

1  =  1  l 


C5) 


„n 


i=l  <  *x 


BE  > 

l 


„Tl 


,  <  BB  .  > 
i  =  l  i 


<6> 


2  ,  <  C  +  CC 

1  =  1  i  i 


2r'  ,  <CC  > 

i*l  l 


<7> 


form  a  straight  line. 
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Rewriting  (4),  <5)>  and  (6)  and  dividing  bg  n  gives 


n 


n 


n 


<8> 


whose  endpoints  also  form  straight  lines. 

When  a  test  is  run*  the  beam  alignment  constant  can  be 
determined  if  the  averages  at  (8)  are  known. 

The  endpoints  of  the  averages  should  represent  two 
straight  lines.  If  they  do  not.  the  difference  between  one 
average  and  the  straight  lines  formed  by  the  other  averages 
is  the  amount  that  that  gage  is  out  of  line  with  the  other 
two  gages.  That  difference  (the  out-of-line  constant) 
should  be  added  to  the  calibration  equation  for  that  gage 
(the  effect  is  to  remove  a  constant  offset  from  that  gage). 

This  method  can  be  extended  to  the  D  gage,  if  the  pave¬ 
ment  is  not  deflected.  In  the  testing  program,  a  lightweight 
tricycle  was  used  for  this  purpose.  It  carried  the  beam 
without  significantly  loading  the  pavement.  The  tricycle  was 
pulled  by  hand. 

The  out-of-line  constants  are  determined  automati cally 
by  the  computer  program  CALC 
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As  mentioned  previously!  equation  (3)  is  not  strictly 
correct.  While  equation  <3)  is  true  for  large  n#  it  is  ren¬ 
dered  invalid  by  special  conditions  at  the  start  and  end  of 
the  test.  An  examination  of  Figure  III — 1  shows  that  A  is 
never  where  C  is  rt  the  start.  Neither  is  C  ever  where  A  is 
at  the  end.  Thus!  a  more  proper  representation  of  equation 
(3)  is 


AA. 

1 


<2> 


This  leaves  the  beginning  and  ending  terms  out.  Those  terms 
are 


CC. 


CC2 


BB, 


BE 


Tl 


Ah 


n-1 


and  AA 


n 


<9> 


If  the  out-of-line  constant  is  calculated  using  equa¬ 
tion  (2)  instead  of  equation  (3)(  no  error  is  introduced. 
The  error  introduced  by  using  equation  (2)  is 


where 


error 


-  <.2B9  -  A9  -  C9> 

Tt 


<10> 


error  =  error  in  the  gage  A  out-of-line  constant 

A9  =  AA  ,  +  AA 
n-1  n 

B9  «  BB,  +  BB 
1  n 


C9  =  CCj  +  CC2 
n  *  number  of  iterations 
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If  the  sum  of  the  AA  terms  (A9)  equals  the  sum  of  the  BB 
terms  (B9)i  and  either  equal  the  sum  of  the  CC  terms  (C9)  , 
no  error  occurs  in  the  calculation  of  the  out-of-line  con¬ 
stant.  .  If  one  of  the  sums  just  mentioned  is  different  from 
the  other  tuo>  an  error  will  occur.  This  error  is  a  function 
of  how  great  the  difference  is.  and  how  large  n  is.  If  the 
d if f erence ( s )  is  small,  and  n  is  large,  the  error  will  be 
small . 


* 

! 


s 

t 
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APPENDIX  F 

List  of  Equipment  Suppliers 


The  following  firms  either  manufactured  or  supplied 
equipment  for  this  research.  Only  major  equipment  is  listed. 

1.  The  computer.  The  computer  was  purchased  from 
Heath  Corporation,  Benton  Harbor,  Michigan.  All 
the  peripherals,  except  the  bus  extender,  were 
also  purchased  from  Heath. 

2.  The  bus  extender  The  bus  extender,  manufactured 
by  Digital  Equipment  Corporation  (Maynard.  Mas¬ 
sachusetts),  was  purchased  from  Hami 1  ton— Avnet 
(Culver  City,  California). 

3.  The  laser  gages  The  laser  gages,  called  Optoca- 
tors,  were  manufactured  and  sold  by  Selective 
Electronic,  Valdese,  North  Carolina. 

4.  The  fifth  wheel  The  fifth  wheel  was  manufac¬ 
tured  and  sold  by  Labeco,  Mooresville,  Indiana. 
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The  following  equipment  appeared  in  the  section  VI* 
Suggestions  for  Further  Research. 

1.  A  pr inter— p 1  otter.  The  Printronix  printei — 
plotter  is  manufactured  and  sold  by  Printronix. 
Irvine.  California. 

2.  Sonic  gages.  The  Polaroid  sonic  gages  are 
manufactured  and  sold  by  Polaroid  Corporation, 
Cambridge,  Massachusetts. 
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APPENDIX  G 

Data  File  Previewer 


This  appendix  describes  what  the  data  file  previewing 
program  LIST  does,  and  why  it  is  useful  to  use  it. 

LIST  converts  the  raw  binary  data  to  ASCII  format  data. 
At  the  users  request,  it  will  list  the  data  on  either  the 
terminal,  the  printer  or  the  disc. 

The  raw  data  consists  of  five  parts  - 

1.  The  station  number,  in  feet. 

2.  the  average  of  the  n  readings  taken  by  each  gage 
over  the  four  inch  distance. 

3.  the  number  of  readings,  n.  taken  in  by  each  gage 
over  the  four  inch  distance. 

4.  the  standard  deviation  of  the  n  readings  for 
each  gage,  and 

5.  the  time  the  readings  started  and  ended,  accu¬ 


rate  to  one-sixtieth  of  a  second. 
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The  advantage  of  looking  at  the  ram  data  before  cali¬ 
brating  and  reducing  it  with  CALC>  is  that  it  can  be  done 
quickly.  The  operator  can  tell  almost  immediately  if  the 
pavement  needs  to  be  tested  again.  By  examining  the  uncali— 
brated  readings*  the  operator  can  see  where  the  gages  were 
out  of  range  for  the  entire  four  inches  (the  averages  are 
zero).  If  there  are  too  many  of  these,  the  operator  may  wish 
to  run  the  test  again.  Zeroes  may  indicate  areas  that 
require  closer  examination  (potholes,  etc.  ). 

The  number  of  readings  is  useful.  If  the  number  of 
readings  is  small,  it  could  mean  that  a  pothole  was  driven 
over.  If  the  number  of  readings  is  different  for  each  gage, 
and  the  average  readings  are  similiar,  it  means  that  the 
pavement  has  small  "holes"  in  it,  where  drop-out  readings 
occur.  Porous  friction  surface  pavements  display  this 
behavior. 

If  the  averages,  the  numbers  of  readings,  and  standard 
deviations  all  compare  well  with  their  respective  counter¬ 
parts,  the  pavement  is  a  typical  one  -  reeltively  flat,  with 
no  "holes". 

This  system  has  no  outward  indicators  that  it  is  work¬ 
ing  properly.  By  using  LIST,  the  operator  can  find  out 
quickly  if  the  system  *'as  failed. 

Table  iy— 1  is  some  raw  data  from  a  porous  friction  sur- 

f  jCC 


